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Effektive Wirkungen fur stark wechselwirkende 
fermionische Systeme 

Zusammenfassung 

Wir vergleichen vcrschicdcnc nicht-storungsthcoretische Methoden zur Beschreibung fermionischcr 
Systeme, die gebundene bosonische Zustande (BBS) und spontane Symmetriebrechung (SSB) 
aufweisen. In einer rein fermionischen Sprachc erfordert das Eindringen in die SSB Phase Tcch- 
niken jenseits von Storungstheorie und Renormierungsgruppengleichungen. Dazu ist eine Beschrei- 
bung, die BBS und elementare Teilchen glcichbcrechtigt bchandclt, besser geeignet. Die "Particlle 
Bosonisierung" fuhrt aber zu einer Willkiir in der Wahl der BBS Folder, da diese durch die klassische 
Wirkung nicht eindeutig festgelegt ist. Die Ergebnisse approximativer Rechnungcn, z.B. mean field 
theory konnen aber von dieser Wahl abhangen. Dies beschrankt die quantitative Aussagekraft. 
Am Beispiel des Nambu-Jona-Lasinio-Modclls zeigen wir, wie diese Abhangigkeit durch geeignet 
gewahlte Approximationen reduziert und manchmal sogar zum Verschwindcn gcbracht wcrdcn 
kann. 

Schwinger-Dyson-Gleichungen (SDE) erlauben eine Beschreibung von SSB ohne Hilfsfelder. Die 
2PI- Wirkung ermoglicht es uns, verschiedene Losungcn der SDE zu vergleichen und so die stabile zu 
linden. Diese Methode wenden wir auf eine sechs-Fermion Wechsclwirkung an, die der drei-Flavor- 
Instantonwechselwirkung der QCD ahnelt. Wir linden einen Phascniibergang erster Ordnung in 
die chiral gebrochene Phase, aber keine stabile Phase mit gcbrochener color-Symmetrie. 

Die Existenz eines elementaren skalaren bosonischen Teilchens im Standardmodcll - dem Higgs 
- fuhrt zu mehreren Fragcn. Die mit fundamcntalen Skalen (~ Mqut) verglichen kleine Masse 
erfordert ein extremes Mafi an Finetuning. AuBerdem ist das </> 4 -Potential moglicherweise nicht 
renormierbar im strengen Sinne. Im Hinblick darauf diskutieren wir die Moglichkeit, dafi das Higgs 
ein BBS aus Fcrmionen ist. 



Effective Actions for Strongly Interacting Fermionic 

Systems 

Abstract 

We compare different non-perturbativc methods for calculating the effective action for fermionic 
systems featuring bosonic bound states (BBS) and spontaneous symmetry breaking (SSB). In a 
purely fermionic language proceeding into the SSB phase requires techniques beyond perturbation 
theory and renormalization group equations. Improvement comes from a description with BBS 
fields and elementary fields treated on equal footing. Yet, "partial bosonization" introduces an 
arbitrariness as the choice for the composite fields is usually not completely determined by the 
classical action. Results of approximate calculations, e.g. mean field theory, may depend strongly 
on this choice, thus limiting their quantitative reliability. Using the Nambu-Jona-Lasinio model as 
an example we demonstrate how this dependence can be reduced, sometimes even be eliminated 
by suitably chosen approximations. 

Schwinger-Dyson equations (SDE) allow for a description of SSB without auxiliary fields. The 2PI 
effective action enables us to compare different solutions of the SDE and find the stable one. We 
apply this method to a six-fcrmion interaction resembling the three-flavor instanton interaction in 
QCD. We find a first order chiral phase transition but no stable phase with broken color symmetry. 

The existence of an elementary scalar boson in the Standard Model the Higgs - raises several 
questions. The smallness of its mass compared to some fundamental scale (~ Mgut) requires an 
extreme amount of fine-tuning. Moreover, its 4 -potential may not be renormalizable in a strict 
sense. In view of this we discuss the possibility of a Higgs as BBS of fermions. 
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Chapter 1 
Introduction 



There are no eternal facts, as there are no absolute truths. 

Friedrich Nietzsche 

1.1 Bosons Made up of Fermions 

Everything is Made of Fermions? 

The search for "fundamental" particles is one of the most ambitious enterprizes 
in physical research. While it has been very successful in terms of uncovering new 
particles it has also led us to question over and over again what "fundamental" really 
means. The notion "fundamental particle" has changed with time. First it was atoms 
then it was electrons and nuclei and later the latter ones were split into protons and 
neutrons. Today, we are quite sure that even neutrons and protons are made up of 
the more fundamental quarks. What once had been fundamental particles became 
bound states. 

Looking at Fig. 11.11 we can see that this evolution is simply like turning up the 
resolution of a microscope. At a low resolution we see nothing but a point (particle) 
while at a higher resolution it exhibits structure, i.e. we can see that it is composed of 
other particles. In particle physics the "microscope" is a scattering experiment and 
the resolution improves with smaller wavelength A ~ ^ , and therefore higher energy 
E of the scattering particle. Consequently, with progress in accelerator physics, it 
might turn out that some (or even all) of the particles of the Standard Model are 
not fundamental, but bound states. 

What has all that to do with fermions? At the moment we are in a very peculiar 
situation. Classifying particles by spin (S) and statistics all particles of the "Stan- 
dard Model" are either vector bosons with S = 1 or fermions with S — |. The only 
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Figure 1.1: With the exploration of smaller and smaller scales it often became ap- 
parent that particles which were thought to be fundamental are instead composed 
of even smaller particles. Nevertheless, for a description at a given scale it is often 
useful to treat bound states as "fundamental" particles, e.g. for a first description of 
water vapor it is a very good approximation to treat the water molecules as funda- 
mental, i.e. we use an "effective theory" in which the water molecules are pointlike 
and have some kind of effective interaction. However, if we want to describe the 
absorption and emission of electromagnetic waves it becomes necessary to consider 
the water molecule as composed of atoms (infrared) and eventually the latter ones to 
be composed of electrons and a nucleus (visible, ultraviolet). Corresponding to the 
scale we probe, we have to consider different effective theories. This kind of scale 
dependent description is one of the main ideas behind the renormalization group 
(RG) (cf. Sect. 12.1. "2^) . The RG provides us with a means to calculate the couplings 
of an effective theory at a given scale from the ones of an effective theory valid at a 
smaller length scale, and ultimately from the fundamental one. 



exception is the Higgs as it is supposed to be a scalar (S = 0) boson. Although, 
recent experiments |l|; 0; 0; 0; 0| at the final runs of LEP may have detected a Stan- 
dard Model Higgs (or a supersymmetric extension of it), the evidence is still quite 
shaky, leaving room for speculations. Being provocative we could state that so far we 
have not yet observed any fundamental scalar boson. Hence, one could conjecture 
that the Higgs might be a bound state of fermions, more explicitly a top-antitop 
bound state p; S B HI Ell ■ I n addition to this phenomenological aspect, the idea of 
"top-quark condensation" might also be a way to circumvent some technical prob- 
lems, like the triviality of 4 -theory (Higgs potential!) 
the way to a "renormalizable" Standard Model. 



11; 12; 13; 14; 15 , paving 



Spontaneous Symmetry Breaking 



The phenomenon of spontaneous symmetry breaking (SSB) and the formation of 
bosonic bound states in strongly interacting fermionic systems are tightly connected. 



1.1. Bosons Made up of Fermions 
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First of all, according to Goldstone's theorem |ly; [Tjj it is inevitable to have 
massless Spin-0 particles, i.e. scalar bosons, if we spontaneously break a continuous 
global symmetry. In a purely fermionic theory these must be bound states. An exam- 
ple for these bosons are the pseudoscalar mesons of QCD, which are the Goldstone 
bosons of chiral symmetry breaking. 

Secondly, it is impossible for a fermionic (Grassmannian) field to acquire a non- 
vanishing vacuum expectation value. Therefore, the simplest possible symmetry- 
breaking term is a bosonic operator made up of two fermions. Thus, even for purely 
fermionic theories, the phenomenon of SSB is characterized by bosons. 

Physically, we can imagine a Mexican-hat-type potential for the bosonic com- 
posite operator. Depending on whether we have a local or a global symmetry the 
Goldstone bosons, corresponding to the angular excitations, may or may not be 
eaten up by gauge bosons. But, in addition to those we always have the radial ex- 
citations, corresponding to (usually) massive bosons. E.g. in the top condensation 
model this would be the Higgs boson. So, at this level, there is really no difference 
to a model with elementary bosons. 

Yet, there is a slight difference in the way the Mexican-hat is generated. For 
elementary bosons we quite often simply choose the potential to be a Mexican-hat, 
whereas for fermions the generic case is that the non-trivial minima are generated 
dynamically by quantum fluctuations. That is why it is often referred to as dynamical 
symmetry breaking. 

Finally, we note that the case where the described radial boson becomes massless 
corresponds to a second order phase transition. With features like universality these 
special cases are especially interesting. 

Last but not least, of course, there can be bound states not directly linked to 
SSB, e.g. the hydrogen atom or positronium. 



A lot of Bound States — Some Models 

Of course, a speculative model of the Higgs and chiral symmetry breaking are not 
the only situations where we encounter the mechanisms described above. 



Color superconductivity: At very high density it is expected that the 




.ensity it is expected tn; 
ground state is a color superconducting phase [H 0; 0; Q 0; H I2J; 
Depending on the specifics of temperature, number of flavors etc. there are several 
different phases. Let us just mention the color flavor locking phase characterized by 
a non- vanishing expectation value of 1 (V'liV'lj) ~ (^ia^bj — ^-^ij^ab) as an example. In 
this phase the S77(3) co i or x SU (3)l x SU (3)r is broken down to a vectorlike SU (3)y . 

: L refers to left handed, R to right handed. 
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Breaking global (chiral) as well as local (color) symmetries, we have Goldstone 
bosons as well as massive gauge bosons (gluons). 



Color symmetry breaking in the vacuum: It was conjectured |27| |28| that 
the QCD ground state at zero density might also have a broken color symmetry. 
In this case the order parameter is (V'LiV'Rj) with the same type of expectation 
value and symmetry breaking pattern as for the color flavor locking phase in color 
superconductivity. It is worth noting that the massive gauge bosons can then be 
associated to the vector mesons of QCD. We will investigate this possibility briefly 
in Chap. |H1 

Chiral symmetry breaking: This is probably one of the most studied cases of 
SSB |29j |. A vacuum expectation value (ipu^Ri) ~ $ab gives a mass to the (nearly) 
massless quarks. The pseudoscalar mesons (pions etc.) bear physical witness of this 
process. 



Superconductivity: Yeah, just plain old ordinary superconductivity |30j] in 
ordinary matter like metals at temperatures of some Kelvin. It comes about due to 
the formation of Cooper pairs and the condensation a non-vanishing (ipip) where ip 
is now an ordinary electron field. This gives us an example where we are at a scale of 
some meV instead of several hundred MeV (chiral SSB) or even GeV (Higgs model). 



1.2 Describing Bound States 



As there are plenty of systems featuring bosonic bound states we better start looking 
for ways to calculate something useful. Interesting quantities are, of course, masses 
and couplings of the bosons, and a potential vacuum expectation value. 

A very useful tool to study such quantities is the effective action 12; EH Elf. 
replacing the action of classical field theory, it allows for a simple description of SSB, 
yet it includes quantum effects. Hence, this is the object we would like to calculate. 



The NJL model 



Many of the problems associated with the description of bound states can be studied 
already in a very simple NJL 2 -type model (for only one fermion species) with a 

2 NJL stands for Nambu and Jona-Lasinio who used this model to study chiral symmetry break- 
ing |29| . Due to its simplicity this and similar models are still very popular, e.g . a sizable part of 
the studies on color superconductivity is based on this model 0; |3 |2{1 El |22| 13 HH ■ 
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chirally invariant pointlike four-fermion interaction: 



d 4 x | ifrifrf) (1.1) 



Depending on the value of A = (A CT , Ay, Aa) we are in a symmetric phase or in 
a phase with broken chiral symmetry. As it turns out, the critical A cr i t separating 
these two phases is an interesting but relatively easy to calculate quantity. 



Arbitrary Parameters — Fierz Ambiguity 

The simplest calculation which comes into ones mind is probably a mean field cal- 
culation. We will apply this method (well-known from statistical physics) to the 
model Eq. (jl.l|) at the beginning of Chap. El We find that there is a basic ambiguity 
connected to the possibility to perform Fierz transformations (FT) on the initial La- 
grangian - we will refer to it as Fierz ambiguity. This Fierz ambiguity can influence 
the value of the critical coupling quite dramatically, severly limiting the applicability 
of MFT " 



The origin of the Fierz ambiguity can be understood quite well when looking at 
the model of Eq. (jl.lj) . Due to the Fierz identity (s. App. El for our conventions on 
7-matrices) 



[{^f _ (^ 7 fy) 2 ] + 2 [{^f - ty^W] = (1.2) 
only two of the quartic couplings are independent and we write 

A CT = K + 27V, Ay = (1 - 7) Ay, A A = 7V- (1.3) 

Where 7 parametrizes the "symmetry" , and the A are "invariant" couplings. 

In a naive way one would like to combine the fermions in the four-fermion inter- 
action of Eq. (jl.lj) into pairs and interpret those as bosons, e.g. one would like to 
take the term multiplying Xa pair it into two i/j^^^i/j ~ A 11 and hence interpret this 
term as an interaction (mass term) for "axial vector bosons". However, using Eq. 
(jl.2|) we can now transform this term to zero, eliminating the "axial vector bosons" . 
As this pairing is the basic idea of MFT, an ambiguity seems inevitable. 
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Auxiliary Fields 

There are several standard methods (perturbation theory, RG equations, Schwinger- 
Dyson equations (SDE)) which allow us to calculate the critical coupling without any 
reference to a pairing into bosons. The results are then naturally unambiguous. Yet, 
these descriptions have their problems, too. Perturbation theory is unable to describe 
SSB at all, the RG calculation cannot be extended into the broken phase without 
considerable calculational difficulty and SDE's beyond the simplest approximation 
become quite difficult to solve, too. At least part of the problem is that we lack 
an intuitive understanding of the momentum dependence of complicated fermionic 
operators. 

Looking at the example of mesons again it is known that an effective theory 
with bosonic meson fields interacting with themselves and the quarks works quite 
well without having very complicated terms in the effective Lagrangian. In addition, 
bosonic fields allow for a good understanding of SSB. Thus, it seems a reasonable 
step to introduce auxiliary fields to describe the bound states. Formalizing the naive 



pairing procedure, partial bosonization |35|; |36f |37|; |38|; |39| (cf. Chap. 0J) leads to a 
model with massive bosonic fields and Yukawa-type interactions, but no four-fermion 
interactions, 

S B = Jd 4 x | i^ + ^J^+^v^ + l^A^ (1.4) 



The identification 

h 2 h 2 h 2 

2 'V 2 U V 2 11 A /i r\ 

/V = My = f i A = Y~ ( L5 ) 

lA a Ay Aa 

makes this model equivalent to the NJL-type model 

However, due to the Fierz identity ()1.2|) the couplings in (jl.4|) are not unique, 
bringing back the ambiguity of MFT. Indeed, MFT appears as a simple approxima- 
tion to this model, neglecting all bosonic fluctuations. 

This is not the only situation in physics where an arbitrary parameter (in our 
case 7) appears in calculations. Prominent examples are the renormalization point fi 
or the gauge fixing a in gauge theories. Ultimately, in an exact calculation physical 
quantities should not depend on such a parameter. Nevertheless, approximate cal- 
culations usually do. Improvement in the approximation often tends to reduce the 
dependence on such parameters. Sometimes, there are even special approximation 
schemes where we can achieve independence of such parameters, e.g. every order of 
perturbation theory provides such a scheme for the gauge fixing parameter. Beside 
the practical advantages, finding such an approximation scheme also shows that 



1.2. Describing Bound States 



7 



there is nothing fundamentally wrong with the method in question. From another 
point of view, considering an approximation which is not independent of the ar- 
bitrary parameter, we can say that the spread of the results under a variation of 
the arbitrary parameter gives us an estimate of the minimal uncertainty of a given 
approximation. 

We will find a very sim ple example of an approximation which is independent 
of the Fierz parameter jifli l4l| . Interestingly, it turns out that this introduces the 
concept of scale dependent degrees of freedom. Technically we use the scale depen- 
dence of our auxiliary fields to keep the form of the (effective action) simple [42]. 
However, on a deeper level we would like to interpret this simplicity as a first step 
to a description with the "right" degrees of freedom at every scale. 

Unfortunately, there is a huge number of possible field redefinitions and in simple 
approximations, e.g. so called local potential approximation (LPA) |43|; |4J; |45( , it is 
a priori not clear which is the "correct" one. A criterium can be obtained only by 
the consideration of terms with derivatives of the fields. In the end this leaves us to 
choose between high algebraic (and/or numerical) difficulty and a "physical guess". 



2PI Effective Action 



In view of the many problems connected with auxiliary fields it seems prudent to 
look for alternatives. One possibility is the 2PI effective action 0;II3;I3. After 



the introduction of sources for the composite operators bilinear in the fields we 
can omit the introduction of auxiliary fields and directly perform an additional 
Legendre transformation with respect to the sources of the composite operators. 
For purely fermionic theories it turns out that this description is redundant and 
the dependence on the propagators is sufficient. Hence, we can omit the dependence 
on the fields, leaving us with a description completely in terms of bosonic variables 
(the propagators). Nevertheless, this method is naturally not bothered by the Fierz 
ambiguity as we do not introduce auxiliary fields. In general the sources in question 
do not need to be local. However, for an interpretation as a potential this is useful, 
but even with this restrictions the 2PI effective action has its own problems. Simple 
approximations are often unbounded from below raising serious questions about the 
interpretation of the 2PI effective action. 

Nevertheless, we do not want to leave this model without at least one not com- 
pletely trivial application. We will investigate an interaction resembling the three- 
flavor three-color instanton interaction of QCD. This six-fermion interaction provides 
for a mechanism to have a first order chiral phase transition. Furthermore, it was 
conjectured to have a color-symmetry breaking vacuum BBS}. 
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Outline 



As the effective action is the central object of our interest, we will briefly review its 
basic definitions and some simple ways to calculate it in Chap.EJ In Chap.|H]we will 
present some very simple calculations of the critical coupling of the model Eq. (jl.ljl . 
In particular, we will perform a MFT calculation and encounter the Fierz ambiguity. 
Chap. 0] introduces the concept of partial bosonization, and clarifies the origin of 
the Fierz ambiguity. A description with scale-dependent bosonic degrees of freedom 
allow us to cure the Fierz ambiguity for the RG calculation in Chap. El The following 
Chap.lBlturns to the concept of the 2PI effective action to circumvent the problems of 
partial bosonization. As an application of the 2PI effective action we investigate an 
instanton-like interaction. In particular, we focus on chiral symmetry breaking and 
a possible breaking of color symmetry by a non-vanishing color-octet condensate. 
Finally, in Chap. [7|we will return to the question of a possible composite Higgs. In 
particular, we are interested in the possibility of a non-perturbative renormalizable 
"Standard Model". We will only sketch some of the ideas and possible problems, 
accordingly this will be more like an outlook. Chap. |H] summarizes and concludes 
this work. 



Chapter 2 

1PI Effective action 



The effective action T [13 tilt Li2f is a very useful tool in quantum field theory 
(QFT). It allows us to calculate interesting quantities like vacuum expectation val- 
ues, propagators and correlation functions more or less by simply taking (functional) 
derivatives. Indeed, we can promote a classical equation to full quantum status by 
replacing the action by the effective action S —>■ T and the fields by their expectation 
values — > ((p). E.g. the equation of motion becomes 

5T 



Knowledge of the effective action is equivalent to knowledge of the full quantum 
theory. From this one can already deduce that calculating the effective action is a 
quite difficult task and can usually be done only approximately. 

Before going into more detail let us briefly review the definition of the 
(1PI) effective action 1 , and some of its basic properties. In the following we will 
write (f) for the fluctuating quantum field and <fi = ((f)). We suppress all indices. In- 
deed, <p might also contain fermionic degrees of freedom. A typical (f> might therefore 
look like (a, a*, V M , A^, . . . ,if)i,if)j . . .) with several bosonic and fermionic species. If 
we keep track of the order of fields and differential operators, no problems arise from 
this notation. 

Moreover, we work in Euclidean space. That is why we have a minus sign in the 
path integrals instead of an % in front of the action. The transition to Euclidean 
time is usually done via a Wick rotation. We do not want to go into detail, here, 
however, for fermions there are some slight difficulties because the action is no longer 
necessarily Hermitian 5(J 5l| . 



1 1PI abbreviates one particle irreducible (cf. Fig. I2.1|l . 
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(a) 



(b) 



Figure 2.1: An example of a diagram which is 1PI 2.1(a)| and one which is not 2.1(b)[ 
The latter one can be split into two by cutting the line between the two bubbles. 



The generating functional (or partition function if one prefers the statistical 
mechanics language) of a quantum field theory is defined by the following functional 
integral 

Z\j\ = I V4>exp(-S{<f>] (2.2) 



Here, S[<f>] is the classical action and j is an external source. We recall that in our 
matrix notation j'0 = j d d xj(x)(f)(x). 

Using the generating functional, expectation values of fields (and products of 
fields) can be calculated by taking derivatives with respect to j 



fV4>4>(g)eM-S[4>}+34>) = _^S_m = 8W\j] 

fV4>exp(-S[4>]+j$) " Z[j}Sj(q) 5j(q) 



with 



(2.3) 
(2.4) 



W\j] = \n(Z\j}). 
Physical values are obtained at vanishing external sources e.g. 
The (1PI) effective action is now the Legendre transform of W, 

and depends on the expectation value of the field. Combining (|2.2jl . ()2.4j) . ()2.5j) and 
shifting the integration variable to <p = <p — <P we obtain the following very useful 
formula 

r[0] = - In / P0exp (-S\(j> + 0] + ^jP-A , (2.6) 



since 



(2.7) 



We note that due to the shift of the integration variable ((f)) = 0. Finally, we 
would like to comment on the notion of one particle irreducibility (1PI). This is 
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m = s[</>] + 




Figure 2.2: Perturbative expansion of the effective action. To be explicit we choose a 
theory which has a quartic interaction. The propagators are propagators in a back- 
ground field <p. For the example of a theory with a interaction the propagator 
in the background field would be (p 2 + m? + 

explained most easily in terms of Feynman diagrams. A diagram is one particle 
irreducible if it is impossible to split it into two parts by cutting an internal line (cf. 
Fig |2.1|) . The effective action is now the generating functional of the 1PI diagrams. 
We stress this point because later in Chap. 16.11 we will encounter a 2PI effective 
action. For a proof of this statement and further details about the 1PI effective 
action we refer to textbooks as e.g. [EH IE2 Eil EH- 



2.1 Calculating the Effective Action 

In the following we will shortly present some of the standard methods to calculate 
the effective action. Before we start let us mention that in this section we will not 
explain mean field theory as it will later (Chap. H} be shown to be a one-loop 
approximation of a modified theory with auxiliary fields. Instead, we will explain 
by example a common method to obtain it in the next chapter and postpone a 
somewhat more thorough discussion until we have introduced partial bosonization 
in Chap. |U 



2.1.1 Loop Expansion 

The loop expansion is a perturbative technique to calculate the effective action. It 
can be shown that the effective action is nothing but the sum of all 1PI vacuum 
diagrams in a background field (f) as depicted in Fig. 12.21 An easy way to obtain this 
expansion is to make a saddle-point approximation of Eq. ()2.6|) . In lowest non-trivial 
order one obtains, 

r[0] = 5[0] + ^STrln(^ 2 )[0]) + ---, (2.8) 



where 



(2.9) 
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In the diagrammatic language (Fig. 12. 2j) the second term of Eq. (|2.8j) corresponds to 
the one-loop diagram, and we have omitted the higher loop diagrams in Eq. ()2.8j) . 
The supertrace (STr) comes around due to our notation where bosonic as well as 
fermionic degrees of freedom are contained in 0. However, its effect is very simple 
as it only provides a minus sign in the fermionic sector of the matrix. 

An advantage of perturbation theory is that it can usually be constructed to 
preserve symmetries order by order in the expansion. However, as we will see in the 
next chapter it has severe shortcomings if we want to go into the non-perturbative 
domain (hence the name) where the coupling is not small. Indeed, we will see that 
it is unable to describe the interesting phenomenon of SSB. 



2.1.2 Renormalization Group Equations 

Originally devised as a tool to hide the infinities of quantum field theory the renor- 
malization group has shown to give us a much deeper insight into the scale depen- 
dence of physics. Especially, the understanding of critical phenomena has profited 
immensely. In addition, in the last fifteen years or so the renormalization group 
has been established as a powerful tool for doing actual calculations in the non- 
perturbative regime. It is the latter aspect on which we want to focus. 

The basic idea of renormalization group equations is to introduce scale dependent 
effective couplings. Roughly speaking we look at the physical system of interest with 
a microscope with varying resolution. E.g. consider a lattice where a particle with 
spin sits at each lattice site. With a very high resolution we can see every particle 
and can measure each spin independently. However, shifting to a lower resolution 
things become a little bit blurry and we can only resolve regions which already 
contain several particles. What we then measure is something like an effective spin 
(more or less the sum of the individual spins) of several particles combined. Indeed, 
the first papers jH; 5(| 57; 58|| on renormalization group equations considered spin 



models like the Ising model. 

There is a wide variety of means for putting this intuitive picture into a mathe- 
matical form. A very convenient picture to do this is the path integral formulation. 
For a theory with an UV cutoff A one can write the partition function as 



n 

p<A 



J dcf>(p)eM-S[<f>)l (2.10) 



where the symbolic notation shall indicate that we integrate only over momentum 
modes with p < A. 

We can now implement the idea of decreasing the resolution of our microscope by 
integrating out modes in a small momentum shell [A' = A — AA, A]. This procedure 
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averages over the small length scales we do not want to see, 

Z= J] / #(p)exp(-5'[0]) (2.11) 

p<A' ^ 

with 

exp(-S'[0]) = J] fd4>(p)exp(-S[<l>]). (2.12) 

A'<p<A^ 

Where the so-called Wilson effective action 5" is now integrated only over a smaller 
range of momenta. Considering a A' which is infinitesimally close to A one can derive 
an evolution equation for 5" in the form of a differential equation. 

Following this general procedure (often followed by a re-scaling p — ► Ap to 
recover the initial momentum range) allows to derive a variety of equations for dif- 
ferent physical quantities like the Hamilton operator, correlation functions, coupling 
constants etc.. In addition, we do not need to restrict ourselves to the sharp momen- 
tum cutoff indicated above. Indeed, a smooth momentum cutoff is often much more 
convenient, as the sharp one has a tendency to introduce non- localities in position 
space. 

All equations derived in this way are exact, i.e. if they could be solved exactly 
they would give the same results for physical quantities. Hence, they are called 
exact renormalization group equations (ERGE). In principle they are all equivalent. 
Nevertheless, in practical computations where we have to use approximations they 
usually differ. For a review and a comparison s. |59| . 

Now, let us go directly to a formulation for the effective action or, more precisely, 
the effective average action as introduced in 0; |6l| EH; |63| 64; 65|. Let us begin 



by noting that for the effective action the running cutoff A' becomes an infrared 
cutoff. This is due to the fact that the effective action contains the modes which are 
already integrated out, thus those in the range [A', A]. A' becomes a lower bound 
for the modes included in T\r and therefore an infrared cutoff. We remark that T^i 
depends both on the IR cutoff A' as well as on the physical UV cutoff A. The first is 
quite clear since A' simply measures how far we have progressed in integrating out 
modes. The dependence on the UV cutoff is actually more physical. E.g. if we want 
to describe interacting particles on a lattice a continuum description might be quite 
good for length scales larger than the lattice spacing. However, roughly speaking 
this model does not contain any momentum modes higher than the inverse lattice 
spacing. Placing the same particles on lattices with different spacing will obviously 
produce different results. Thus, we have a dependence on the UV cutoff (for some 
more details s. also App. EJ). Indeed, only in very special theories is it possible to 
send the UV cutoff to infinity and still obtain finite results. These are, of course, the 
renormalizable theories. In this case we can have A = oo. 
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A big advantage of the formulation in terms of the effective action is that IV 
has a direct physical interpretation. Having included all quantum fluctuations above 
A' we can now view IV as the "microscopic action" on the scale A' where we have 
averaged over volumes of a size ~ A'(~ d ). Hence, the name effective average action or 
coarse-grained effective action. In an ideal description we would be able to observe 
the change of the relevant degrees of freedom from one scale to another, e.g. at very 
small scales we would have quark and gluon degrees of freedom while at a larger scale 
we observe mesons and nucleons, and at even larger scales we would have atoms or 
molecules, putting the basic idea of the renormalization group into full effect. 

Let us now get started and derive an explicit equation. To establish more clearly 
its function as an IR cutoff we write k instead of A' or use the convenient t = \n(k). 

In order to achieve a suppression of the low-momentum modes in the functional 
integral we add an effective momentum-dependent mass term 

A^#] = \Jf{-v)Rk(p)^P) (2-13) 

to the initial action. The idea is to add a high mass to the momentum modes p < k 
and a small or zero mass to those with p 3> k, therefore effectively removing (or at 
least suppressing) the low-momentum modes in the functional integral. To render 
this more precise we demand the following constraints for the function Rk{p), 

1. lim R k (p)>0, 2. lim R k (p) = 0, 3. lim R k (p) oo. (2.14) 

p 2 /k 2 ^0 p 2 /k 2 ^oo k^oo 

The first condition is the statement that we want to suppress the low momentum 
modes by an additional mass term. Having a mass term for the zero-momentum 
modes has another very nice effect as it removes all IR divergences produced by 
massless particles. The second condition ensures that the high-momentum modes 
(high compared to k) are not suppressed and that the cutoff is removed in the limit 
k —>■ 0. As we will later see when we have the explicit expression for the flow equation 
it is useful to choose a cutoff that vanishes sufficiently fast (e.g. exponentially) in 
the UV to avoid UV divergences. The third condition ensures that at k — > oo no 
modes are integrated out. However, a comment is in order as we want to send k to 
infinity. Originally, we wanted to suppress all modes with p < k. But, if we have a 
finite physical cutoff A there are no modes with p > A. Thus, one might want to 
rewrite condition 3 as 

3'. lim Rdp) — > oo, 

and indeed this is a quite common form given and used e.g. in . However, 

both conditions are equivalent as we can use every bijective function k'(k) which 
maps [0, A] into [0, oo] to re-scale k such that either 3 or 3' is fulfilled. Using a smooth 
cutoff the notion of modes being included or not included is anyway somewhat blurry. 
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With a little bit of discretion in the choice of k'(k) we do not distort the picture of 
integrating out modes down to the scale k very much for k < A, making the choice 
between 3 and 3' a matter of convenience 2 . In this work we will choose the former 
since it makes some of the analytical expressions easier. More important for the 
interpretation of k as the coarse graining with modes p < /c not yet being integrated 
out is that the shift between small values of Rk and large values of Rk occurs roughly 
at p ~ k. This is not a necessary condition to define the flow equation but useful for 
the interpretation. We leave it at the somewhat rough statement 

4. Rk{p) large for p < k, Rk{p) small for p > k. (2-15) 

Having talked so much about the cutoff let us finally give an explicit example of a 
very convenient one introduced in 69l |. 



R k ( p ) = Z k (l-p 2 )<d(l-?Pj. (2.16) 

Here we have included a factor of the wave function renormalization Z k . This choice 
allows us to write 

Z k p 2 + R k {p)=Z k P{p), (2.17) 

and guarantees that the reparametrization invariance of physical quantities <fi a<p 
is respected. 

As discussed in the last few paragraphs the term ()2.13|) acts more or less like 
an additional mass term. In consequence we might get concerned about violating 
symmetries like chiral or gauge symmetry. In principle there are two possible ways 
to tackle this problem. The first is to construct a cutoff function Rk{p) which fulfills 
Eq. ()2.14l) but does not violate the symmetry. This is relatively simple for chiral 
fermions |61: Tot 71; z2|. An example is 




R k {p) = Zki>\\j , --l\Q{l-^), (2.18) 

which is more or less ()2.16|) adapted to chiral fermions. The second strategy (usually 
used for gauge theories) is to accept the fact that during the flow the symmetry might 
be violated by the cutoff and is completely symmetric only at the endpoint when 
the cutoff vanishes. Nevertheless, the symmetry is onl y h idden during the flow and 
reveals itself in modified Ward-Takashi identities 0; H El E3 E3 • 

Taking these 



into account it is possible to construct an invariant flow. 

Adding the cutoff to the initial action we get a set of different actions 
parametrized by the scale k, 

Sk[<t>] = S[<j>] + AS k [<j>]. (2.19) 

2 Some additional details concerning different possible UV regularizations are provided in 
App.O 
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Correspondingly, we obtain 

W k [j] = HZ k [j}) = J V4>exp(-S k [4>} + J4>) (2.20) 

and _^ 

<P={4>) = ^W k \j}. (2.21) 

Now, we want to introduce the effective average action by a modified Legendre 
transform 

r fe [0] = -W k \j[<f>]] + j[<j>]<f> - &5k[<f>], (2-22) 

where we have substracted ASk[(f>] in order to remove the cutoff effects from IV 
This is particularly clear in the case k — » oo. Consider the formula analogous to Eq. 

T fc = -ln yp0exp^-5[0 + 0] + ^l0-A5 fc [0]^ . (2.23) 

Due to the third condition in (j2.14j) the additional term exp(— A [</>]) in the func- 
tional integral acts like a functional <5(0)-function for k — > oo and we obtain 

r oo [0] = S[0], (2.24) 

the microscopic action without the cutoff. 

Since the cutoff vanishes for k — > (second condition) we have in addition 

r o [0] = r[0]. (2.25) 

Thus, the effective average action interpolates between the classical or bare action 
and the full effective action (cf. Fig. 12 .4jl . 

Now, let us come to the final piece, the ERGE which governs the evolution from 
k = oo to k = 0. Taking a derivative with respect to k, 

— (r fc [0]) = -(d k W k )\j]-{d k j)j,W k V] + {d k j)<f>-d k AS k [<l>] (2.26) 
= (d k AS k [4>]} - d k AS k [<j>] 

= ^(($ T (d k RS)-(4> T )(d k R k )(4>) 



i(sTr[((00 T )-(0)(0 T ))9 fc i? fc ]) 



Expressing this with 
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gives 



Since 



we have 



with 




d k T k = l -$Ti{W { k 2) d k R k ). (2.28) 



£ 

Jk = (T k + AS k )— (2.29) 

00 



1 =l^ T )(^]=Wf , (rf + ^) (2-30) 



^ - (2 ' 31) 

Inserting this into Eq. ()2.28|) yields, finally, the flow equation 

d k T k = isTr {(rf + R^dkR,} . (2.32) 

Defining the operator 

d t =(d t Rh)-£ r (2-33) 
we can obtain an even more compact form 

d t T k = ~STr[d t ln(rf + R k )}. (2.34) 

This looks quite similar to Eq. ()2.8j) . and indeed if we neglect the change of T k on 
the right hand side, integrate and use ()2.24j) at k — oo, t = — oo, we recover ()2.8j) . 

In short, substituting rjjj^ for and writing it in a differential form turns 
the one-loop expression ()2.8|) into an exact equation. The one-loop form of Eqs. 
421321), (f2~34]l is depicted in Fig. ESI 

In addition to being exact, Eq. ()2.32|) has two more nice features. First, due to 
the presence of R k the expression is IR finite. Second, for R k decreasing sufficiently 
fast in the UV, d k R k provides an UV regularization. The flow equation is therefore 
completely finite. Of course, those divergences must still be included in some way. 
While the IR divergences might reappear in the integration of the flow equation for 
k — > 0, the UV divergences have been absorbed in the initial condition for some 
finite k. Specifying the initial conditions of the flow at some finite k gives a special 
regularization scheme called the ERGE scheme (cf. App. IC.3|) . 



Having a flow equation is only part of the game. Since it is a functional differential 
equation, it is in most cases impossible to solve it analytically, and we remember 
once again that this would amount to solving the quantum field theory in question. 
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Figure 2.3: Depiction of the flow equation ()2.32)1 . The line with the shaded circle is 

(2) 

the full field dependent and IR regularized propagator 1^. + Rk- The dot denotes 
the insertion of dkRk- Taking functional derivatives with respect to the field <fi adds 
external legs, i.e. we obtain flow equations for the propagators and vertices. An 
example for the flow of the propagator is shown on the right side. The shaded circle 
denotes the full fc-dependent vertex. 

So, it will be impossible or at least difficult to avoid using approximations in most 
of the physically relevant cases. 

A consistent and systematic approach is the use of truncations. In a truncation we 
restrict the space of all possible actions (T k ), spanned by all possible combinations of 
field operators compatible with the symmetries to a (very often finite dimensional) 
subspace given by a subset of operators. The approximate flow equation now is the 
projection of the flow onto this subspace (cf. Fig. 12. 4jl . From this we can calculate 
flow equations for the coefficients (generalized couplings) in front of the operators. 
We stress that the approximate flow is only driven by the operators in the subset. 
An easy, nevertheless usually quite tedious, way to improve the approximation and 
to check for errors is to enlarge the subspace. Doing this successively we may find 
a "convergence" of the results, and we may be tempted to interpret this as the 
approach to the right result. Still, we should be careful with this as it may well 
be, that we have indeed convergence, but convergence to the wrong result. This is 
usually the case when we have missed a relevant operator. As the number of all 
possible linearly independent operators is infinite we can always add operators but 
still miss the relevant one. 

From a physics point of view it is clear that a good approximation should include 
all relevant degrees of freedom, i.e. the corresponding operators. However, as we 
discussed in the introduction the relevant degrees of freedom can change with scale, 
making it necessary to adapt the description during the flow. For our case of interest 
a step in this direction was taken in |42j and we will discuss it at length in Chap. El 

Finally, let us come to the role the cutoff plays in approximations of the flow 
equation. By construction, an exact solution has a r = T independent of the cutoff. 
But, the trajectory is not independent of the cutoff. As depicted in Fig. l2.4l it may 
well be that a certain choice of IR cutoff may bring the real trajectory closer to the 
subspace defining our truncation, therefore usually improving the approximation. A 
systematic study to exploit this possibility has been put forward in 0; [Til 0; 
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Figure 2.4: RG flow in the space of all action functionals. The thick line is the exact 
flow in the (here 3- dimensional) space, the thin line is the approximate flow in the 
2-dimensional truncation. We point out that the approximate flow does, of course, 
not coincide with the projection of the exact flow on the subspace of the truncation 
(dotted line), as the latter is driven also by the operator in the third direction. The 
red set of lines, shows the same but for an "optimized" IR regulator. The exact flows 
coincide only at the start- and endpoint. The optimized flow is generally closer to 
the plane of the truncation, and the approximate flow is improved. 

A more traditional approach would be to use the IR cutoff dependence as a 
measure of uncertainty, s. e.g. jj^ . 



2.1.3 Schwinger-Dyson Equations 



Schwinger-Dyson equations (SDE) |83|; |84j (for a review and some applications s. 
(HI; l86t were one of the first really non-perturbative tools in quantum field 
theory. For a theory with polynomial interactions up to <fi m they provide an (infinite) 
hierarchy of equations which connect a 1PI Greens function of order n (nth derivative 
of the effective action) on the one hand with a set of 1PI Greens functions up to 
order (n + m) on the other hand. 

In principle they are a consequence of the fact, that the functional integral over 
a total derivative vanishes, if the functional vanishes at the boundary (like in normal 
calculus) , 

= J X?0Xexp(-S[0]+j0) (2.35) 
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Figure 2.5: SDE for the propagator (in 4 -theory). Full propagators and vertices are 
depicted with a shaded circle while bare quantities are represented by a dot. This 
equation is exact. However, it involves the full 4-point function which is given by 
another SDE. A simple approximation would be to neglect the last diagram on the 
right hand side giving us a closed equation. 




where + is for bosons and — is for fermions, respectively. Examining the last line 
in Eq. ()2.35|) it becomes clear that the SDE's are the Euler Lagrange equations of 
quantum field theory. 

Due to the appearance of S on the right hand side of Eq. (|2.H5jl we have always 
exactly one bare vertex in every expression contributing to the right hand side. A 
typical SDE therefore looks like in Fig. 12.51 

In its basic form the SDE allow us to calculate only the derivatives of the effective 
action. At first sight this may not seem like a major weakness, but in situations 
where we have multiple solutions for the SDE we have no way to compare them 
without knowledge of the value of the effective action. However, the case of multiple 
solutions is one of the most interesting ones, as it usually signals the possibility for 
spontaneous symmetry breaking. In Chap. El we will discuss the 2PI effective action 



I46t 14 7t I48| as a remedy to this problem. 



As in the case of the ERGE most of the time we are unable to solve the complete 
set of SDE. A popular approximation scheme is to neglect all 1PI Greens functions 
starting from a certain order. This gives a closed set of integral equations. More 
generally, similar to a truncation for an ERGE, one can restrict the space of all 
possible T and hence its derivatives to a subspace. 

On the exact level the RG and SD approaches are equivalent in the sense that the 
propagator and higher N-point functions calculated using the flow equation ()2.32|) 
are also solutions of the SDE [lH • Nevertheless, once approximations are used 
the results will, in general, differ. 



Chapter 3 



A Simple Example: The NJL 
Model 



In this chapter we want to get a grasp of the problems associated with the intro- 
duction of bosonic composite fields by studying the model, Eq. in some very 
simple approximations. In particular, we use this as an opportunity to introduce 
MFT. 



3.1 Critical Couplings from Mean Field Theory 

A mean-field calculation treats the fermionic fluctuations in a homogenous back- 
ground of fermion bilinears = (ip (^j-^j 4>), 4>* = — (i> \ ^2j Vn = (4 , 1^i J ) 

and = (■07 / /y 5 ^>). It seems straightforward to replace in the four-fermion inter- 
action in Eq. (jl.lj) one factor by the bosonic mean field, i.e. 

(^/>) 2 -(^7W -> 20^(1 +7 5 )^-20>(l- 7 5 )^, 

(^7m7 5 ^) 2 - 2A^Yl 5 ^- (3.1) 
The partition function becomes then a functional of 0, V^, A^, 

Z$,V,A)= /DWexp^Sft^^y.i]), (3.2) 
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where S is given by (jl.lj) . with the replacements (j3.1|) . Self consistency for the 
expectation values of the fermion bilinears requires 

= \m-^)= l -K l ^Z, (3.3) 

% = (^HA^^lnZ, 

and similar for the other bilinear A^. Chiral symmetry breaking by a nonzero 
requires that the "field equation" (j3.3|) has a nontrivial solution. We note that 
Z[(f),V,A] corresponds to a one- loop expression for the fermionic fluctuations in 
a bosonic background. With T\ = — In Z the field equation is equivalent to an 
extremum of 

T (F) = J d A x ^2X a (f)*(f) + ^\ V %V» + \\ A A^ + rS F) . (3.4) 

A discussion of spontaneous symmetry breaking in MFT amounts therefore to a 
calculation of the minima of r( p ). 

As we already noted in the introduction this calculation can be done equivalently 
in the Yukawa theory (jl.4j) . (|1.5j) . The mapping of the bosonic fields reads = 
(h a / 0, = (hy/Hy) Vft, A^ = (Ha/ Ha) A*- Keeping the bosonic fields fixed and 
performing the remaining Gaussian fermionic functional integral yields precisely 
Eq. (|3.4jl . Mean field theory therefore corresponds precisely to an evaluation of the 
effective action in the partially bosonized Yukawa model in a limit where the bosonic 
fluctuations are neglected. 

We want to compute here the critical couplings (more precisely, the critical line 
in the plane of couplings X a , Ay) for which a nonzero expectation value 0^0 
indicates the onset of spontaneous symmetry breaking. For this purpose we calculate 
the mass term ~ 0*0 in r^ F ^ and look when it turns negative. This defines the critical 
couplings. We assume here a situation where the expectation values of other bosonic 
fields like or A^ vanish in the relevant range of couplings. It is then sufficient to 
evaluate r( F ) for V^ = A^ = 0. 

In a diagrammatic language Gaussian integration over the fermionic variables 
corresponds to evaluating the diagram of Fig. 13.11 We define our model with a 
fixed ultraviolet momentum cutoff q 2 < A 2 , such that the MFT result becomes 
(v 4 = 1/(32tt 2 ), x = q 2 ): 

T ( p = -Av A dxxln(x + hl<j)*<f)). (3.5) 
Jo 

From this one finds the mean field effective action 

r( p ) = r( F) + rS F) (3.6) 

= (/i 2 - Av 4 h 2 a A 2 ) 0> + const + O 



3.1. Critical Couplings from Mean Field Theory 



23 




Figure 3.1: Bosonic mass correction due to fermion fluctuations. Fermionic lines are 
solid with an arrow, bosonic or "mean field lines" ~ {ip V') are dashed. For 

use in Chap. |S] we have indicated an external momentum p. A MFT calculation 
corresponds to an evaluation for p = 0. 

where we have expanded in powers of <fi. The mass term turns negative if 

< (3 ' T) 

As it should be this result only depends on the ratio h 2 a / ^ = 2X a . 

We now want to determine the critical line in the plane of invariant couplings 
A CT , Ay from the condition ()3.7j) . i.e. 

Using the relation (jl.3j) we infer a linear dependence of A^. nt on the arbitrary un- 
physical parameter 7 whenever Ay 7^ 

(For numerical values see Tabs. l3~T1 and l3.2)l . This dependence is a major shortcoming 
of MFT. We will refer to it as "Fierz ambiguity". The Fierz ambiguity does not 
only affect the critical couplings but also influences the values of masses, effective 
couplings etc.. 

The origin of the Fierz ambiguity can be traced back to the treatment of fluc- 
tuations. A FT of the type (jl.2j) changes the effective mean field. In a symbolic 
language a FT maps (ip a ip a )(ip b i(:b) — > (ip a ipb)(ipbi J a) where the brackets denote con- 
traction over spinor indices and matrices ~ 7^ or ~ 7 s are omitted. A mean field 
i/j a ip a , appears after the FT as ip a ipb- From the viewpoint of the fluctuations one 
integrates out different fluctuating fields before and after the FT. It is therefore no 
surprise that all results depend on 7. 
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Approximation 


fit 

Ghap. 


7 = 


0.25 


0.5 


0.75 


1 


MFT 


i3~n 


39.48 


38.48 


37.48 


36.48 


35.48 


Ferm. RG 




41.54 


41.54 


41.54 


41.54 


41.54 


Bos. RG 




36.83 


36.88 


36.95 


37.02 


37.12 


Adapted Bos. RG 


rm 


41.54 


41.54 


41.54 


41.54 


41.54 


SD 


m 


37.48 


37.48 


37.48 


37.48 


37.48 



Table 3.1: Critical values A CT for Ay = 2 and for various values of the unphysical 
parameter 7 (with A = 1). In anticipation of Chaps. El we give also results for 
the (adapted) bosonic RG. Progressing from MFT to the bosonic RG and adapted 
bosonic RG the dependence on 7 decreases as more and more diagrams are included. 
The Schwinger-Dyson result is independent of 7 but contains no vertex corrections 
in contrast to the RG calculations. 



Approximat ion 


Chap. 


7 = 


0.25 


0.5 


0.75 


1 


MFT 




39.48 


29.48 


19.48 


9.48 


-0.52 


Ferm. RG 


run 


14.62 


14.62 


14.62 


14.62 


14.62 


Bos. RG 


ipi 


15.44 


13.39 


13.45 


15.55 


19.46 


Adapted Bos. RG 





14.62 


14.62 


14.62 


14.62 


14.62 


SD 


Ed 


19.48 


19.48 


19.48 


19.48 


19.48 



Table 3.2: The same 1 as in Tab. 13.11 but with Ay = 20. 



3.2 Perturbation Theory 



In order to cure the unpleasant dependence of the MFT result on 7 we will include 
part of the bosonic fluctuations in Chaps. |U and El Some guidance for the level of 
approximations needed can be gained from perturbation theory in the fermionic 
language. Since the four-fermion vertex is uniquely characterized by X a and Ay 
the perturbative result must be independent of 7 at any given loop order. The 
lowest-order corrections to the four-fermion couplings are obtained by expanding 
the one- loop expression for the effective action Eq. ()2.8)1 2 



Ar (i-ioo P ) = l gTr 

1 The negative sign for the critical coupling at 7 = 1 in the MFT calculation means that the 
system is in the broken phase for any positive value of \& in this calculation. 

2 We remember that in the full S*- 2 -* matrix we have a term from the S 2 /SipSip derivative (S^) 

and a term from S 2 /5ip8ip, accounting for a factor of 2 in the language with the normal trace. 
Moreover, the trace includes momentum integration and summation over internal indices. 



In (S< 2 > 



-Tr 



In S 



i(2) 
FF 



(3.10) 
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Figure 3.2: Perturbative correction to the four-fermion interaction. Solid lines with 
an arrow denote fermionic lines. The letters in the diagrams are given to visualize 
the ways in which the fermionic operators are contracted, e.g. the first diagram in 

the second row results from a term [(V , a VO(V'cVO][(V , cVO(V , b'06)]- Evaluating the 

i i 

diagrams for fc-dependent "full" vertices and IR regularized propagators the above 
set of diagrams gives us the flow equation of Sect. 13.31 



up to order (V^) 2 - For this it is useful to decompose according to 

S (2) =V + jr 



(3-11) 



into a field-independent part V (inverse popagator) and a field-dependent part T. 
The RHS can then be expanded as follows, 



O 



-STr 



(3.12) 



V 



T 



+ 



This amounts to an expansion in powers of fields and we can compare the coefficients 
of the four-fermion terms with the couplings specified by Eq. (jl.ljl . Only the second 
term on the RHS contributes to order (ipif)) 2 . The corresponding graphs with two 
vertices are shown in Fig. 13.21 From 



Ar (l-loop) = ^ A 2 | [4A 2 _ 4A(j(Aa _ 2Xy)] |^)2 _ ^ 5 ^y 

+ [ 2X a Xy + 4(A A - X V )X V ] [(^ 7 ^) : 
+ [-Xl + 2X a XA + 3X 2 v -2X A X v -X 2 A ] ($n^y 



(3.13) 



we can read off the corrections AA CT , AAy and AAa to the coupling constants. In 
order to establish that our result is independent of 7 we use the freedom of FT to 
bring our results into a standard form, such that = j^. Inserting next the 
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invariant variables (jl.3j) leads us to: 

A\ = 4v 4 A 2 (\l + 4\X v + 3X 2 v ) } (3.14) 
AX V = 2v i A 2 (X a + \ v ) 2 . 

In contrast to MFT the result does not depend on 7. 

The perturbative result, Eq. ()3.14|) . always leads to finite corrections to the 
coupling constants. Remembering that in the fermionic language the onset of SSB is 
marked by a divergence of the coupling constants, it becomes clear that we will never 
get SSB in perturbation theory. No critical couplings can be calculated. This is a 
severe shortcoming of perturbation theory which cannot be overcome by calculating 
higher loop orders. Only an infinite number of loops can give SSB. In the next section 
we establish how a renormalization group treatment can overcome this difficulty 
without encountering the Fierz ambiguity of MFT. A calculation of the critical 
coupling becomes feasible. Nevertheless, even this RG treatment has its limitations 
once the couplings diverge. In particular, it does not allow us to penetrate the phase 
with SSB. In Chaps. H] and 03 this shortcoming will be cured by a RG treatment in the 
partially bosonized language. In particular, we will see in Chap. El which diagrams 
are needed in order to maintain the independence of results on 7 in analogy to 
perturbation theory. 

3.3 Renormalization Group for Fermionic Inter- 
actions 

Let us now turn to an RG equation. More explicitly the ERGE for the effective 
average action discussed in Chap. 121 Neglecting all change on the RHS leads to the 
perturbative result of the previous section. Consequently, in this approximation we 
cannot observe SSB. For a better approximation we restrict to the terms specified 
in Eq. (jl.lj) but take all couplings to be explicitly fc-dependent. In the action (jl.lj) 
we have only local interactions. Expressed in momentum space the four-fermion 
interactions have no momentum dependence. This is often referred to as the local 
potential approximation (LPA) jH lil EsT]. 

Performing the decomposition into a field-dependent part and a field-independent 
part as in the previous section, 



r fc + Rk — Vk + Tk 



(3.15) 
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we obtain an expansion of the flow equation 1)2.34)1 

This is in complete analogy to the previous section, only written in a differential form 
and with /c-dependent vertices. We obtain a set of ordinary differential equations for 
the couplings: 

d t \k = -toj{ F) '\s)k*(? atk + 4A CT ,^ + 3Ay ife ), 

d t X Vtk = -Av A l i f )A {s)k\\ a ,k + \v,k)\ (3.17) 

in agreement with where the same model has been studied. The threshold func- 
tions Z^' 4 (for our conventions cf. App. IC~2l or [67] ]) originate from the momentum 
space integration over the IR regulated propagators and replace the factor A- in Eq. 
()3.14p . For our actual calculation we use a linear cutoff 3 [78] and adapt the thresh- 
old functions to our setting with fixed momentum cutoff q 2 < A 2 in App. 1(121 The 
dependence on s = k 2 /A 2 becomes relevant only for k > A whereas for k < A one 
has constants Z^' 4 = 1/2. Although useful for the comparison of the RG with MFT 
or perturbation theory, the explicit /^-dependence of the threshold functions makes 
the fixed momentum cutoff somewhat cumbersome. An alternative is to use the so 
called ERGE scheme for the UV regularization. The basic idea (for details s. App. 
IG.3)1 is to use standard threshold functions without a UV cutoff in the momentum 
integral and implement the UV regularization by specifying the initial conditions for 
Tk at some finite k — A. This has the advantage that threshold functions without a 
finite UV cutoff are not explicitly fc-dependent, greatly simplifying numerical calcu- 
lations. The price to pay is that it is in general impossible to compare non-universal 
quantities like critical couplings for different choices of the IR cutoff function Rk{p)- 

The fermionic flow equations 4 ()3.17)) do not depend on 7. In a diagrammatic lan- 
guage we again have evaluated the diagrams of Fig. 15.21 but now with /c-dependent 
vertices. In the RG formulation we only go a tiny step Ak, and reinsert the re- 
sulting couplings (one-loop diagrams) before we do the next step. This leads to a 
resummation of loops. Since Eq. ()3.17)) is now nonlinear (quadratic terms on the 
RHS) the couplings can and do diverge for a finite k if the initial couplings are large 
enough. Therefore we observe the onset of SSB and find a critical coupling. Since 

3 The threshold functions depend on the precise choice of the cutoff. For the very simple trun- 
cation used in this section this dependence can actually be absorbed by a suitable rescaling of k, 

cf. App.ES 

4 As discussed above, the perturbative result Eq. (|3.14|l can be recovered from Eq. i|3.17[) if we 
neglect the fc-dependence of the couplings on the RHS and perform the ^-integration. 
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Eq. ()3.17|) is invariant this critical coupling does not depend on 7! Values for the 
critical coupling obtained by numerically solving Eq. (|3.17j) can be found in Tabs. 
IOandl3~2l 

The next step in improving this calculation in the fermionic language would be 
to take the momentum dependence of the couplings into account (e.g. jsjj) or to 
include higher orders of the fermionic fields into the truncation. This seems quite 
complicated and at first sight we have no physical guess of what is relevant. The 
renormalization group treatment of the bosonic formulation in Chap. 0] seems much 
more promising in this respect. 



3.4 Gap Equation 

Let us finally turn to the SDE as the last method discussed in Chap. 121 For the 
model Eq. (jl.lj) the SDE, approximated to lowest order, is depicted in Fig. 13.31 It is 
a closed equation since only the bare four-fermion vertex appears. (Only higher order 
terms involve the full four-fermion vertex.) We write the full fermionic propagator 
Gp as 

G F \p) = G F ^ P ) + J: F ( P ) (3.18) 

with the free propagator Gfo and self energy E F . Using this one obtains a gap 
equation for the self energy which can be solved self consistently. To simplify the 
discussion we make an ansatz for the self energy: 

£ F = M f7 5 , (3.19) 

where the effective fermion mass M F obeys the gap equation 

M F = 8t> 4 [Ao- + Ay] / dxx F (3.20) 
Jo x + Nl F 

The onset for nontrivial solutions determines the critical couplings: 

F* + vL lt = ^. (3.21) 

This result is shown in Tabs. 13. 1| 13.21 and does not depend on 7, as expected for a 
fermionic calculation. We observe that the MFT result for the A CT coincides with 
the SD approach for a particular choice 7 = 1/2. However, in general MFT is not 
equivalent to the lowest-order SDE. This can be seen by computing also the critical 
coupling for the onset of SSB in the vector channel. The MFT and SD results do not 
coincide for the choice 7 = 1/2. This becomes evident if we use a vectorlike ansatz 
instead of Eq. (|3.19|) for the self-energy: 



S F = f. 



(3.22) 
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Figure 3.3: Diagrammatic representation of the lowest order Schwinger-Dyson equa- 
tions for the fermionic model Eq. 1)1.1)1 . The shaded circles depict the full propagator. 

Using this ansatz for the self energy we find for the onset of non trivial solutions 

[A CT + 3V] c y rit =2^X*- ( 3 - 23 ) 
This is different from the MFT result 

obtained from Eq. (j3.3j) by setting = = 0. Again it would be possible to find a 
choice for 7 which makes both results equal. But, in general, this will not be 7 = 1/2. 
We note that in this channel the dependence of the MFT result is even worse than 
in the scalar channel. 

At last, let us note that MFT and the SDE gap equation share a nice feature, 
both allow us to proceed into the region of broken symmetry. The self consistency 
conditions (j3.3j) and (j3.20j) are valid for non-trivial solutions and hence finite values 
of the bosonic condensate, whereas the simple flow equation ()3.17j) breaks down (the 
couplings become infinite) when we approach the phase with broken symmetry. 



Chapter 4 

Partial Bosonization I: Basic Idea 



The MFT calculation introduces "mean fields" composed of fermion - antifermion 
(or fermion - fermion) bilinears. This is motivated by the fact that in many physical 
systems the fermions are not the only relevant degrees of freedom at low energies. 
Bosonic bound states become important and may condense. Examples are Cooper 
pairs in superconductivity or mesons in QCD. For a detailed description of the 
interplay between fermionic and composite bosonic fluctuations it seems appropriate 
to treat both on equal footing by introducing explicit fields for the relevant composite 
bosons. This will also shed more light on the status of MFT. 



4.1 Calculating the Bosonized Action 

A method for introducing the desired composite fields is partial bosonization j35|; |3_ 
37 ; Hi); 3^], sometimes also referred to as a Hubbard Stratonovich transformation. 
Regardless of the name, in principle it is nothing else but the insertion of a nicely 
written factor of unity in the functional integral for the partition function. 



4.1.1 A Toy Model Bosonization (d = 0) 

Before going into the details, let us demonstrate the idea on a 0-dimensional toy 
model. In dimensions the field variable <j)(p) is replaced by a simple c- or Grassmann 
number x depending on whether we deal with fermions or bosons 1 . To be explicit 
let us consider an "action" 

S(x) = ^x T Mx + ^(x T Sx) 2 (4.1) 

■""Of course, if x is a Grassmann number it is necessary to have several different components Xi 
in order to have non-trivial interactions. 
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where M and S are matrices simulating the mass and kinetic terms (M), and the 
spin structure of the interactions (S). Since we live in dimensions the momentum 
or position space integral is trivial, i.e. absent. This action models a massive field 
with an inverse propagator M and a quartic interaction with coupling strength A. In 
the following we want to introduce an auxiliary "composite" field y for the operator, 



h 



O s (x) = -x'Sx, (4.2) 

m z 

where we have introduced the for the moment arbitrary constants h, m for later 
convenience. To study a composite operator Os(x) it is useful to introduce a source 
term kOs{x) for this composite operator in addition to the ordinary source term jx. 
The functional integral over the field variables becomes an ordinary integral over x, 
and the partition function reads, 

Z{j,k) = J dxexp(-S(x) + jx + kO s (x) + % 2 ), (4.3) 

where we have used the freedom to add a field independent term quadratic in k to 
the action. 

Using the translation invariance of the integral we can obtain the following rather 
trivial identity, 



1 = N J dyex V (-—y 2 ) =N J dyex V (- — (y - O s (x) + d) 2 ), (4.4) 

where iV is nothing but a normalization constant N = (^J dy exp(— ^y 2 ) j and d 

is for the moment arbitrary, but will be determined later. Inserting this under the 
integral in ()4.3j) yields, 

Z(j,k) = N J dxdy exj>(—S(x,y) + jx + ky), (4.5) 

S(x, y) = Vmi + ^(x T Sx) 2 + k-^x T Sx + °-k 2 + ^—y 2 

1 h 2 m 2 d 2 
+hyx T Sx H -(x T Sx) 2 + m 2 dy + hdx T Sx H h ky. 

2 m 2 2 

The first line looks promising, as it is the partition function for a theory with two 
"fields" x, y and an action S(x,y). The second line is still a mess which over and 
above depends explicitly on the sources k and k. However, remembering that we have 
introduced several arbitrary parameters we can choose those to our convenience, 

1 k 

a = -, d= -, k = k, (4.6) 

m z m z 
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simplifying 



1 1 / h 2 \ m 2 

S(x, y) = -x T Mx + - A + — (x T Sx) 2 + —y 2 + hyx T Sx. (4.7) 

2 2 V m z I 2 



Finally, employing the choice 



h 2 

= -\ (4-8) 



m 2 



cancels all quartic interactions of x, leaving us with a mass term for the "composite" 
field y, a Yukawa-type interaction between y and x in addition to the propagator 
term x T Mx for the "elementary" field x, 

A 1 Tfl 2 

S(x, y) = -x T Mx + —y 2 + hyx T Sx. (4.9) 

This also explains why we have introduced the constants h and m in the normaliza- 
tion of O s [x], Eq. flOp . 

Having accomplished the "partial bosonization" of our O-dimensional model we 
would like to comment on some rather technical points: 

1. Physically it is clear that with respect to the symmetries y should have the 
same transformation properties as the composite operator Os{x). From a more 
technical point of view this is necessary as we would otherwise be unable to 
perform the shift in the integration variable in Eq. (j4.4j) . 

2. In the derivation given above we did not specify if x is bosonic or fermionic. 
We can use the same procedure to introduce composites made up of fermions 
or bosons. However, we should be careful. If the integral in ()4.3j) is fermionic, it 
is convergent for all possible choices of A, because of the rules of Grassmannian 
integration. The integral over the auxiliary field y, Eq. ()4.4|) . is only convergent 
for m 2 > 0. This gives us the condition that A < 0, in order to render every- 
thing finite. For bosons, however, a A < leads to a divergence in Eq. ()4.3j) . 
So, naively our bosonization procedure works only for fermions and a certain 
region of the coupling constant. Although it is possible to circumvent these 
naive arguments by an integration along the complex axis, any interpretation 
of y as a bound state is still awkward. Therefore, we will restrict ourselves to 
stable potentials of the composite field (s. below), i.e. the integration over the 
composite is convergent. 

3. It is not necessary that the integration over the auxiliary field is Gaussian as 
in Eq. ()4.4|) . Indeed, we can replace the term — ^y 2 by any potential — V(y) 
as long as, 

V(y) 

V(y)>c, lim — ^-^oo, (4.10) 
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i.e. V(y) is bounded from below and grows sufficiently fast for y — > oo. This 
allows us to absorb also higher order interaction as e.g. a term — k(x t Sx) 4 by 
a term Kj^y 4 in V(y). In general, 

-F(O s (x)) - -F{O s (x)) + F(O s (x) - y), (4.11) 

such that the purely fermionic terms are cancelled. Of course, this leaves us 
with non-linear Yukawa couplings, as e.g. ~ y 3 x 2 and other complicated inter- 
actions ~ y 2 x 4 or ~ yx e . 

4. As can be seen from (j4.11j) we can also treat terms linear in the composite op- 
erator Os(x), removing e.g. all parts ~ x T Sx from ^x T Mx. A typical example 
for this would be the translation of a fermionic mass term mipifj into a source 
term j(p for a boson corresponding to i/ji/j. 

5. Although it is the most common case, it is not necessary that the composite 
operators that we want to bosonize are made up of exactly two field operators. 
In principle, they can contain an arbitrary number of fields. The composites 
can even be fermionic. 



Using the translation invariance of the integral as in Eq. (|4.4|) is the simplest 
but not the only possible way to obtain an identity useful for the introduction 
of composite fields. In general, any identity 

1 = Nexp(F(O s (x))) [ dyexp(-V(x,y)), (4.12) 



can be used to cancel a part F(Os(x)) in the initial action. However, the direct 
interpretation y ~ Og(x) will usually be lost. Of course the V(x, y) in Eq. 
()4.12j) is far from unique. One possibility is always V(x, y) = F(Os(x) — y) as 
obtained in (j4.11j) . In practice it quite difficult to find a V(x, y) with a suitably 
simple form like V(x, y) ~ xy + V'(y). 

7. Sometimes, it might seem useful to add some form of interaction, e.g. ~ y 4 , 
between the composite fields to the bosonized action ()4.9|) . We can then use 
the argument of HI in a backward way to determine what kind of (higher order) 
interactions this would introduce into the initial unbosonized action. 



8. We can recover the initial action by performing the integration over the aux- 
iliary field y, 

exp(-S(x)) = J dyexp(-S(x,y)). (4.13) 
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4.1.2 The Fierz Ambiguity 

Already in this simple model we can get a grasp how the Fierz ambiguity arises. Let 
us take a look at the "four-fermion interaction" in Eq. 1)4.1)1 . 

~ {x Sx) —Si 1 i 2 Si 3 i 4 Xi 1 Xi 2 Xi 3 Xi 4 — (S* ® S~)i 1 i2i 3 i 4 Xi 1 Xi 2 Xi 3 Xi 4 (4.14) 

1 A 

n '1'2'3'4 '1 W >3 14) 

where the LHS defines A. We can now permute the Xi, e.g. let us exchange Xi 2 and 

Ajj 121314 ^iv^i% ^13 A-il 1213 lA^i\ ^ ii^iz"^i2 ^1112*314 "^iu^iz J (4-15) 

where 

KhisU = ± A iiW2, (4-16) 
and the sign is + for bosons and — for fermions, respectively. Now, let us assume 
that A' has a decomposition (this assumption is nearly always fulfilled), 

A' = Xtt'(T <g> T') 7^ X(S <S> S). (4.17) 

Accordingly, we would bosonize the RHS of Eq. ()4.15j) with fields corresponding to 
the operators Ot{x) and the coupling strengths Xtt 1 - Hence, we obtain a different 
set of composite fields and coupling strengths for the identical action. In general, 
we can perform an arbitrary permutation of the X{, and we can obtain not only two 
but several different bosonized actions. This is even worse for higher order (e.g. x 6 ) 
interactions. 

Comparing Eq. (j4.17j) with the Fierz identity Eq. ()B.5j) it becomes clear that an 
exchange of fields like in Eq. 1)4.15)1 is a Fierz transformation. Since the bosonized 
action may look quite different, it is no big surprise that simple approximations 
might yield different results. This is what we call "Fierz Ambiguity" . 

One might wonder about the fact that different A, A' describe the same (un- 
bosonized) action. For fermions this is quite easy to understand. Due to the Grass- 
mann identity only the completely antisymmetric parts of A give 

non-vanishing contributions to the action. Hence, all 

A = A + S (4.18) 

yield identical action, as long as S is a sum of terms which are symmetric in at least 
two indices. If we want, we can choose A to be completely antisymmetric. Any Fierz 
transformation described above can be obtained by adding a suitable chosen E. The 
problem is that a non-vanishing E usually does not give a vanishing contribution in 
the partially bosonized action. As there is great freedom in choosing E we can get a 
nearly arbitrary bosonized action. For bosons the story is essentially the same, only 
one has to replace symmetric by antisymmetric and vice versa. 
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4.1.3 The Case d > 0, MFT Revisited 



It is straightforward to generalize the procedure described in the previous section to 
the case of d > 0. Indeed, the change is more or less only a matter of semantics, as 
we replace functions by functionals and integration by functional integration, 



d^V. 



(4.19) 



To demonstrate this let us repeat the procedure for the action (jl.lj) 2 . Introducing 
bosonic fields for the composite operators corresponding to scalar, vector and axial 
vector bosons, 



2/4 
we obtain, 

Z =J V^ViPexp(-S 

with 



/4 



(4.20) 



V^V^V^VV^VA^M^MvHa^V (-SM) (4-21) 



My 
M A 



exp 
exp 
exp 



+ ^(l + 7 5 )^) (0-^(l-7 5 )^ 



2 V 



A„ 



(4.22) 



Collecting the terms in the exponentials and using Eq. (jl.5|) as the equivalent to Eq. 
(|4.8jl . the four-fermion interaction is cancelled. As expected, it is now replaced by 
mass terms for the bosons and Yukawa couplings between bosons and fermions as 
given by the expression 1)1.4)1 . We note, that the bosons do not yet have a non-trivial 
kinetic term and the propagator is simply -y. 

Having arrived at the partially bosonized action (jl.4)l for our model (jl.l)) we 
can use it to gain new insight into MFT. The action Eq. (jl.4j) is quadratic in the 
fermionic fields, hence the functional integral over the fermionic degrees of freedom 
is Gaussian and can be done in one step. As we have seen in the previous chapter 
this leads exactly to the MFT results. More precisely, we understand now that for 
different choices of 7 the MFT treatment leaves out different bosonic fluctuations. 

In this context we note that the condition (j4.1()j) restricts the possible couplings 
to X a , Xy, A a > 0. In the invariant variables this restriction translates to A CT , Ay > 
and for 7 it implies < 7 < 1. 



2 For simplicity we skip the introduction of bosonic sources. 
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Incidentally, we note that in the case of a four-fermion interaction we have an 
alternative to Eq. (j4.13j) to recover the initial action from the partially bosonized 
one. Instead of integrating over the bosonic auxiliary fields we can simply solve the 
classical field equations for the bosonic fields in terms of the fermionic fields and 
reinsert them in the partially bosonized action. Starting from Eq. (jl.4j) this returns 
us to Eq. (fl~T|) . 

The crucial advantage of the bosonic formulation is that it can easily be gener- 
alized. For example, the bosonic bound states become dynamical fields if we allow 
for appropriate kinetic terms in the truncation, i.e. 

Ar kin = / d'xiZ^dy+^V^V^ + ^A^ + p^id.V^ + i^id.A^ 2 ) 

J 4 4 lay J 

(4.23) 

with 

= dpVy - d u V^ A^ = d^A v - d u A fl . (4.24) 

Also spontaneous symmetry breaking can be explicitly studied if we replace fi^4>*<j) by 
an effective potential U(<j)*<f)) which may have a minimum for <fi ^ 0. This approach 
has been followed in previous studies [^3 : IS We remark that for those 
terms to be present in the effective action it is not necessary for them to be present 
in the initial (bosonized) action. They naturally receive non-vanishing corrections 
by loop diagrams. E.g. the kinetic terms ()4.23|) get a correction from the diagram 
depicted in Fig. l3.1l with non-zero external momentum. Nevertheless, it is instructive 
to investigate how such terms would look like in the unbosonized language. For the 
potential terms this has already been discussed in Sect. 14.1 .11 i.e. Eqs. (j4.11|) . (j4.12j) . 
for the kinetic (derivative) terms we will do this in the next section. 



4.1.4 Beyond Pointlike Interactions 

So far our bosonization procedure seems relatively simple, and it is. However, we 
should mention that above we have bosonized only the very special case of a pointlike, 
i.e. local four-fermion interaction 3 , 



\V(x)V(x)V(x)V(x) = I \^{pi)^(p2)^(p4)^ip3)S(pi+P2+P3+Pi)- (4.25) 

Pl,P2,P3 ,P4 



°This section discusses some details needed in Sect. 15.41 and can also be read then. 
In this section we suppress all internal indices. In particular indices distinguishing between 
ip and if). A is a matrix with four such indices. All problems connected with the internal indices 
are completely analogous to the previous sections. Therefore, we allow ourselves to be somewhat 
sloppy concerning the internal indices, simplifying the notation. 
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Figure 4.1: Typical diagram contributing to an effective four-fermion interaction 
like in Eq. ()4.25|) in QCD. The solid lines denote fermions with propagator Pp , 
the wiggled lines gluons with propagator Pg . The labelled arrows denote the mo- 
mentum flow. The diagram suggests that the contribution ~ f Pp 1 (q) Pp 1 (q + Pi — 

P3) Pb 1 (q + Pi) Pr 1 (Q ~~ V2) to the effective four-fermion interactions is not a constant 
but depends, at least on some combination of the external momenta pi, P2, P3, Pa- 

This is by no means the most general form of a four-fermion interaction. Giving up 
locality, A can become an arbitrary 4 function of the four momentum variables, 

A ->• \(pi,P2,Pa,Pi)- ( 4 -26) 

At first, giving up locality sounds like a big step not to be treated lightly. How- 
ever, we should remember, that we frequently use those four- and multi-fermion 
interactions not as a fundamental interaction but to model an effective interaction 
at some intermediate scale. An example are the four-fermion interactions used to 
model QCD at low energies. As an example, a diagram contributing to lowest order 
is depicted in Fig. 14.11 

So, what can we do about the bosonization of those awfully complicated inter- 
actions? In the previous section we have considered local operators of the form (C 
is a constant matrix in the space of internal indices), 

0[*](ar) = ^ T (x)C^(x) = I ^ T \y)C^{z)5{x - y)5{x - z) = f 0[*](p) exp(ipx), 

J x,y Jp 

(4.27) 

with 

0[*](p)= I V T (qi)Cy(q 2 )5( qi + q 2 -p). (4.28) 

J Ql,Q2 

Keeping in mind our physical picture of a bound state, we find that ()4.27|) is 
very restrictive. Indeed, for a physical bound state we would expect that the 
particles of which the bound state is composed are usually not at the same 

4 Of course, it must have the right transformation properties under Lorentz transformations. 
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place, but smeared out over a certain region of space. Therefore, we replace 
8(x — y)S(x — z) — > g(x — y, x — z) in Eq. ()4.27jl . i.e. the "elementary particles" need 
no longer be located exactly at x but they can be somewhat distributed around x. 
In a momentum space formulation we find, 

Om(p) = hG(p) [ ^ T (p 1 )C^(p 2 )g(p 1 ,p 2 )5(p 1 +p 2 -p). (4.29) 



establishing that the momentum p of the composite operator is the sum, p = pi +p 2 , 
of the momenta of the "elementary particles" , as it must be for a bound state. The 
function g{pi,p 2 ) is the so called (amputated) Bethe-Salpeter wave function 
The function hG(p) is a generalization of the factor -\ in Eq. (|4.2jl . and it serves 
the same purpose namely, to obtain a simple form with Yukawa coupling ~ h while 
keeping the direct relation (ft=0 between the bosonic fields (ft and the composite 
operator 0[\1/], as we will see below. 

Proceeding along the lines of the previous sections we insert a functional integral 

1 = M J VcfteM- J p \f(-v)G-\p 2 )(ft{p))- (4-30) 

Shifting the functional integral by the operator ()4.29|) we find that we can absorb a 
four-fermion interaction of the following form, 

KPuP2,P3,Pi) ~ h 2 g(p u p 2 )G(s)g(-p 3 , -p 4 )C T <g> C, s = (p x + p 2 f = (p 3 +p 4 ) 2 , 

(4.31) 

in a contribution 

S x [%<l>] = [ [hg{p 1 ,p 2 )^{p l )C<p T {-p 1 -p2)V(P2) + <ft T (pi)G-\pj) ( f ) (p 2 )5(p 1 +p 2 )] , 

Jpi,P2 

(4.32) 

to the partially bosonized action. It is clear that G(s) plays the role of the bosonic 
propagator. Hence, G(s) should have an appropriate pole structure in the complex 
s-plane, e.g. 

G(s) ~ -i— . (4.33) 
m 2 + s 

Let us summarize this in the following properties: 

1. Eq. (|4.31j) with the pole structure given by ()4.33j) is the most general momen- 
tum structure for the four-fermion interaction ~ A we can absorb in a single 
bosonic field and an action quadratic in the bosons. We can only bosonize 
four-fermion interactions factorizing into two pairs of momenta. Usually, con- 
tributions to A like the one depicted in Fig. 14.11 do not factorize completely, 
therefore bosonization is usually only an approximation. On the other hand 
factorization of the four-fermion interaction signals the onset of physical bound 
states and can be checked numerically [64| . 



4.1. Calculating the Bosonized Action 
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2. Permutation of the fields (corresponding to Fierz transformations) permute 
the momentum variables on the RHS of Eq. (J4.31)) . This allows us to 
absorb momentum structures with poles in the t- respectively w-channels 
{t = (pi - p 3 ) 2 = (P2 - Pi) 2 , u= (pi - Pa) 2 = {p 2 - p 3 ) 2 , s, t, u are the Man- 
delstam variables). 



3. Turning the argument of El around, we can determine the "correct" Fierz 
transformation by an examination of the momentum structure (poles!). 



Since it might help us to resolve the whole mess of the Fierz ambiguity, let 
us illustrate the third point by calculating an example. In addition this will also 
demonstrate how a momentum dependence of the four-fermion interaction and a 
wave function renormalization for the composite bosons are connected. 

Starting from the action (jl.lj) extended by the kinetic terms f!4.23|) let us calculate 
the corresponding purely fermionic action. For simplicity we take Zy = Za = hy = 
Ha = 0, i.e. we have no vector and axial vector bosons. In momentum space we then 
have, 



Si 



b / S(qi + q 2 ) 
91,92 



(/4 + Z a ql) ( jf{q 1 ) ( j>{q 2 
1 + 7 5 



(4.34) 



'pi,P2 



^(p2)(j>(p2)S(qi - Pi -P2) 



ha(p*(qi)ip(p4) ( n 7 ) ip(-Ps)S(q 2 -P3 -Pi) 



P3,P4 
2 1 ry 2\ 



S(qi + q2) (Vv + Zvq-L 



<t>*(qi) + 
x ( - 



h. 



1 + 7 5 

2 i>(-pi) ( — n — ) V'(P2)*(gi -pi -P2) 



/<3,P4 ^ct + ^<tQ2 

8(pi + P2 -P3-P4) 



1p(p 4 



2 

1 -7 5 

~~ 2 

hi 



^(-P3)%2 -P3 -P4) 



Pl,P2,P3,P4 



2(^2+^(^+^)2) 
X [^{-Pi)^{P2)^{Pa)^{-P 3 ) - ^(-Pl)l 5 ^(P2)^(p4h 5 ^(-P3)]- 



After the usual shift in the integration variable we can perform the Gaussian inte- 
gration over the bosonic fields, removing the first term on the RHS of Eq. 1)4.34)1 . 
From the second term we can read of the four-fermion interaction, 



hi 



hi 



hi 



2{n* + Z a ( Pl + p 2 f) 2(/4 + Z a (p 3 + p 4 )2) 2(/i2 + Z a s) • 

(4.35) 
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In particular, the four-fermion interaction A only depends on the Mandelstam 
variable s, while it is constant in t and u. 

An FT permutes — P2 and P3. After relabelling the integration indices we find, 

K(P1,P2,P3,P4) = 0, (4.36) 

h 2 h 2 
4(//2 + ^(pi - p 3 ) 2 ) 4(/i2 + Z a t) 

First of all these coupling do not factorize in functions depending only on pi,p2 and 
P3,P4 respectively. Secondly, they have a pole in t which cannot be directly absorbed 
by bosonization. Therefore, this is not the "right" FT. A similar calculation for the 
vector and axial vector bosons would have resulted in four-fermion interactions de- 
pending on s in the vector and axial vector channels, respectively. While after an FT 
we would have interactions in all channels, scalar, vector and axial vector, but again 
with the "wrong" dependence on t which cannot be absorbed into bosons. Roughly 
speaking we have the following recipe, if the four-fermion interaction depends on s, 
bosonize, else if it depends on t FT exactly once and then bosonize. 



4.2 Bosonic RG flow 

Having talked at length about how we can obtain the partially bosonized action 
Eq. (jl.4j) it is time to put it to some use. Let us start with an RG calculation in 
a very simple truncation. The flow equations in the bosonic language are obtained 
in complete analogy with the fermionic formulation. In this section we restrict the 
discussion to a "pointlike" truncation as given by Eq. (jl.4|) with /c-dependent cou- 
plings. We will see (Chap. EJ) that we reproduce the result of the last section in this 
approximation if we take care of the fact that new fermionic interactions are gener- 
ated by the flow and have to be absorbed by an appropriate fc-dependent redefinition 
of the bosonic fields. 

It is instructive to neglect in a first step all bosonic fluctuations by setting all 
bosonic entries in the propagator matrix V~ l equal to zero. This removes all di- 
agrams with internal bosonic lines. Among other things this neglects the vertex 
correction Fig. 14.21 and therefore the running of the Yukawa couplings. Indeed, Fig. 
13.11 is the only contributing diagram and we recover MFT. One obtains the flow 
equations 

dtAj. = Sh% k v,k 2 l[ F) '\s), d t » 2 v>k = 8h 2 v>k v 4 k 2 l[ F) ' A (s), (4.37) 

9t^A, k = %h 2 Atk v A k 2 lf )A (s), 

d t h a>k = 0, d t h V)k = 0, d t h A ,k = 0. 



4.2. Bosonic RG flow 
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Figure 4.2: Vertex correction diagram in the bosonized model. Solid lines are 
fermions, dashed lines are bosons. There exist several diagrams of this type since 
we have different species of bosons. The momentum configuration indicated by the 
arrows is such that it gives a contribution ~ <P{—pi)'ip{—P2)il , {p3)- The pointlike limit 
employed in this section corresponds to an evaluation for p, = 0. 



As long as we do not consider the wave function renormalization (J4.23)) for the 
bosons, the flow can completely be described in terms of the dimensionless combi- 
nations 

~ ^V,k ~ ^\,k , . QO s 

'V,fc ft a A,k K 



Due to the constant Yukawa couplings we can integrate Eq. (|4.37jl . We find 
critical couplings: 



j-L^ fJj^ fJ^ 

' crit = 4^4, ToT? I crit = 4t>4, TTT? I crit = 4^4- (4.39) 



h 2 A 2 ' l-t ^ ^ 

These are, of course, the results of MFT, Eq. (|3.7j) . We note that in Eq. (|4.37jl the 
equations for the different species of bosons are completely decoupled. The mass 
terms do not turn negative at the same scale for the different species. Indeed it 
is possible that the mass of one boson species turns negative while the others do 
not. Such a behavior is expected for the full theory, whereas for the fermionic RG of 
Sect. 13. 31 all couplings diverge simultaneously due to their mutual coupling. However, 
no real conclusion can be taken from Eq. ()4.39|) because of the strong dependence 
on the unphysical parameter 7. 



Now, let us also take into account the bosonic fluctuations. This includes the 
vertex correction Fig. 14.21 and the flow of the Yukawa couplings does not vanish 
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anymore 



dthl k 



-32vj[ F ^(s)k>hl k 
Av4 n \s)k 2 h 
Avd[ F)A {s)k 2 h 



h 2 
u V.k 



h 2 
u A.k 



^V,k ^A,k 



(4.40) 



2^2 
V,k 



+ 2 



h 2 
t^V.k 



+ 



h 2 

u A,k 



2 U 2 
A,k 



hi 



Am 



+ 2 



h 2 

^V,k 



h 2 

a A, k 

^\k 



Using the dimensionless e's we now find: 



9ticr,k 

dt£v,k 

dt€A,k 



-2e CT , fc + 8 
-2iv,k + 8 
-ZtA.k + 8 



1 + 4e CT , fc 

2 + e v ,k 

2 — €A,k 



' 1 


1 


\^V,k 




1 








2e,7,fc 


£A,k 




1 


2e CT ,fc 





^i F) ' 4 (,), 



(4.41) 



The onset of spontaneous symmetry breaking is indicated by a vanishing of e for at 
least one species of bosons. Large e means that the corresponding bosonic species 
becomes very massive and therefore effectively drops out of the flow. 

For initial couplings larger than the critical values (see Tabs. IH.ll and 1.12)1 both 
e^fc and iy,k reach zero for finite t. Due to the coupling between the different channels 
they reach zero at the same t. At this point eA,k reaches infinity and drops out of 
the flow. This is quite different from the flow without the bosonic fluctuations where 
the flow equations for the different species were decoupled. The breakdown of all 
equations at one point resembles 5 now the case of the fermionic model discussed 
in Sect. The 7-dependence of the critical couplings is reduced considerably, 
as compared to MFT. This shows that the inclusion of the bosonic fluctuations is 
crucial for any quantitatively reliable result. Nevertheless, the difference between 
the bosonic and the fermionic flow remains of the order of 10%. 



4.3 Gap equation in the Bosonized Language 

Next, we turn to the SDE for the bosonized model (jl.4|) . They are depicted in Fig. 
14.31 We will make here two further approximations by replacing in the last graph of 
Fig. 14.31 the full fermion-fermion-boson vertex by the classical Yukawa coupling and 
the full bosonic propagator by /j,^ 2 . We remain with two coupled equations. 

This is an artefact of the pointlike approximation. 



4.3. Gap equation in the Bosonized Language 
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Figure 4.3: Diagrammatic representation of the lowest-order SDE for the partially 
bosonized model (Eq. (|1.4j) ). The shaded circles depict the full propagator, the circle 
with the cross is the expectation value of the bosonic field and the empty circle is 
the full Yukawa vertex. 



In a first step we approximate this equations even further by neglecting the last 
diagram in Fig. 14.31 altogether. Then no fermionic propagator appears on the right 
hand side of the equation for the fermionic propagator which only receives a mass 
correction for (0) ^ 0. Without loss of generality we take real such that Mp = h a (p 
and 

G F \q) = -h + K<p 1 \ (4.42) 



Inserting this into the equation for the expectation value we find 



4w 4 



a 2 i 

Pa Jo 



h 2 A 



x + ti^cj) 2 

For the onset of nontrivial solutions we now find the critical value 



(4.43) 



1.2/4 



crit 



8w 4 A 2 



crit 



(4.44) 



which is the (ambiguous) result from MFT given in Eqs. ()3.7|) and ()3.9|) . This is 
not surprising since this exactly is MFT from the viewpoint of Schwinger-Dyson 
equations. Indeed, Eq. ()4.43|) is precisely the field equation which follows by differ- 
entiation of the MFT effective action (|6.19|) with respect to 0. 



r( p ) = /i^ 2 - 4v 4 [ dxx \n(x + h 2 a (j) 2 ) 
Jo 



(4.45) 



In a next step we improve our approximation and include the full set of diagrams 
shown in Fig. 14.31 Using the same ansatz as before, the self-energy Sp now has two 
contributions 

S F = M F7 5 = /i CT 7 5 + Am F7 5 . (4.46) 

The first one is the contribution from the expectation value of the bosonic field 
whereas Amp is the contribution from the last diagram in Fig. 14.31 given by an 
integral which depends on Mp. Both in the equation for (0) and in the equation for 
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the fermionic propagator only Mp appears on the RHS. Inserting (0) in the graph 
Fig- ED one finds a gap equation which determines Mp: 



Once more, this can be expressed in terms of the invariant couplings and again we 
arrive at Eq. (JHUDD - 

We point out that, in order to recover the result of the fermionic SDE we have 
started with MFT and added diagrams. Therefore, the fermionic SDE (or the bosonic 
SDE with the extra contribution from the fermionic mass shift diagram) sums over a 
larger class of diagrams which contains MFT as a subset. This is evident in the lan- 
guage of statistical physics: MFT is the Hartree approximation, while the fermionic 
SDE is Hartree-Fock. 

Looking more closely at the two contributions to Mp we find that alone neither 
the contribution ~ <fi (which amounts to MFT as we have discussed above) nor the 
"fermionic contribution" Amp are invariant under FT's. Only the combination Mp, 
which is the fermion mass and therefore a physical quantity, is invariant. Indeed, 
changing the FT amounts to a redefinition of the bosonic fields. This allows us to 
choose bosonic fields such that Am F = 0. Taking 7=1/2 gives us such a choice of 
the bosonic fields. This explains why MFT gives the the same result as the purely 
fermionic calculation in this special case. 

In the next chapter we want to adapt this idea of a redefinition of the bosonic 
fields to the RG calculation, i.e. we want to do it continuously during the flow. 




(4.47) 



Chapter 5 



Partial Bosonization II: Scale 
Dependent Degrees of Freedom 



In the last Sect. 14 .31 we found that in the SDE formulation we can complete MFT by 
adding the mass-shift diagram for the fermions (cf. Fig. 14.3)1 . The mass-shift diagram 
is a contribution to the purely fermionic part of the effective action. Therefore, it 
makes sense to look for purely fermionic contributions in the RG, too. 



5.1 New Four-Fermion Interactions 



In our truncation of Chap.EJthe bosonic propagators are approximated by constants 
/i^" 2 . The exchange of bosons therefore produces effective pointlike four-fermion in- 
teractions. One would therefore suspect that this approximation should contain the 
same information as the fermionic formulation with pointlike four-fermion interac- 
tions. An inspection of the results in Tabs. 13.11 13.21 shows, however, that this is 
not the case for the formulation in the RG context. In particular, in contrast to 
the fermionic language the results of the bosonic flow equations still depend on the 
unphysical parameter 7. 

In fact, even for small couplings A the bosonic flow equations of sect. 14.21 do 
not reproduce the perturbative result. The reason is that at the one-loop level new 
quartic fermion interactions are generated by the box diagrams shown in Fig. 15.11 
A straightforward inspection shows that they contribute to the same order A 2 as 
the diagrams in Figs. 13. ll and 14.21 Even if we start from vanishing quartic couplings 
after partial bosonization, such couplings are generated by the flow. The diagrams 
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Pi 




P3 














(a) 



Pi 



!>■< 



Pi 



Pi 



(b) 



Figure 5.1: Box diagrams for the bosonized model. Again, solid lines are fermions, 
dashed lines bosons and vertices are marked with a dot. The diagrams generate 
new four-fermion ~ , 4 , {—Pi)'4 , {P2)' l l , {Pi)'4 ) {~Pz) interactions even for the model (jl.4j) 
without direct four-fermion interactions. 



in Fig. 15.11 yield 



dt\A,k 



(3 Xa = -v 4 l[ F)A (s)k 2 



A*CT,& ^A,k 

2Av^\s)k^^ 

/V.fc ^A,k 



2k- 2 ^(k) 



h 4 h 4 h 4 



^V,k 



^A,k 



(5.1) 



+ 2k~^(k) 



Here, 7(fc) is an in principle arbitrary function of scale determining the choice of FT 
for the generated four-fermion interactions. In other words, ^{k) allows for the fact 
that we can choose a different Fierz representation at every scale. We will make a 
special choice of this function (similar to the one made in sects. IH.2I and lH.Hj) namely 
we require 

Im = ^1_ Vfc, (5.2) 
1-7 



with e given in Eq. (|4.H8j) . The resulting equation d t (iv,k/^A,k) = fixes j(k). An 
improved choice of j(k) can be obtained once the momentum dependence of vertices 
is considered more carefully (cf. |42| and Sects. l4~l~4l I5.4J1 . 



5.2. Adapted Flow Equation: Solving the Ambiguity 
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5.2 Adapted Flow Equation: Solving the Ambi- 
guity 

An inclusion of the couplings into the truncation of the effective average action 
does not seem very attractive. Despite the partial bosonization we would still have 
to deal with multi-fermion interactions and the bosonic formulation would be of 
even higher algebraic complexity than the fermionic formulation. A way out of this 
has been proposed in There, it has been shown that it is possible to reabsorb all 
four-fermion interactions generated during the flow by a redefinition of the bosonic 
fields. In the following brief description of this method we use a very symbolic 
notation 1 . In Sect. 15.41 we will add some details on the momentum dependence of 
field redefinitions. 

Introducing an explicit /c-dependence for the definition of the bosonic fields in 
terms of fermion bilinears, the flow equation Eq. (|2.H2jl is modified 2 : 

d t T k = dtTklfa + -—d t (j)h- (5.3) 

Here d t Y k \^ k = d t T k \ is the flow of the effective average action at fixed fields. Shifting 
by _ 

d t <pk = (ipip) d t uj k (5.4) 

we find 

d t /j, 2 = <9 t /i 2 |, d t h = d t h\ + [i 2 d t u k , d t X = d t X\ - hd t u k (5.5) 
and we can choose u k to establish: 

d t \ = 0. (5.6) 

Instead of including running four-fermion couplings explicitly we therefore have to 
use only adapted flow equations for the couplings contained in Eq. (jl.4j) . 

Let us now apply this method explicitly to our model. Shifting 

_/l_ 7 5\ _/l+ 7 5 \ 

d t <f> = -4> ( — - — J i'dtu^k, d t <p = i) ( — - — J 4>d t u a)k , (5.7) 



^^We refrain from explicitly returning to the 0-dimensional toy model of Sect. E!T~T1 as we hope 
it is clear from Sect. 14.1.31 that for pointlike interactions the case d > involves no additional 
difficulties. Nevertheless, let us note for completeness, that the replacements S — > d, 4>k —> Uk, 
tpt/j — > x T Sx would bring us back to the toy model. 

2 It has been pointed out by Jan Pawlowski that after the appropriate modification of the infrared 
cutoff for the scale-dependent fields 0] > the flow equation Eq. (|2.32(l does not give the exact flow 
for dtT k \- However, in the simple approximation of this section the bosonic fields do not yet have 
an infrared cutoff. Therefore, we can still use l|2.32|l . 
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we have 
Requiring d t X - 



= d t X a ,k\ — h a ^dt^a,k, 

= d t \v,k\ — 2hv,kdtU>v,k, 9 t \A,k = 9 t \A,k\ — 2hA, k d t UJA,k- 
for all A's we can determine the functions uj: 



d t UJ a . 



dtUA,k 



(5.8) 



(5.9) 



h a ,k %hv,k 2hA,k 

with the /^-functions given in Eq. (j5.1|) . This yields the adapted flow equations for 
the Yukawa couplings 

d t K >k = d t h^ k \ + nl k d t uj a:k , (5.10) 
dthy,k = d t h V)k \ + fiy k d t u Vt k, d t h A , k = d t tiA,k\ + A,h^t^A,k- 
Combining Eqs. l|OT) . (JSU, (521), Q, (j5ZEDD determines 7(fc) 



7(A;) = 2^f ) ' 4 ( S ) 



1 



=2 

fc V,fc 



+ 



(4e c 



4e 2 ~ 



(5.11) 



Having fixed the ratio between i V)k and €A,k we only need two equations to 
describe the flow. We will use the ones for e a;k and Iv )k = (1 — 7)eV,fc 



^e CT fc = -2e CT , fc + 4 



(1 + 7 ) - 4(-2 + 7 + 2 7 2 )^ + 4(3 - 7 7 + 4 7 3 )^ 



22 



d t e v .h 



2e y , fe +4 



2^ 



(27 - i; 



e V,fc 



-V,fcJ 



(5.12) 



These equations are completely equivalent to the fermionic flow Eq. (|3.17j) . In order 
to see this we recall that the simple truncation of the form ()1.4|) is at most quadratic 
in the bosonic fields. We can therefore easily solve the bosonic field equations as a 
functional of the fermion fields. Reinserting the solution into the effective average 
action we obtain the form (jl.lj) with the fc-dependent quartic couplings 
— 1 11 

X ^ = 777^ 2 7t^> X v,k = t^=—- (5-13) 



2k 2 e a>k k 2 e v ,k ' k 2 e V)k 

Inserting this into Eq. ()5.12|) we find Eq. (|3.17j) . establishing both the exact equiv- 
alence to the fermionic model and the 7-independence of physical quantities. On 
a numerical level, we can see the equivalence from the critical couplings listed in 
Tabs. EH EH 

On this level of truncation the equivalence between the fermionic and the adapted 
bosonic flow can also be seen on a diagrammatic level. As long as we do not have a 
kinetic term for the bosons the internal bosonic lines shrink to points. On the one- 
loop level we find an exact correspondence between the diagrams for the bosonized 
and the purely fermionic model summarized in Fig. 15.21 This demonstrates again 
that one-loop accuracy cannot be obtained without adaption of the flow. 



5.3. Trouble With the LPA, an Example 
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c 




c 






Figure 5.2: Summary of all diagrams encountered in the previous sections. There 
is a one to one correspondence between the diagrams of the bosonized model (first 
row) and the purely fermionic model (second row). Solid lines with an arrow denote 
fermionic lines. The letters in the diagrams are given for visualizing the ways in 
which the fermionic operators are contracted, e.g. the first diagram in the second row 

results from a term [(V'aVO(Vv0c)][(Vv0c)(V'&V'b)]- Shrinking bosonic lines (dashed) 

i i 

to points maps the diagrams in the first row to the second row. In the approximations 
of sect. 14.21 only the first or the first two diagrams are taken into account. 

5.3 Trouble With the LPA, an Example 

So far, everything seems quite satisfactory. In our simple truncation we have been 
able to solve the problem of the Fierz ambiguity for the partially bosonized language. 
Moreover, the adaption of the bosonic flow is quite intuitive as it implements the idea 
of scale dependent degrees of freedom. But, not everything is as rosy as it seems. 
Looking a little bit more closely we notice, that although the critical coupling is 
independent of 7, the values of the individual bosonic couplings e CT , ey, Za are not. 
Only two of them are fixed by the flow equation (j5.12|) . while the third one can 
be chosen freely. In particular, we cannot determine from this truncation which 
type of boson will condense. This is completely analogous to the purely fermionic 
description, and in view of the equivalence of both descriptions not too surprising. 
It seems evident that this is a shortcoming of our present truncation. 

Thinking about enlarging the truncation, two possibilities come to mind imme- 
diately, a more complicated bosonic potential and a non-zero kinetic terms for the 
bosons. In this section we will consider the first possibility. As it turns out this does 
not solve the problem, i.e. we still cannot decide which type of bosons will condense. 
Hence, we turn to the alternative of kinetic terms for the bosons in Sect. 15.41 
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5.3.1 The Gross-Neveu Model 



To get a first impression it is instructive to study an even simpler model than Eq. 
(jl.lj) . the so called N = I 3 Gross Neveu model !95j in three spacetime dimensions, 

Ss#, V>] = J d 3 x + ^(^) 2 | (5.14) 

Actually, it is indeed the same model, only the number of spacetime dimensions has 
been reduced by one. The simplification as compared to Eq. (jl.lj) lies in the fact 
that in three dimensions we can use spinors with only two components (s. App. |B|) . 

Despite its simplicity it is still interesting in its own right. In particular it has 
a parity-like symmetry ip(x) — > —ip(—x), ip(x) — > ip(—x) which is expected to be 
spontaneo usly broken by a non-vanishing vacuum expectation value (ipip) for large 
enough G \9(t\ . 

Following the lines of Chap. 0] it is straightforward to obtain the equivalent 
partially bosonized action, 

— f ( — - m 2 } h 2 

S B [ip, $ = / d 3 x I tyfal) + ihcj)^ + — 2 L G = —. (5.15) 

We now want to study this model in a truncation which includes an arbitrary local 
potential V(4>) but no kinetic term ~ 9^09^0 (in the following we will suppress the 
integration over the spacetime coordinates), 

r k = tfiw + ih k <fiM> + v k (4>). (5.16) 

The parity like symmetry translates into <f>(x) — > — 4>{— x) in the bosonic language, 
restricting any bosonic potential V(4>) to even powers of <fi. 

Let us, for the moment, assume that during the flow we generate fermionic 
~ (tf)ip) n , bosonic ~ (f) n and mixed ~ (jf-^ipip)" 1 interactions. We neglect all other 
contributions as they lie outside of our present truncation. Hence, we write for the 
flow at fixed fields, 

^tT k \=^ t U k { ( t )k ^) = Y,<( ( t>k)Ui' + ^^\ , (5.17) 
where we have expanded the RHS in powers of ipip about the point — V ^ , and 



5 N is the number of different fermion species. 



5.3. Trouble With the LPA, an Example 
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Vki&k) = Allowing for a scale dependence of the field <fi k we obtain, 

d t T k = d t T k \ + 5 -^d t( t> k (5.18) 
= d t U k ((f) kj tfjip) + (ih k iptlj + V k ((f) k ))d t (f) k 

= al(cj )k ) + J2^k)^ + ^^) +(ih k ^ + V k \cj )k ))d t cj )k . 

n>l ^ ' 

Now, let us make a non-linear field redefinition 4 , 

ih< k \ %h k J 

n>l 7 



Inserting this into Eq. ()5.18)1 yields a simplified flow which only affects the potential 
in <p k , 



«9 4 r fe = a£(0 fc ) = d t U k ^ fc ,M^J . (5.20) 

At this point one might wonder why h k does not receive any corrections as in ()5.10|) . 
The reason for this is that we have used the additional freedom to rescale (f) with a 
contribution ~ (f) in dt<f), transforming any change in h k into a change of V(<f)). 

So far this looks quite appealing, as we succeeded in including all the complicated 
interactions contained in dtUk^^ip) into a simple truncation which includes a 
potential depending only on <fi. However, let us now show that we can simplify 
our result even further by adding a suitably written to d t U k {(j) k) %ljip). Indeed, 
starting from the action ()5.15|) it is possible to reduce V k ((p k ) to a mass term ^y<t> k - 

If V k (4> k ) = ^Y<p k at a given point, we have V'(4> k ) = ml(j) k , and the flow of the 
potential reads, 

dtV k {(j) k ) = a°M = d t U k ($ k , i^fo) = ( 5 - 21 ) 

Since we do not consider gravity, field independent terms in the effective action are 
of no physical significance. Dropping those and remembering that V((p) can contain 
only even powers of <f), it is no restriction to write 

d t U k L,i^<t> k ) =5>r0 2 ". (5.22) 

^ 7 n>l 



4 One may wonder if such a field redefinition has the right symmetry properties. A thorough 
inspection tells us that the given field redefinition is of the type dt4>k = dx ^^J''"^ ; where X 
is a singlet under all symmetries (the dots denote any type of additional fields). A derivative of 
this type belongs to the conjugate representation of <f> k and thereby to the same as cf>k since (j> k is 
self-conjugate. This also provides us with a recipe how we can obtain field redefinitions respecting 
the symmetries for more complicated fields, e.g. vector fields. 
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Now, let us use that the spinors have only two components, i.e. at any given point 
x we have only four different Grassmann variables, ipi2( x )> ^1,2 (%)■ Consequently, 
we obtain 

(${x)i){x)) n = 0, Vn > 3, (5.23) 
by use of the anticommutation relations. In particular, we have, 

(■1 \2n 

^ (^) 2n = d t U k (<p k M)- (5-24) 
, _ m kJ 

Performing field redefinitions as above, but for d t U instead of d t U yields, 

dtVM = d t U U k , i^A = tifrl, (5.25) 

i.e. our potential remains a mass term for all k as we have claimed above. 

At first, this seems very strange. Yet, as we know from Chap. H|we can, at least 
formally, remove the bosonic fields by performing the appropriate functional integral. 
Since we have no kinetic term for the bosons, and V(<ft) is local, this results in a 
completely local fermionic interaction. However, Eq. ()5.23|) tells us that the highest 
order local and purely fermionic interaction is (ijjip) 2 . Therefore, any action of the 
form (j5.16J) is equivalent to (j5.14|) and in consequence also to Eq. (J5.15)) . as long as 
we choose G and ^2 correctly. In other words, if we do not consider a non- vanishing 
kinetic term for the bosons, an inclusion of a full bosonic potential does not give us 
any more physical information than a simple truncation to a mass term or a purely 
fermionic calculation with a local four-fermion interaction. In particular, we cannot 
proceed into the SSB regime. 



5.3.2 General Discussion 

It is straightforward to extend Eq. 1)5.21)1 to more general cases with several differ- 
ent bosonic composite operators and corresponding bosonic fields, e.g., induced by 
fermions with more components. In fact, using our symbolic notation of Chap. |U 
the generalization looks like Eq. (|5.21|) . The discussion leading to Eq. ()5.25)) can be 
generalized, too, but it is somewhat more complicated, as we have to use the inverse 
function V 1 ' 1 of V to write down dtU . Assuming the existence of V'~ l we find, 

dAi^O^) = dtUk&iiOip) - J2 b k [V'~\-ih k ^Oil))\ n , (5.26) 

\n\>N 

where, for simplicity, we have employed a symbolic notation with the components 
<pi corresponding to the operators Oj. The 6's are defined as in Eq. f)5.24|) . but 
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n = (jii, . . .) is now a multi index. N is the number of components of the spinor 
if) and determines the highest-order monomial in if) which does not vanish by the 
anticommutation relations. E.g. for a Dirac fermion in four dimensions (four spin 
indices) with three colors we have N = 4 x 3 = 12. Using d t U to define the field 
redefinitions we can force all contributions n , |n| > N to vanish as in Eq. (|5.25|) . 

Typically, a potential up to order 4 is sufficient to describe at least the basic 
features of SSB. As most models have A > 4 we might be tempted to conclude that, 
with the exception of some special cases, the LPA works. However, this is not the 
case. We have already seen that we can modify d t Vk(<f>) by adding a conveniently 
written zero to d t U ((f) , ifjip) . But, Eq. ()5.23|) and its generalizations to the case of 
more spinor components are not the only way we can write a zero. Fierz identities 
like Eq. (jl.2)l provide another one. As we can see from the example of Eq. (jl.2)) they 
allow us to find non-vanishing c n such that 

c n $Oip) n = 0, m < N. (5.27) 

\n\=m 

An addition of this to dtU((f>, if>Oif>) can be used to eliminate terms with n and 
\n\ < N in the potential. The c n in Eq. ()5.27|) are not all independent, but typically 
we can eliminate at least one species of bosons completely from the potential. In 
our model and truncation this freedom is reflected by the arbitrariness of j(k) in 
Eq. (|5.1j) and 7 in Eq. (j5.12j) . respectively, e.g. choosing 7 = in Eq. (j5.12j) yields 
fiA = and effectively removes the axial vector bosons. 

Physically, our findings in this section tell us that in the LPA without any kinetic 
terms for the bosons we simply cannot decide which type of boson will condense. 
To do this we need additional information. Therefore, we will investigate (simple) 
momentum dependent terms in the effective action in the next section. 



5.4 Going beyond the LPA 

In Sect. 14. 1.41 we have seen that a simple kinetic term in the bosonized action gives 
a momentum dependent four-fermion interaction in the purely fermionic language. 
Moreover, we could absorb this momentum dependent four-fermion interaction into 
a boson only if we chose a certain Fierz transformation. This provided us with the 
information to decide which FT is the "right" one. Therefore, let us add the kinetic 
terms specified in Eq. f!4.23j) to our pointlike action (jl.4j) . Finally, let us impose for 
simplicity one more restriction, ay = a a = 1 on our effective action. This simplifies 
the expressions for the vector-boson propagators, similar to Feynman gauge in gauge 
theories. 
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Let us first outline the following complicated calculation, using a symbolic no- 
tation. In a next step we then discuss some points hidden in the notation and some 
details of the employed approximations. Finally, we comment on a few properties 
and give numerical results, while the whole set of equations and a more explicit and 
step-by-step calculation is given in App. |D] 



5.4.1 Adapting the Flow 

We restrict our truncation to lowest non-trivial order in p 2 , and expand all couplings 
up to this order. We write, 

» 2 (p)=/J 2 +p 2 Z + 0(p A ), <V(p) = d tf i 2 + p 2 d t Z+ (5.28) 
h(p) = h+ p 2 h^ + 0(p 4 ), d t h(p) = 4 0) + p 2 k~ 2 p {2) + ..., 
A(p) = A<°> + p 2 \^ + CV), d t X(p) = k~ 2 (3f } + p 2 k~ 4 p {2) + .... 

At first sight, in the partially bosonized language it seems reasonable to take 
fi 2 (p) as in (|5.28|) but restrict h(p) = h and \(p) = 0. However, in Sect. 15.21 we have 
seen that at least the latter is not a good approximation because the term A^ -* is 
crucial for restoring the invariance under FT of the initial (pointlike) interaction. 
Furthermore, we found that the flows of fx 2 , h and X^ all contribute to the same 
order to the effective four-fermion interaction (after integrating out the bosons). It 
seems natural that this is also true for the terms of order p 2 : Z, h^ 2 \ A 1 - 2 **. Hence, 
we consider all these terms on equal footing. 

Having chosen our truncation, we can calculate the flow equations. Since we 
did not add higher powers of field operators, the remaining task is to evaluate 
the diagrams depicted in Figs. 13. 1| 14. 2\ 15.11 The only difference to our previous 
calculations is that we have to consider non-trivial external momenta. 

As in Sect. 15.21 we want to keep the desired simple form of the effective action, 
i.e. h{p) = h and X(p) = 0, by choosing appropriate field redefinitions and neglecting 
terms of order 0(p A ). To do so, we shift 

dt<fik(q) = (i>^)(q)d t uj k (q) +d t a k (q)(j)k(q), (5.29) 
and, as in Sect. 15 .2\ we employ 5 , 

d t T k = d t T k \ + / J£*d t <t> k (q). (5.30) 

5 As mentioned in Sect. 15.21 after the appropriate modification of the cutoff, the flow equation 
(|2.32(l does not give the exact flow d t Tk\- Since we now have a cutoff for the bosonic fields, the use 
of Eq. H2.32J) is really an approximation. However, numerical tests performed in the appendix of 
|42| for a very similar case suggest that it is a very good approximation. So we will use it without 
further comment. 
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This results in the following changes of the flow equations, 

d t ^ 2 (q) = d t ^(q)\+2fi 2 (q)d t a k (q), (5.31) 
dth(q) = d t h(q)\ + h(q)d t a k (q) + [i 2 (q)d t uj k (q), 
d t X(q) = d t X(q)\ - h(q)d t uj k (q). 

We can now use the freedom in choosing the field redefinitions to enforce, 

d t X(q) = 0, d t h{q) = dth®, d tf i 2 (q) = d t fi 2 + 0(p*). (5.32) 

Roughly speaking we absorb the four-fermion interactions in the masses and Yukawa- 
couplings, and the momentum dependence of the latter ones in the wave function 
renormalizations for the bosonic fields (dtZ = 0). In particular, we keep the simple 
form with a momentum independent Yukawa coupling and no four-fermion interac- 
tion. 

Using the fact that we start with a constant Yukawa coupling and A = we can 
solve the equations ()5.32|) and find, 

V = 2d t a(0) = - dj ^ + 2 ~Y [dth {2) \ + V 2 dA {2) \ + Zd t \ {0) \] • (5.33) 

We use this equation to define the anomalous dimension. From Eq. (|5.29J) it is 
clear that d t a(0) modifies the overall normalization of 0, hence the wave function 
renormalization. Moreover, setting &th®\ = 0, dtX(q)\ = 0, it coincides with the 
original definition r\ = — ^jjr. Using this, we get pretty much the standard equations 
for the flow of the mass and Yukawa coupling, 

<V = W 2 + d t » 2 \, (5.34) 
dth = ^ v h + d t h\ + ^-d t \W\. 

We remark that together with the initial condition Z = 1 the conditions ()5.32j) 
automatically guarantee that the couplings are renormalized. 

5.4.2 Choosing the Momentum Configurations 

So far everything seemed relatively straightforward. However, looking more closely, 
we soon find that d t h(p) can actually depend on two and d t \(p) even on three 
momentum variables. This is in contrast to fi 2 (p) which depends only on p 2 . As we 
ultimately want to absorb those momentum dependencies in dtfi 2 (p) it is clear that 
we have to make an approximation such that d t h{p) and d t X(p) depend only on one 
momentum squared. To decide which of the possible momenta to choose, we look 
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at our example of Sect. I4.1.4| in particular at Eq. 1)4.34)1 . From this we can read of, 
that the bosons have the form, 

<!>{q) = f{q) ^{Pi)^{-P2)S{q-pi-p 2 )- (5.35) 

Jp\,P2 

Since we want to keep the simple form of the effective action, it is clear that the 
bosons have to keep this form, too. 

d t (f>k(q) = d t a{q)(j) k (q) + d t u k (q) J ip(p 1 )'4)(-p 2 )S(q - Pi - P2), (5.36) 

is the most general form that accomplishes this. With Eq. ()5.36)) we can give meaning 
to our notation, 

W0(?) = / 4>(pi)il>(-P2)S{q-pi-P2)- (5.37) 
Inserting this into Eq. (J5.30)) we find that the most general structures we can absorb 



arc, 



F{p x ) (f)(-pi)ip(p2)ip(-P3)S(pi +p 2 + Ps), (5.3* 

P1,P2,P3 



/ G((pi + p 2 ) 2 ) i>{-P\)il>{p2)i>{Pz)il>{-PA)KPi +P2-P3- Pi), (5.39) 

J Pl,P2,P3,Pi 

where F and G are arbitrary functions which can then be expressed in terms of dta k 
and d t uj k . 

Comparing the vertex correction ~ <P{~Pi)' l l , {~P2)'4 , {Pz) depicted in Fig. l4.2l with 
(15.38)1 it is clear that we have to restrict the momentum dependence to p\. A suitable 
configuration for the evaluation then is, (pi,P2,£>3) = (p, \p, \p)- 

Recalling that a Fierz transformation for the four-fermion interaction ex- 
changes P2 and — £>3 we can absorb either a dependence on s = (pi + P2) 2 
or one on t = (pi — p^) 2 (corresponding momentum configurations would be e.g. 

(Pi,P2,P3,Pa) = |(P,P,P,P) and (pi,P2,P3,pd = |(P,-P, ~P,P), respectively, 
cf. Fig. 15.1)) . Therefore, we have to ask which gives the better approximation. In 
principle, we would have to calculate d t X(s,t) (or even better d t X(pi,p2,P3,P4)) and 
test at every scale whether d t X(s, t) w d t \(t) or d t \(s, t) w d t X(s) is a better approx- 
imation. Analytically as well as numerically this is rather complicated. Therefore, we 
have adapted a much simpler scheme: we have always absorbed the dependence on t, 
i.e. we have evaluated the diagrams in Fig. 15. II and Fierz transformed the resulting 
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interaction once. There are two reasons why we believe that this is a reasonable 
approximation. First of all, in the pointlike limit (at the beginning of the flow), i.e. 

h 2 

/j 2 — > oo, h 2 — > oo, — = const, (5.40) 

one finds d t X(s,t) = d t X(t) exactly. Secondly, we have checked for various combina- 
tions of h and \x that in the vicinity of (s, t) = (0, 0) the dependence on t is usually 
(but not always) stronger than the dependence on s. 



5.4.3 Initial Flow and Numerical Results 

To get an impression of what we have achieved by all this let us take a look at the 
effective four-fermion interaction at the beginning of the flow 6 . 

As we want to start with action (jl.ljl where the four-fermion interaction is point- 
like, the kinetic terms in the partially bosonized action vanish and the renormalized 
couplings obey Eq. ()5.4()j) where the constants are given by Eq. fll.5jl . 

In Sect. 14.1.41 we have already calculated, 

h 2 h 2 h 2 

Hp) = ^— 2 = - 2 - V + -. ( 5 - 41 ) 

jj z + p z /i z /I 4 

where fi 2 and h 2 are constants in momentum space. Using our flow equations given 
in App. [0] an d using the properties of the threshold functions given in IC~2~3l we 
find (after some algebra) Eq. (j5.12|) for the flow of the p°-terms. For the flow of the 
p 2 -term in A CT we find, 

d t (ji) = lfirf ),4 (0)(A, + V) 2 , (5.42) 

and similar expressions for Xy and A^ which are invariant under Fierz transforma- 
tions, too. This shows that at least at the beginning we have no Fierz ambiguity up 
to order p 2 . 

Now, let us come to the numerical results. Numerically it is impossible to employ 
the pointlike limit exactly, therefore we have started with large values of /i 2 ~ 10 5 . 
In addition to the results of the full set of flow equations we have given some results 
for more simple approximations in Tabs. l5~T| 15.21 The first approximation, (1) + (3), 

6 Actually, in this section we make one more approximation: we have ignored the anomalous 
dimensions in the arguments of the threshold functions. In the pointlike limit this neglects terms 
of order A 3 . The flow equations in App . iDl include those terms, but they make the evaluation much 
more difficult. 
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Approximat ion 


Chap. 


7 = 0.1 


0.25 


0.5 


0.75 


0.9 


MFT 


im 


78.56 


77.96 


76.96 


75.96 


75.36 


SD 


urn 


76.96 


76.96 


76.96 


76.96 


76.96 


(1) 




76.32 


76.36 


76.42 


76.49 


75.54 


Ferm. RG = (1) + (2) 


I3~3l 


86.15 


86.15 


86.15 


86.15 


86.15 


(1) + (3) 


El 


76.43 


76.45 


76.49 


76.55 


76.58 


(1) - (3) 


EM 


83.97 


83.95 


83.92 


83.89 


83.87 


(1) - (4) 


EM 


86.17 


86.18 


86.20 


86.21 


86.22 



Table 5.1: Critical values X a for \y = 2 and for various values of the unphysical 
parameter 7 (with A = 1). To keep the table of manageable size we have abbreviated: 

(1) the pointlike contibutions to the mass and the Yukawa coupling (Figs. I3~T] 14.21) . 

(2) the pointlike contributions from the box diagrams (Fig. l5.ljl . (3) the contribution 
to the WFR from the purely bosonic diagram (Fig. 13. lj) and (4) the contribution to 
the WFR from the momentum dependence of the diagrams 14.21 and 15.11 We point 
out that differing from Tabs. 13. lj 13.21 we have employed a UV regularization by 
the ERGE scheme (cf. App. IC.3J) with the linear cutoff Eq. ()C.3j) which is better 
suited for numerical computations. This is why the values for the critical coupling 
are roughly twice of those given in Tabs. 13.11 13.21 since the critical coupling is not a 
universal quantity, and therefore scheme dependent. 

corresponds to the naive approach to the bosonized model, where all contributions 
from the four-fermion interactions and the momentum dependence of the Yukawa 
coupling are neglected. In the next step, (1) — (3), we have included the box diagrams, 
but only its constant parts, not the terms of order p 2 . However, indirectly we have 
included some knowledge of the momentum dependence as we have chosen the same 
FT as for the truncation (1) — (4) where we have included all terms of order p 2 . 

Moreover, we notice that the values for the critical coupling in the pointlike 
approximations are roughly twice of those given in Tabs. 13. 1\ 13.21 This is due to 
a difference in the UV regularization. In this section we have employed the ERGE 
scheme described in App. IC.3I Non-universal quantities can and do depend on the 
choice of UV regularization. In the pointlike approximation this yields exactly a 
factor of two in our case (for the pair of couplings (A CT , Ay)). In the more involved 
approximations this is not necessarily so, but the factor will still be somewhere 
around two. 

Aside from this, there is nothing new in the first four lines of Tabs. 15.11 15.21 
Comparing the pointlike truncations for the RG with the improved approximations 
of this section we find that the effect is of the order of 10%. Moreover, comparing the 
different non-pointlike approximations we find that the differences between them are 
of the order of 5% — 10%, too. While most of the Fierz ambiguity (more important 
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Approximation 


Chap. 


7 = 0.1 


0.25 


0.5 


0.75 


0.9 


MFT 


im 


74.96 


68.96 


58.96 


48.96 


42.96 


SD 


13.41 


58.96 


58.96 


58.96 


58.96 


58.96 


(1) 




53.16 


52.93 


53.32 


54.64 


55.88 


Ferm. RG = (1) + (2) 




58.83 


58.83 


58.83 


58.83 


58.83 


(1) + (3) 


E31 


53.89 


53.66 


54.00 


55.23 


56.37 


(1) - (3) 


EM 


58.14 


58.04 


57.88 


57.73 


57.64 


(1) - (4) 


EM 


61.60 


61.69 


61.82 


61.91 


61.94 



Table 5.2: The same as in Tab. 15. II but with Ay = 20. 

for large values of Ay, Tab. 15. 2 j ) is eliminated by including the pointlike contributions 
of the boxes (in the "right" FT), (2), for the absolute values of the critical coupling, 
the momentum dependence of the Yukawa coupling and the box diagrams is not 
negligible. 

On the more qualitative side we have checked for various values that for all 
truncations which include kinetic terms (last three lines in the tables) and values 
of X a slightly larger than the critical A^ 1 only the renormalized scalar boson mass 
turns negative, while the renormalized vector and axial vector boson masses remain 
positive. This allows the conclusion that the scalar boson will condense first, and we 
have a phase where only chiral symmetry is broken (at least in our approximation). 



Chapter 6 

Bosonic Effective Action (2PI) 



In the last two chapters we have mainly worked on improving the RG description 
in the partially bosonized language. It turned out, that it is necessary to include a 
wave function renormalization (WFR) for the bosons. Without a WFR we cannot 
determine the type of the bosonic condensate (e.g. if it is a vector or scalar conden- 
sate). Yet, consistent inclusion of a WFR leads to high algebraic complexity. This 
might be appropriate for a quantitative description. However, if we want to get a 
first, more qualitative, overview of a physical system this seems to be a little bit 
excessive. For this purposes the SDE or MFT approaches seem much more suitable. 

Both methods allow for a computation of the order parameter in systems which 
exhibit spontaneous symmetry breaking (SSB). However, while the SDE approach 
leads directly to the gap equation the MFT approach provides naturally a free- 
energy functional for the bosonic composite degrees of freedom introduced by partial 
bosonization via a Hubbard Stratonovich transformation (s. Chap. EJ. The field 
equation for this functional corresponds to the gap equation. Knowledge of the free- 
energy functional becomes necessary if the gap (or field) equation allows for solutions 
with different order parameters and the free energy for the different solutions has to 
be compared. The reconstruction of the free-energy functional from the gap equation 
is not trivial and the method used in j9?J for the case of color superconductivity may 
not always work. 

From this it seems that MFT is superior to SDE. Unfortunately, as we have seen 
in Chaps. El IH it has a severe disadvantage: partial bosonization is not unique and 
the results of the MFT calculation depend strongly on the choice of the mean field. 
Moreover, MFT only includes a subset of the SDE diagrams. 

Hence, we want to find a functional which has the SDE as its equation of motion, 
and which can be interpreted as a free energy. Such a functional is given by the 2PI 
effective action 0; E3 E3 
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In general, the 2PI effective action is a functional of fields and propagators 
r( 2P7 )[0, G\. However, for a purely fermionic system, all the information is already 
contained in r( 2P/ )[0, G\. r( 2P/ )[0,G] depends only on the bosonic variable G, and 
therefore we will call it Bosonic Effective Action (BEA) [q8( | . 

In this chapter, we want to calculate a simple approximation of the BEA for a 
general local multi-fermion interaction. Already for a four-fermion interaction the 
lowest non-trivial contribution to the BEA is of two-loop order. For a general n- 
fermion interaction we have an —loop structure. However, we will show that this 
can be reduced to a one-loop expression at the solution of the SDE, allowing for a 
comparison to MFT. 

As an application we want to study an interaction resembling the six-fermion 
inte raction generated by instantons in the case of three flavors and three colors 
93; Eq3 ESI Eol Eo3- In QCD this interaction i s of special interest as it is £/(!)- 



anomalous and solves the famous [/(l)-problem |104| . In the simpler case of two 
flavors instantons mediate a four-fer mion interactio n which has been investigated in 
works on chiral symmetry breaking |l05l Il06l ; Il07l ] and color superconductivity e.g. 

BE}. 

The effective interaction generated by the instantons does not only lead to inter- 
actions between color singlet effective quark-antiquark degrees of freedom (— > chiral 
symmetry breaking) but also between octets leading to the possibility of octet con- 



densation and spontaneous color symmetry breaking |2j| UB, |49j. In the following 
we will consider both possibilities. 



6.1 Bosonic Effective Action (2PI) 

To simplify the presentation we summarize all indices of the fermionic field in ip a . 
The index alpha contains all internal indices (spin, color, flavor etc.) as well as 
position or momentum. Furthermore it also differentiates between ip and ip. 

The partition function reads 

z [vJ] = J T>^exp(r] a ijj a + ^jap^a^p - S int [ip)) (6.1) 

where we treat all quadratic terms as a bosonic source term. 
We specify the interaction as 
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The usual generating functional of 1PI Greens functions in presence of the bosonic 
sources j is defined by a Legendre transform with respect to the fermionic source 
term rj: 

r F [,J>,j] = -W[ri,j]+V«il>« (6-3) 

where 

dW 

W = lnZ[ V ,j], ^ a =(^ a ) = — (6.4) 
Tp can also be obtained by the following functional integral 1 : 



rp#,j] = _m y Dip exp(ri a ip a - Sj[tjj + tp]), (6.5) 
1 > „ 

S 3m = —^japtpalpp + Si nt [lp). 

This form is especially useful to derive the SDE. Taking a derivative with respect 
to ip we find 

= -Jf3a 2 Va 2 (6.6) 

+ £ -^ff^ Q2 { (if ^ ■ ■ ■ (rP)-l ian + z a3 ... an + om j 

where 

F n = (n-2)(n-4)---2 (6.7) 

and Z summarizes all terms containing third and higher derivatives of V. These are 
terms which have at least two vertices. Taking another derivative with respect to ip a 
and evaluating at ip — we find the SDE: 

fPu = -u + zZ^ff^^Pr.L ■ ■ ■ FPr«La n + m 

In this chapter we are only interested in the lowest order. Therefore, from now on, 
we neglect Z, i.e. terms with at least two vertices. 

The "Bosonic Effective Action" (BEA) |98| , is defined by another Legendre trans- 
form with respect to j: 

T B [G] = -W[0,j]+jG, (6.9) 

G a p = 77— = (ry)^, 7777— =Ja/3- (6.10) 
OJ a/3 OLr a p 



1 Note that in this formula tp is shifted such that (tp) = 0. 
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Since IV is an even functional of tp the BEA contains the same information as IV. 
Indeed it is related to IV by means of functional differential equations like Eq. 1)6.10)1 . 
Using this relation we can conveniently write the SDE ()6.8)) as 

A (n) 

G a = -ja/3 + a °^'" an G a3aA ■ ■ -G an ^ ian . (6.11) 



Using (j6.10j) we obtain a differential equation for IV 
fir 

B — i \^ a/3a s ...a n n n ^ 19 \ 

oc/3 "r p ^0304 ' ' ' ^Qn-lQn- l^U.lZj 



dG Ql 3 

it 

which we can integrate 2 to obtain 

1 A^ n) 
IV = + ai F a " G aia2 ---G an _ ian , (6.13) 



the BEA at "one- vertex order". 

It is sometimes convenient to introduce an auxiliary effective action 

t[G,j] =r B -^j a pG a p (6.14) 

such that the physical propagator corresponds to the minimum of Y (cf. Eq. (|6.10)1 ). 



6.2 BEA for Local Interactions 



In the following we want to consider local interactions. For clarity we now write 
x (or momentum p) explicitly and use latin letters for the remaining indices. The 
standard procedure would be the insertion of the ansatz G~ h l (x,y) = —j a b( x ,y) + 
A-ab{%)3{ x — y) into Eq. ()6.11jl to obtain the SDE for the local gap A. Since the 
BEA Eq. (j6.13j) is related to the SDE (j6.11j) by differentiation with respect to G it 
is not clear that an effective action functional depending on A can be obtained by 
integration with respect to A. Instead wc want to follow the construction presented 
in 



9a | and start directly from the approximate BEA Eq. ()6.13|) . With 

g ab (x) = G ab (x, x) (6.15) 



2 Note that in our notation = 5 a ^5ps — SasSp^ 
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we have 

» 



i i r \W 

f = -Tr\ a G+-Tr(Gj) + J^^r L 9a^(x)---9a n - 1 a n (x). (6.16) 

For this relation it is essential that the interaction is strictly local. Furthermore, 
we can use the locality of the interaction to write (|6.11|) in the form of a local gap 
equation 

G ab( x 'V) = ~jab(x,y) + A ab (x)5(x - y). (6.17) 

We will evaluate the functional T[G] for G Q/ g corresponding to Eq. ()6.17j) . This is 
actually a restriction to a subspace of all possible G. However, locality tells us that 
the extremum (solution of the SDE) is contained in this subspace. 

Using j = — G~ l + A we find (up to a shift in the irrelevant constant and using 
A ab (x, y) = A ab (x)S(x - y)) 

T[g,A\ = ~Tr\n(-j + A) - ^ J A ab (x) 9a ,(x) 

+ ^Er L 9a 1 a 2 (x)---g an _ ian (x). (6.18) 



For the search of extrema of f it is actually convenient to treat A and g as inde- 
pendent variables. The extremum of F[g, A] then obeys 

gf^A] dT[g,A] 

-^A~ = °' -^g— = °- (6 ' 19) 

Evaluating the derivative with respect to A we recover the inverse of Eq. (|6.17|) for 

x = y, 

g ab {x) = (-j + A)- b 1 (x,x)=g[A(x)] (6.20) 

Inserting this functional relation into Eq. (J6.11|) leads to a gap equation for A. In 
case of a six-fermion interaction this takes, however, the form of a two-loop equation. 

For n-fermion interactions with n > 4 it is more appropriate to go the other way 
around and first take a derivative with respect to g. We obtain 

A ab (x) = ^ ab " 3 - Q " g aaa4 (x) ■ ■ ■ g an . ian (x) 

= A ab [g(x)], (6.21) 
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which is precisely the value of the gap in Eq. (jfi.lljh Inserting A[g] into ()6.18j) we 
find the effective action depending on g 

f\g] = ~Tr\n(-j + A[g})-±jA ab [g}(x)g ab (x) 

+ J Yl 9aia2 ^ " ' 9an - ian ^ " ^ 6 - 22 ^ 



Searching for an extremum yields 

*M = { (- j+ Ab])-'-*}^ = 0. (6.23) 

For j=f ^ Eq. (l6~2~3j) indeed corresponds to the SDE (IFTTTI) . i.e. 

g ab (x) = (-j+A[g])- b \x,x) (6.24) 

This will be our central gap equation. We should point out that possible extrema 
of T[g] corresponding to ^- = are not solutions of the gap equation ()6.1H) and 
should be discarded. Finally, we also have 

dt[g] df[g,A[g\] df[g,A[g\] + df[g,A[g\] dA[g] 



dg dg dg dA dg 

df[g,A[g}}dA[g] 
OA dg ' 



(6.25) 



Only as long as ^ is fulfilled we can conclude that a solution of ()6.23|) fulfills 
both extremum conditions (j6.19|) . 

Our procedure is quite powerful if Trln(— j + A) can be explicitly evaluated as 
a functional of A. Then T[g] allows not only a search for the extremum (discarding 
those with ^ = 0) but also a simple direct comparison of the relative free energy 
of different local extrema. This is crucial for the determination of the ground state 
in the case of several "competing gaps" . 

This "one-loop" form of the equation of motion but also of the effective action 
itself ()6.22|) is very close to what we would expect from MFT (cf. also the next 
section). In contrast to the standard SDE, which is an equation of motion, we can 
use Eq. ()6.22j) to compare the values for the effective action at different solutions of 
the equation of motion ()6.2Hj) . providing us with information about the stability. 

Nevertheless, we have to be careful when considering Eq. ()6.22j) at points which 
are not solutions of Eq. ()6.23|) . Going step by step through the procedure above, 
we find that if we are not at a solution of (j6.23J) we do not necessarily fulfill the 
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ansatz (jfi.l 7)1 . Therefore, at these points we are mathematically not allowed to insert 
the ansatz into Eq. (JHHSJ). So, strictly speaking ()6.22j) only gives the value of the 
effective action at the solution of the equation of motion 3 . Although, this is already 
more than we get from the standard SDE we would like to interpret (I6.22|) as a 
reasonable approximation in a small neighborhood of the solution to the equation of 
motion. Remembering g(x) = (ijj(x)ip(x)) it is suggestive to interpret g as a bosonic 
field. Eq. fl6.21|) gives the (non-linear) "Yukawa coupling" of g to the fermions, i.e. the 
relation between the gap and the bosonic field. The Tr In is the contribution from the 
fermionic loop in a background field g. The remaining terms can then be interpreted 
as the cost in energy to generate the background field g. This interpretation allows 
us to use ()6.22j) to calculate the mass and the couplings of the bosonic field g. 



6.3 Comparison with MFT 

From Chap. 0] we know that partial bosonization is not restricted to four-fermion 
interactions. In particular Eq. (|4.11|) provides us with the means to calculate a par- 
tially bosonized action for an arbitrary local multi-fermion interaction. Associating 
<t>ap{x) — (i J a{x)'4'i3(x)) , the partially bosonized form of Eq. ()6.2j) becomes 

Sirt[<f>,i>] = J2 m *X...anbn\^b 1 (x)---<f>anbM) (6.26) 

Jx n ? I 

j r Tl ~\ ~ I 

+ _ [^a 1 (x)ij bl (x)^ a2b2 (x)---(f) a „ b „( y x)+peTm. {1, ... ,-}J +C(^ 3 )|, 

<U = -^ + ^... an - (6.27) 

Where £ is a sum of terms which are symmetric in at least one pair of indices 
(cf. Sect. 14.1. 2j) . The condition ()6.27|) ensures that the partially bosonized action is 
equivalent to the original fermionic one. 

Neglecting the terms 0(ip 3 ) and performing the functional integral over the 
fermions provides us with the MF effective action: 

r MF M = ~'&M-J + ^])+/E m S..6n^«i6i---^, (6-28) 

9[<f>]aibi = Z [)| m< ai...b^ 2 b 2 ---<Panbn + perm, {l, ... ,|} 



3 An alternative would be to choose the gap A as the "bosonic field" . Inserting Eq. H6.17JI into 
Eq. (|6.13|l we could calculate a functional T[A]. However, as one can check there are two drawbacks. 
First, even for four-fermion interactions, T[A] is usually unbounded from below when considering 
A — > oo. Second, in the case of a large four-fermion coupling the "stable" solution of the field 
equation is usually a local maximum. 
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By construction this is a one- loop result. Moreover, it is strikingly similar to Eq. 
Ijfi.22j) . However, the coefficients differ. As discussed in Sect. l4.1^2"l the Fierz ambiguity 
is reflected by the presence of a nearly arbitrary E in Eq. ()6.27|) . Results usually 
depend on the choice of E. Of course, further considerations as e.g. the stability of the 
initial bosonic potential might reduce the freedom of E somewhat. But, sometimes 
this is not even enough to get qualitatively unambiguous results [^H EH- 

Eqs. (j6.18J) . ()6.22|) . ()6.23|) do not suffer from such an ambiguity since in the 
derivation of the SDE ()6.8|) the coefficients become antisymmetrized and symmetric 
terms drop out. In Sect. 14.31 we demonstrated that the inclusions of certain diagrams 
cures the Fierz ambiguity for four-fermion interactions and leads to the SD result. 
We believe that this holds for higher fermion interactions, too. Thus, we propose 
ClS db natural generalization of (j6.28|) . 

Finally, let us stress the similarity of both approaches by noting that g(x) = 
(i/j(x)i/j(x)) is exactly what we had in mind as a "mean field". 



6.4 Wave Function Renormalization 

In the partially bosonized language (cf. sect. 16. 3 j) it is possible to calculate a wave 
function renormalization for the bosons. Allowing not only for constant but also for 
a slightly varying <p (p small) 

(j)(x) = 0(0) + 5<pexp{ipx). (6.29) 

Using this it is still possible to perform the fermionic integral. Expanding in 
powers of the momentum up to the p 2 -term we can read off the wave function 
renormalization. 

The same can be done for Eq. ()6.18|) 

g(x) = g(0) + Sg exp(ipx). (6.30) 

Expanding again in powers of momentum we interpret the p 2 -term as the wave 
function renormalization for the boson corresponding to the field g. As in MFT the 
only contribution to the wave function renormalization comes from the Trln and 
therefore from a simple one-loop expression. 

Knowledge of the wave function renormalization together with the second deriva- 
tive of the effective action for constant fields allows us to compute the mass of the 
boson. 

Again, this calculation is unambiguous. This is in contrast to a calculation in 
the partially bosonized language where we again have the problems with the Fierz 
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ambiguity. On the other hand we do not want to hide the fact that by considering 
values of g which do not coincide with the solution of the SDE we have left the solid 
ground of a direct computation from the BEA given in Eq. ()6.16|) . 



6.5 Chiral Symmetry Breaking from a 3-Flavor 
Instanton Interaction 



In this section we want to use the method described above to study chiral symmetry 
breaking in an NJL-type model with a six-fermion interact ion mode l ling the Q CD- 
instanton interaction with three colors and three flavors 9^; lOOl : 101 ; 1021 : Hoi . The 



three flavor instanton vertex can be written in the following convenient form |4 

'S'illStl^] g / ( 'I ( l.ilVj/;;; 



-{R^L)\ (6.31) 



where X z are the Gell-Mann matrices corresponding to the SU (3) c color group and 
the brackets ( ) indicate contractions over color and spinor indices. 

The coupling constant ( can be calculated in terms of the gauge coupling. How- 
ever, it involves an IR divergent integral over the instanton size. Therefore, one 
needs to provide a physical cutoff mechanism. To avoid this difficulty we treat ( as 
a parameter. 

Inspection of ()6.3H) tells us that this interaction is U(l) anomalous with a resid- 
ual ^3-symmetry. This is important because we cannot restrict ourselves to real 
condensates from the start. 

In order to extract the interaction matrix A we have to antisymmetrize over 
flavor indices (a = 1 ... 3), color indices (i = 1 ... 3) Weyl spinor indices (a = 1, 2), 
chirality indices (x = 1, 2 = L, R) and the indices distinguishing between ip and ip 
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(s = l,2). 
A, 



x mi...me 



A 



m\...m§ 



i£aia2CL3^CL4CL5CLQ 8a\ 04 8&2 (25 8, 



a 3 a 6 

1 1 2 1 



P{A mi ... m6 }, (6.32) 
12 

^ [^1*4^2*5^326 

X [^Xi 1^X21^X31^X42^X52^X62 - ^xi2^X22^X32^X4l^X5l^X6l] 
x[5 Sl2 5 S2 25 S3 2^4l^5l^6l + ^il^ 2 1^3l^ 4 2^ 6 25 S6 2]. 

Here P denotes the sum over all 6! permutations of the multiindices mj = 
(cij, ij, ctj, Xj, Sj), j = 1 . . . 6, with minus signs appropriate for total antisymmetriza- 
tion. 

As a first example we consider a flavor singlet, color singlet scalar chiral bilinear 



(a = i^ L r R , = - 

9mn 9 'aia\s ,bj firt 
1 



(6.33) 



6 



5ab5ij5 a /3 



o-(S xl 8 T2 S s2 Sti - 8 x2 8 rl 8 sl 8 t2 ) - o-*(8 x2 8 rl 8 s2 8 n - 8 xl 8 T2 8 s i8 t2 ) 



We evaluate 



A 



(n) 

mm; m 3 nms mg 



" ' 9m%m*,9m$mQ 



X (n) _ x (n) + X (n) 



(6.34) 



9m,2m 5 9' 



m2m 5 ym 3 m 6 9m2rri( i 9m 3 m 5 



and 



U 



— A Tnl n I A mi .. -m6 
2 li»J mn9mn i 



9mim2 9m,3m4,9msms 



' r) .^mi...m 6 9mim29rn 3 rn i 9rri 5 9rn e , 



(6.35) 



where we have used the fact that A is symmetric under permutations of the three 
ipip bilinears. Exploiting the flavor, spin and color structure for (j6.34j) yields 



AM, 



— — (8 a b8ij8 al 3 

y 



a 2 (8 x i8 T2 8 s2 8 t i - 8 x2 8 Tl 8 s i8 t2 ) 



-a* (8 x2 8 T i8 s2 8ti - 8 x1 8 t2 8 s i8 t2 ) 



(6.36) 
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and 

20 

U[a] = -av 3 + v* 3 ). (6.37) 

Evaluating for a a constant in space and pulling out a volume factor we obtain the 
effective potential and the relation between condensate and fermion mass, 

f[<j] = — 36t>4 J dx x[\n(x + \m\1)} + U[a], 

m a = ^Ca 2 + m°. (6.38) 

Here m° is a current quark mass which we take to be equal for all quarks. The 
integral in (j6.38|) is, of course, divergent. Our UV regularization is simply to cut it 
off at A 2 . Measuring all quantities in units of A we can put A = 1. 



6.5.1 The chiral limit m a = 

Let us now look at the field equation or equivalently search for extrema of T[a]. 
Since d/ ^J^ ^ for all a ^ we do not need to worry for non-trivial solutions to be 
spurious. In addition a = is always a solution in the chiral limit. 

Inspection of f [a] tells us that it is invariant under the combined operation 
( — ► — £, a — > —a. This allows us to restrict our analysis to positive £. 

In the chiral limit it is useful to parametrize 

a = \a\ exp(ia). (6.39) 

Form this one finds 

40 

r[|<7|, a] = -CM 3 cos(3a) + /(H), (6.40) 

where / is a function determined by the integral in Eq. ()6.38)1 . We can see that the 
only a-dependence comes from cos(3a) which is the explicit manifestation of the 
i^-symmetry. 

It is clear that extrema can only occur at a = y, n G Z. Using the ^3-symmetry 
we can restrict ourselves to a = 0, n or restrict ourselves simply to real a. 

Taking all this into account we find up to three solutions (cf . Fig. 16.1)1 . As already 
mentioned a = is a solution for all values of the coupling. Going to larger couplings 
we encounter a point Ccrit where we have two solutions. For even larger couplings 
there are three solutions = a < ai < a 2 . We know that T[ui] > T[o"o = 0] 
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Figure 6.1: Plot of the BEA for various values of the coupling constant increasing 
from the topmost line ( = 3000 to the lowest line ( = 4200. The second line 
(long dashed) is for ( crit rs 3350, the critical coupling for the onset of non- vanishing 
solutions. The third is for ( = 3600 while the fourth (short dashed) is the for the 
onset of SSB (ssb ~ 3900. The horizontal line indicates the value of T[a] at the 
trivial solution o = 0. 

therefore <j\ is not the stable solution. As can be seen from Figs. 16.11 16.21 there is a 
Ccrit < C < Cssb where there exist non-trivial solutions to the SDE but there is still 
no SSB because r[cr 2 ] > r[a = 0]. We point out that in order to calculate Cssb we 
need to know the value of T, i.e. information beyond the SDE. 

In Fig. l6.2l we have plotted the mass gap versus the six-fermion coupling strength. 
Looking at Fig. 16.21 we observe a first order phase transition. From Eq. ()6.22|) we 
actually expect this quite generically as long as we have only n-fermion interactions 
with n > 6. However, this might also be an artifact of the "one-vertex" approxima- 
tion. 

Finally, we would like to remark that in general it is not enough to simply inte- 
grate the SDE with respect to the gap A reconstruct the effective action functional. 
This functional, let us call it T[A] is in general not equal to the the BEA T[G], not 
even at solutions of the SDE. Indeed, as can be seen from Fig. I6.2[ the results for 
physical quantities like the effective fermion mass can differ. The underlying rea- 
son for this is that the gap A is not the correct integration variable. The SDE is 
obtained by a G-derivative of the BEA functional T[G}. Therefore, in order to recon- 
struct T[G], we have to integrate with respect to G. As can be seen from Eq. (|6.21j) 
A is, in general, not even a linear function of G. Simple integration with respect 
to A therefore neglects the Jacobi matrix, which is a non-trivial function of A for 
interactions more complicated than a four-fermion interaction. 
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Figure 6.2: Fermion mass \m\ a (gap) versus the strength of the six-fermion inter- 
action ( in the chiral limit m° = 0. The thick line corresponds to the solution 
with smallest action. The dashed line is the largest non-trivial solution. Finally, the 
thin line is obtained by minimizing an "effective action" T[A] obtained by direct 
integration of the SDE with respect to the gap A. We find three special couplings, 
Ccrit for the onset of non-trivial solutions to the SDE, £ssb for the onset of SSB, 
i.e. a non-trivial solution has lower action than the trivial solution and ( SSB where 
the lowest extremum of T[A] becomes non-trivial. We point out that these three 
"critical" couplings differ. Moreover, to calculate Cssb we need to know the effective 
action. In our approximation we obtain a first order phase transition. 

6.5.2 Non- vanishing current quark masses m® ^ 

The non-vanishing current quark mass explicitly breaks the residual i? 3 -symmetry. 
r[|er|, a] does no longer depend on cos(3a) only, and we have to look at the complete 
complex cr-plane for possible extrema. 

Moreover, for m° a = 0, T[a] is completely symmetric under ( — > — (, a — > — a. 
Therefore, we could restrict ourselves to ( > 0. For ^ we need to add the 
transformation m° — > -m". We can still restrict ourselves to positive ( but we need 
to consider both positive and negative m° a . 

In the case of m° a ^ we still encounter an extremum of f [a] at a — 0. However, 
in this case it is not a solution of the SDE. It is a spurious solution due to = 0. 
The difference to the chiral limit is that the derivative of the effective potential f [a] 
now has only a simple zero while in the chiral limit it is a twofold zero. After dividing 
the field equation by a simple zero remains, giving a solution of the SDE in 
the chiral limit. 

Although, chiral symmetry is now broken explicitly we can still observe a first- 
order phase transition signaled by a jump in the fermion mass. The critical coupling 
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Figure 6.3: Dependence of the critical coupling on the current quark mass m° a (thick 
solid line). Jump in the fermion mass at the phase transition (thick dashed line). 
We observe that there exists a critical m° « 0.076 above which there is no phase 
transition. 



for the phase transition depends on m° as depicted in Fig. I6.3I The critical line ends 
at m° = m|J crit « 0.076, i.e. for m° > m" crit we have no first order phase transition 
in our approximation. 



6.6 Color-Octet Condensation 



In the last section we have considered only one direction in the space of all possible 
g resulting in a phase diagram for chiral symmetry breaking. Let us now consider 
the more general case where we also allow for a non-vanishing expectation value in 
the color-octet channel, more explicitly in the color-flavor locking direction. 



Qmn 9aict X s,bj /3rt 

1 

6 
1 

12 



(6.41) 

cr(S x iS T2 S s2 S t i - S x2 S T i5 s i5 t2 ) - cr*(S x2 S T i5 s2 5 tl - 5 xl 5 T2 5 sl 5 t2 
x{^xi^r2S s2 S t i - S x2 S Tl S s iS t2 ) - x*{^x^TiS S 2^n - S x iS T2 S sl 8 t2 ) 



Following the outline of the previous section we obtain 



AM, 



-C 



5 2 i a t 2\- - ^ 



1296 

(X <-> r, a a*, x ->• X 



v 432 
5 x i5 T2 5 s2 Sti - 5 x2 5 Tl 8 s i5 t2 



X(6cr + x)SaiSbj 
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and 



U[a, X ]=C 



20 



-a 



27 



243 



X 



c.c. 



(6.42) 



f[<7, x) = — 4t> 4 / dx x[8 ln(x + |m| 2 ) + ln(x + \m\ 

J X 



)} + U[*,x) 



,C(288a 2 + 6a X -7x 2 )+m° 



1296 

^C(18a 2 -3a X - X 2 )+m°. 



(6.43) 



In the chiral limit every point on the line x = ~ 6a (m a = and m x = 0) has 
the same value of T, and both derivatives with respect to a and x vanish. But, in 
this direction the derivative of A with respect to the condensates vanishes, too, and 
is indeed the null function in this direction. Therefore, on this line only the point 
(a, x) — (0, 0) is a true solution to the SDE. 

Restricting both a and x t° be real we have not found a solution with x / 0. 
Thus, we have not identified a solution which breaks color symmetry but not parity. 

For the most general case of complex a and x things are considerably more 
difficult since we now have to search for an extremum of a potential which depends 
on four real parameters. We checked several values of the coupling constant. So far 
we have not found a solution which has lower action than the lowest one with x = 0. 

Still, we would like to point out that the potential is unbounded from below. In 
various directions including those wi th y 7^ 0, T — > —00. Therefore, a physical cutoff 
mechanism like the one discussed in |49( or a better approximation which makes the 
potential bounded from below may provide additional solutions. 



Chapter 7 



Outlook: Quest for a 
Renormalizable Standard Model 



The "Standard Model" (SM) of particle physics is probably one of the most widely 
studied physical theories. It describes a wide range of physical situations with a 
satisfactory amount of accuracy. Yet, there are still open questions. In particular, 
we want to concentrate on "Renormalizability" and the "Hierarchy Problem", as 
those are problems tightly connected to the existence of an elementary scalar boson 
- the Higgs - in the SM. This returns us to the speculation of a composite Higgs 
which we mentioned in the introduction. 

Before going into details, let us briefly outline those two problems. 

Renormalizability: If asked whether the SM is renormalizable many physicists 1 
would answer this question positively. So, why search for something we have already 
found? Well, while this answer is certainly correct for most practical purposes, i.e. 
when the UV-cutoff scale is of a reasonable size A < ITeV, it is not certain that this 
is so if we send the cutoff to infinity. 

To get a grasp of the problem let us look at the example of 4 -theory (a model 
for the Higgs potential). Roughly speaking renormalizability means that the physics 
at short distances does not really matter, and therefore we can send the cutoff to 
infinity without changing the results. Yet, a straightforward calculation gives the 
(1-loop) running four-boson coupling as follows, 



°This is work in progress in collaboration with Holger Gies and Christof Wetterich. 
X A (not representative) survey in the "Graduiertenkolleg: Physical Systems with many Degrees 
of Freedom" of the university of Heidelberg resulted in 70% positive answers. 
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where q 2 is the (Euclidean) momentum scale and A is the coupling at the scale fi 2 . 
From this we can read off, that the coupling constant grows with increasing momen- 
tum scale. It is plausible that this does not fit into our picture that long distance 
physics (small momenta) decouples from short distance physics (large momenta), as 
the coupling is strong in the latter regime. 

This gives us an intuitive understanding that problems might arise when we 
want to send the UV cutoff to infinity. Those theories are not renormalizable in a 
strict sense. Requiring that a "fundamental" theory should be renormalizable such 
theories cannot be "fundamental", i.e. valid at all scales. They are valid only as 
effective theories up to a certain cutoff scale A where "new" physics will set in. 



This [H EES HQ and a similar problem in QED (Landau pole!) [Toil: flOQ] . 



hint to serious trouble in the Higgs respectively ?7(l)-sector of the SM. This is often 
referred to as "triviality" because when starting with some bare interactions and 
sending the cutoff to infinity, the renormalized coupling will be strictly zero. 

Hierarchy Problem: In contrast to chiral fermions, where chiral symmetry 
prevents the fermion mass from acquiring large quantum correction, the mass of a 
scalar boson is not protected against such corrections. Thus, assuming that the RG 
flow of the SM is "released" at, say, the GUT 2 scale, enormous fine-tuning of the 
scalar initial conditions is required to separate the electroweak scale from the GUT 
scale by many orders of magnitude. 

To illustrate this let us do a very simplified calculation. In a crude approximation 
the Higgs mass runs as follows (c 2 = const. = 0(1)), 

mWi)-mW2)=c 2 (q!-q 2 2 ). (7.2) 

Defining the dimensionless Higgs mass fi 2 (q 2 ) = mH q [ q ~ this can be rewritten as, 

A£) = c 2 - 4 (c 2 + ■ (7-3) 

So far so good, but let us now consider two very different scales, e.g. the GUT 
scale q\ = Mq UT ~ 10 30 (GeV) 2 and the electroweak symmetry breaking scale 
q 2 ~ 10 4 (GeV) 2 . Inserting values for the Higgs mass at the electroweak scale we 
find that the brackets on the RHS are of order 0(1). Hence, 

/i 2 (M P ) = c 2 - e>(l(T 26 ), (7.4) 

the dimensionless mass at the GUT scale must be fine-tuned to c 2 by an incredible 
amount. Although not excluded, such an amount of fine-tuning seems unnatural. 



2 Grand Unified Theory. 



7.1. UV Fixed Points and Renormalizability 
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7.1 UV Fixed Points and Renormalizability 

Let us now address the problem of renormalizability on a more formal level. 



7.1.1 Non-perturbative Renormalizability 

Commonly, theories are considered to be renormalizable if they have a (small cou- 
pling) perturbative expansion and we can absorb the infinities by a finite number of 
counterterms. In consequence, they only have a finite number of couplings (masses 
etc.) since it requires physical information to fix a counterterm unambiguously. For 
instance, the value of a coupling or mass at some scale cannot be determined in 
the theory itself, but must be measured. All other quantities can be calculated from 
these couplings. From this renormalizable theories derive their predictive power. In 
particular, the renormalization procedure allows us to take the continuum limit, 
i.e. to send the UV cutoff to infinity. In this sense we can consider renormalizable 
theories as fundamental. A prominent example of such a theory is QCD. 

It comes as no surprise that not all theories have this property of "perturba- 
tive renormalizability" (PR). Gravity is probably the most well known example, 
but there are many more: as theories which contain a coupling of negative mass 
dimension are usually not PR. The NJL model (jl.ljl has a coupling ~ (mass)~ 2 
and consequently is not PR, either. Moreover, naively renormalizable theories with 
dimensionless coupling constants may also face difficulties spoiling PR, as discussed 
at the beginning of this chapter. Usually, theories which are not PR are thought to 
be effective field theories valid only up to a finite UV scale A. 

Yet, PR is not the end of the game. Using a little bit of imagination it is obvious 



that a logical alternati ve is a "non-perturbative" renormalizable (NPR) jllOt llll 



1121 list 114 ; 11.4 llfij | theory. So far so good, but does such a thing exist and what 



is it? Fortunately, the answer to the first question is yes. Among the examples are 
various versions of t he Gross-N eveu model. More recently, it has been proposed that 
gravity is NPR, too 117; 118l |. But now, let us find out what hides behind NPR 3 . 



The space of all possible action functionals can be parametrized by an infinite 
number of dimensionless couplings (if necessary we use a suitable rescaling with the 
scale k). As discussed in Sect. 12.1.21 the RG describes a trajectory in this space (it 
is parametrized by the scale t = \n(k/ka)), 

d t gl = (3\g\g 2 ,...). (7.5) 

Starting point for the construction of a NPR is the existence of a non-Gaussian fixed 
point (FP) g* = (gl, g\, . . .) with at least one g\ 7^ 0, in the RG flow, 

P\gl,gl...) = VI (7.6) 



In this brief description of NPR we follow [11 
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In the setting of statistical physics such a FP is exactly what we would associate 
with a second order phase transition. The dimensionless constants do not change 
with the scale (typically the scale can be associated with the difference from the 
critical temperature ~ \T — T c \), accordingly the dimensionful quantities are simply 
powers of the scale 4 . This gives the typical power laws of the critical behavior near 
a second order phase transition. Moreover, we would like to remark that in this 
context a PR is just the special case of a Gaussian, i.e. g± = 0, FP. 

In the vicinity of this FP we can linearize the RG equations, 

<hg (7.7) 



where 

dgi 



(7.8) 



a* 

is the Jacobi matrix of the /9-functions. This is a linear differential equation and the 
general solution reads, 

9l = 9l + T, ClV '{j) 61 > ( 7 - 9 ) 

where the V 1 are right eigenvectors of B with corresponding eigenvalues 5 — 9 / , 
BV 1 = — O^V 7 , and the C 1 are constants determined by the initial conditions. Set- 
ting C 1 = if the corresponding eigenvalue 6 7 < we can safely take the UV limit 
k — > oo. In other words this gives us an RG trajectory which ends in g* for k — > oo. 
Any such trajectory defines a theory with a meaningful UV limit. The space of all 
such trajectories is a submanifold S of dimensionality A, given by the number of 
eigenvalues Q 1 > 0. We can specify a trajectory in this space by giving the values of 
integration constants C 1 . Thus our theory has A "renormalizable couplings", which 
have to be taken from experiment. In particular, A must be finite (and preferably 
small) in order for our theory to have predictive power. 

Stated differently, all trajectories in the submanifold S are attracted toward g* 
for increasing scale k, consequently, for decreasing k they are repelled. Therefore, 



4 It is often useful to include the wave function renormalization into the couplings. This modifies 
the naive power laws by anomalous dimensions. 

5 The Q 1 are not necessarily real as B is not necessarily symmetric. For simplicity let us pretend 
that they are real, anyway. The general case can be recovered by replacing, 9 > — > Re{Q) > 
etc. in our argumentation. 

6 For Q 1 = it depends on the details if the UV limit is finite or not. Correspondingly, it might 
or might not be necessary to set C 1 = 0. Prominent examples for both cases are the </> 4 -theory 
where the UV limit is not finite, and QCD where the limit is finite and which is therefore strictly 
renormalizable. 
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the A parameters specifying the trajectory are the relevant parameters to describe 
physics in the fixed point regime. The remaining "non-renormalizable couplings" 
corresponding to Q 1 < are irrelevant in the sense that starting with a finite value 
at some scale k < oo they are attracted toward the sub manifold S for decreasing k. 
Thus it does not really matter if we give them some finite value at a large scale 7 k. 

Finally, let us give a naive argument why the submanifold S should be finite di- 
mensional. To obtain the dimensionless coupling we have to rescale the dimensionful 
coupling constant G % by an appropriate power of k, g l = k~ d 'G. The /9-function has 
now the form f3 l = —d l g l + • • • where the dots denote the loop corrections. In the 
Jacobi matrix this gives a contribution B l i = —d l 5 l i + • • • . The eigenvalues get a 
contribution Q l — d l + • ■ ■ . Constructing a higher order operator, we usually add 
derivatives or field operators, thus decreasing d % . Therefore, we should only have a 
finite number of 0* < 0. Yet, this argument is based on the assumption that the loop 
corrections are small. Hence, it might be not valid in the non-perturbative regime 
with strong coupling and large loop corrections. Still, it seems reasonable to first 
look for a possible FP and establish that it is not an artifact of the approximation, 
and then worry about the dimensionality of S. 

7.1.2 A Toy Model 

Let us illustrate this idea for a simple NJL model in <i-dimensions, a four-fermion 
interaction with dimensionless (positive) coupling constant A = k d ~ 2 \, and the flow 
equation (which, for the moment, is assumed to be exact), 

d t X = {d-2)\-C\ 2 . (7.10) 

For d > 2 the Gaussian FP A = is UV repulsive, i.e. A is an eigenvector with 
eigenvalue — O = (d — 2) > 0. The model is not PR for a non-vanishing A. Moreover, 
for any value A < ^ we approach the Gaussian FP in the IR, A is "irrelevant" . By 
contrast, for A > ^ the coupling A grows and eventually diverges. This typically 
signals some kind of symmetry breaking, A is not so irrelevant any more. Looking a 
little bit more closely we notice that A c = is a FP. The eigenvector (A — A c ) has 
an eigenvalue —0' = —(d — 2) < 0, i.e. it is a relevant parameter in the fixed point 
regime. On the other hand, the FP is UV attractive and we have a meaningful UV 
limit even for a (small) non- vanishing (A — A c ). The theory with this parameter is 
NPR. 

Finally, let us remark that the ad hoc flow equation of our toy model resem- 
bled a typical flow equation for an NJL-type model in a truncation to four-fermion 

7 We will see below that this statement should be taken with some care. More precisely, this is 
true only if we are close enough to the fixed point. 
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interactions, e.g. with C = 4(iV — 2)t> 3 Zj and d = 3 the flow equation for the 
Gross-Neveu model in three dimensions. 



7.1.3 Manifestation of the Hierarchy Problem 

In this setting we also have the tools available for discussing the hierarchy problem. 
So far we have already seen that a renormalizable coupling is relevant in the infrared, 
i.e. starting with a small deviation from the fixed point in the far UV this deviation 
soon grows large. This is quantified in Eq. (|7.9j) . Taking k to be some large UV 
scale k = A ~ Mp and G ~ 0(1), the deviation grows very fast with (4) e when k 
is lowered. Turning the argument around, we have to fine tune the initial conditions 
(choose the initial C very small) at the scale A in order to achieve a value comparable 
to some scale k -C A for the coupling at the scale k. 

From this it may seem that renormalizability and a solution to the hierarchy 
problem more or less exclude each other. However, there is a way out, we simply 
need to make small or zero. At first this may sound simply like another type of 
fine-tuning, but it is not necessarily so, as the eigenvalues are a prediction of our 
fixed point solution and not a parameter. To understand this we can once more 
look at our simple toy model. For d = 2 we have exactly the case of a vanishing 
eigenvalue. Solving the flow equation (j7.10j) . we find 

A(t = ln(fc))= A ° (7.11) 
1 + CA In | 

where Ao = A (to = ln(A). As A depends only logarithmic on the scale we can have 
very different scales without fine-tuning. 

Finally, let us remark, that Eq. (j7.11j) also gives us an understanding of what 
can happen when = 0. Only for C > we approach the fixed point A* = in the 
far ultraviolet. Only in this case our simple toy model is renormalizable. Moreover, 
our d = 2 flow equation for our toy model has the same form as the lowest order 
terms of the d = 4 flow in QED and QCD around their Gaussian fixed points. In 
this language QCD has C > and is renormalizable while QED has C < and is 
not strictly renormalizable. 



7.2 One more NJL Model 



Inspired by our toy model let us once again investigate an NJL model but now with 
Nf fermion species and an S'f/(N c )-gauge interaction. 
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7.2.1 Truncation and Flow Equations 



A simple truncation including all possible pointlike the four-fermion interaction 
reads, 



T fe = / + \f^F, u + ~ [K (S-P) + A_ (V-A) + A + (V+A) (7.12) 

+A^ (S-P) N + K (S-P) N + A VA (V-A)r 



with the covariant derivative, 



(7.13) 



The color and flavor singlets are 



(S-P) = W) 2 - (V76<0) 2 , 

(V-A) = (^) 2 + (^7m75^) 2 , 
(V+A) = (^) 2 - (V^TsVO 2 , 



where color (z,j, ...) and flavor (a, 6, . . . ) indices are pairwise contracted, e.g., 
(■0V) = i^i^i)- The operators of non-trivial color or flavor structure are denoted 
by, 



(s-p)n = mrf - (wwfc) a = (^V;) 2 - G5?7bV£) 2 , 

(S-P) N = (f^ 6 ) 2 - (?W) 2 = @M) 2 - (^75^) 2 , 
(V-A) N = (V^^i) 2 + tyfttfd'jY = (v-a) n . 

The last equation holds because of a Dirac Fierz identity (cf. App. IB.3|) . The trun- 
cation ()7.12j) is symmetric under a simultaneous exchange of N c <-> Nf, A^ <-> A£, 
(. . . ) N «->(.. . ) N . However, it is not invariant under S , Z7(Nf)] J x S'J/ (N f ) R . We can ob- 
tain the part invariant under this additional symmetry by setting A CT = X c a = 0. The 
full action Eq. ()7.12|) is invariant only under the subgroup SU (Nf)v of simultaneous 
right and left handed rotations. 

Following along the lines of Sect. 13.31 in particular using Eq. ()3.16|) . the cal- 
culation of the flow equations is straightforward. Using the dimensionless coupling 
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constants A = k 2 \ we find, 



dtX-i 



d t X 



t^VA 



2X a - r2g 2 [%-^\ a + \l)v 4 



(7.14) 



8v/P' 4 [2N c N f A 2 - 2X a [A_ + 3A+ - 2(N f A= + N C A^)] + 3A C CT A^ , 

(12 + N 2 (3 + af) gW ( ™ )A + (2A_ - ^ (a_ - N C A VA ) „ 4 ) 



(F),4 



N f N c (A 2 _ + A 2 ) + A 2 - 2(N C + N f )A_A V A 



+A + (A CT + N C A^ + N f A^) + 2X\ A - \X 2 a - fX c X 
L (12 + N 2 (3 - «(6 + a))) gW^ + + A + 



8N 



- 8v/f )A - 3X 2 + - 2N c N f A_A + - 2A+(A_ + (N c + N f )A VA ) 

+A_(A CT + N C A C CT + mi) + A VA (A C CT + X{) + \{Xl + A C CT 2 + A^ 2 ) 

;(F),4 " 



S^/P 2N f A^ - 2A CT A VA - 2A_A C CT - 2N C A^A VA - 6A+A C CT - A CT A* 

N 2 . - 1 

4N r 



3# = -L (24 - N 2 (3 + a) 2 ) ,V/£ B) ' 4 + (2A f CT + 24, 2 (a + - , 4 ) 



(F),4 



2N C A CT 



f 2 



6A+A^. — A CT A^. 



2A(jAva — 2A_A^. — 2N f A^AvA 
gi- (24 - N 2 (3 - a(6 + a))) sVi?' 4 + ^A VA - ^(Ava - N C A_)«A 4 ) 
- 8t; 4 /! F) ' 4 [ - (N c + Nf)Ay A + 4A_A VA - |(N C A° 2 + N f A^ 2 ) - iA CT (A< + A f j] , 
where the gauge fixing parameter is denoted by a. 



7.2.2 Many Fixed Points but no Solution to the Hierarchy 
Problem — the Case of Vanishing Gauge Coupling 



To get a first insight we have numerically solved the FP equation, 

d t X = 0, 



(7.15) 



for vanishing gauge coupling g and N c = Nf = 3. It turns out that not only do we 
find a solution, but quite a few, 64 = 2 6 to be exact (one is always the Gaussian 
FP). Looking more closely, only 44 are real. 
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Let us, for the moment postpone the question, if the FP are physical or an artifact 
of the truncation. Having found so many FP we could become kind of greedy and 
ask if among all those 44 FP there is one which in addition to providing us with a 
renormalizable theory, could solve the hierarchy problem. However, calculating the 
eigenvalues of the stability matrix we find, that the largest eigenvalue max > 2 for 
all non-trivial fixed points. Or, stated more physically, one direction is at least as 
unstable as a scalar boson mass. 

Is this a very special property of our choice of Nf and N c ? No, but it is a property 
of our truncation to pointlike four-fermion interactions and vanishing gauge coupling. 
More precisely, in this truncation there is always an eigenvalue O = [d — 2). This 
can be seen by the following argument. In a four-fermion truncation we can write 
the flow equations in the form, 

d t X i ={d-2)X i + X k A i kl X h (7.16) 

where A 1 is a symmetric matrix i.e. A kl = A\ k . The stability matrix is then given 
by, 

B i3 = = (d - 2)5 V + 2X k A i kj . (7.17) 

Let A* be a solution of the fixed point equation, 

d t X-=(d-2)Xt + XiAl l X* = V*. (7.18) 
Acting with \\* on A* ^ we have, 

B H I a* A* = (d-2)A* + 2A*A} fc A* (7.19) 
= -(d-2)X* + 2(d-2)X* + 2X*A%Xt 
= -(d - 2)X* + 2((d - 2)A* + X*A) k X* k ) 
= ~(d-2)Xt, 

where we have used the fixed point equation ()7.18j) in the last step. 

This shows that A* itself is an eigenvector of the stability matrix with the eigen- 
value — (d — 2), hence O = (d — 2). Therefore, in this truncation there cannot be an 
infrared stable fixed point beside A = 0. 

7.2.3 Non- vanishing Gauge Coupling 

Having not found the desired properties for the eigenvalues of the fixed point, let us 
look if a non- vanishing gauge coupling can stabilize the fixed point. 
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Figure 7.1: Largest eigenvalue 0™ ax for the most "stable" fixed point (cf. Eq. ()7.20j0 
depending on the gauge coupling. At a realistic value g « 1 for the gauge coupling 
Qmax > 5 Giving us no solution to the hierarchy problem. 

A measure for the stability of a fixed point A*, x = 1 . . .44 is the size of its 
largest eigenvalue Q™ ax . When G™ 3 ^ is smaller the fixed point A* is more stable. 
Thus, we have searched for the smallest, 

Q™ ax = minQ™ ax . (7.20) 

X 

We have plotted this eigenvalue as a function of the gauge coupling in Fig. 17.11 To 
get an impression of a reasonable value of the gauge coupling, let us take the strong 
gauge coupling. At a scale of ~ lOOGeV we have g s « 1. At larger scales the gauge 
coupling is even smaller. As can be seen from Fig. 17.11 Q™ ax > 1.5 in this range, 
bringing us nowhere near the desired G ~ 0. 

7.3 The Future 

So far the results of our toy model can be summarized as follows. It looks as if there 
are many possibilities to have fixed points, but it seems difficult to find one which 
has very small positive eigenvalues G. But, as we have already seen in Sect. 17.2.21 
this may be an artifact of our truncation to four-fermion interactions. Therefore, a 
next step is certainly to enlarge the truncation. 

However, there is another interesting direction we can take. So far, in our simple 
approximations we have neglected the influence of the four-fermion interaction on 
the running of the gauge coupling. Instead we have treated the gauge coupling as 
coming from the outside. Yet, at a strong-coupling non-Gaussian fixed point the 
flow of the gauge-boson-fermion vertex is certainly influenced by the presence of the 
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Figure 7.2: Correction to the gauge-boson-fermion vertex (fermions solid with arrow, 
gauge boson wiggled) in the presence of a four-fermion interaction. 

four-fermion interaction. The lowest order contribution is depicted in Fig. 17.21 If we 
define the gauge coupling by this vertex, the running is modified by a contribution 

(7.21) 
with 

First of all this may have interesting consequences for properties like asymptotic 
freedom, because at a non-Gaussian fixed point the terms ~ g 2 X will dominate for 
very small gauge coupling. Therefore, the fixed point in A determines if the gauge 
interaction is asymptotically free or not 8 . 

Secondly, in Fig. 17.21 we have concentrated on the gauge-boson-fermion vertex. 
But, in non-abelian gauge theories alternative definitions of the gauge coupling are 
by the three- and four-gauge-boson vertices. Some thought reveals that those do not 
get a direct contribution from the four-fermion interaction. This raises the question 
how this can be reconciled with the fact that, at least naively, gauge invariance tells 
us that both definitions agree. These interesting possibilities are subject to future 
investigations. 



d t g 2 



-Sv 4 g 



Zi F) ' 4 (A ff - 2A_ + NfA^ - 2N f A VA ) + M 



8 As /3q > it is sufficient for asymptotic freedom that the sum of the four-fermion couplings 
in the brackets is positive. First numerical checks indicate that regardless of which of the 44 fixed 
points (N c = Nf = 3) we choose, asymptotic freedom is ensured. This even seems to apply to 
various other combinations of Nf, N c . 



Chapter 8 

Summary and Conclusions 



We encounter strongly interacting fermions in many situations, ranging from color 
superconductivity and chiral symmetry breaking (~ lOOMeV) to ordinary supercon- 
ductivity (~ meV). Typical features of such systems are the formation of bosonic 
bound states and spontaneous symmetry breaking (SSB). Non-perturbative tech- 
niques are essential as SSB cannot be described by (standard) perturbation theory 
in these systems (cf. Sect. 13 .2j) . Therefore, all methods discussed in the following 
correspond to non-perturbative resummations of perturbative diagrams. 

We have started this work with a comparison of various standard methods used 
for non-perturbative calculations in this kind of systems. We have calculated the 
critical coupling for the onset of spontaneous chiral symmetry breaking in a simple 
NJL model. In particular, we have used mean field theory (MFT), a renormalization 
group (RG) equation with a truncation to pointlike four-fermion interactions (from 
now on referred to as fermionic RG) and the lowest-order Schwinger- Dyson equation 

(SDE). The results for the critical coupling A™* and two different values of Ay can 
be found in Tabs. 15.11 15.21 Since the most characteristic features and problems of 
the different methods are most clearly seen when the couplings A CT and Ay are of 
similar size we concentrate on Tab. 15.21 

Both MFT and the lowest order SDE sum only over fermionic fluctuations in 
presence of a bosonic background. They include, in principle the same type of dia- 
grams, Fig. 13.11 The MFT result depends strongly on the choice of the background 
field. This "Fierz ambiguity" (FA) is expressed by the dependence on the unphys- 
ical parameter 7 in the tables. No such ambiguity appears in the SDE approach 
which therefore seems more reliable. We note that for a particular choice of 7 the 
MFT and the SD approaches give identical results - in our case 7 = 1/2. This has 
led to widespread belief that MFT and SD are equivalent if the basis for the Fierz 
ordering is appropriately chosen. However, this is not the case, as can be seen by 
calculating also the critical coupling where spontaneous symmetry breaking sets in 
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in the vector channel (in absence of other order parameters). There is again a value 
7 = —(A CT + Ay)/(2Ay) where MFT and SD give identical results, but it differs from 
7 = 1/2 as encountered in the scalar channel 1 . We conclude that there is no possible 
choice of 7 where both critical couplings for SSB in the scalar and vector channels 
are identical in the MFT and SD approaches. 

For a better understanding of the FA of MFT it is instructive to consider MFT 
on a more formal basis as a simple approximation, taking into account only the 
fermionic fluctuations, in a "partially bosonized" language. Partial bosonization is a 
powerful tool for an understanding of strongly interacting fermionic systems beyond 
the level of MFT or SDE. It allows us to treat the bosonic fluctuations in an explicit 
manner, treating them on equal footing with those of the elementary particles, and 
provides for a rather simple framework for the discussion of SSB. Most importantly, 
it permits the direct exploration of the ordered phase which is, in practice, almost 
inaccessible for the fermionic RG. Yet, already on the level of the classical action 
we can get a grasp of the origin of the FA. Partial bosonization is not unique. From 
one fermionic action we can obtain a whole family of equivalent bosonized actions, 
related to each other by redefinitions of the bosonic fields. 

In order to permit a simple comparison with the fermionic RG we have used a 
rather crude approximation for the purely bosonic sector by retaining only a mass 
term and neglecting bosonic interactions as well as the momentum dependence of 
the bosonic propagator. In this approximation the effect of the boson exchange 
between fermions does not go beyond pointlike fermionic interactions. Taking into 
account only the running of the Yukawa couplings (Fig. 14. 2|) in the bosonic RG of 
Sect. 14 .2] we observe already a very substantial improvement as compared to MFT. 
The dependence on 7 is greatly reduced and the numerical value of the critical 
coupling comes already close to the result of the fermionic RG. These features can 
be compared to the inclusion of higher loop effects in perturbation theory in particle 
physics: they often reduce the dependence of the results on unphysical parameters, 
like the choice of the renormalization scale. 

Using SDE in a partially bosonized language we fared even better. Again, MFT 
appears as the approximation which includes only the fermionic fluctuations in a 
bosonic background. Adding the mass-shift diagram (cf. Fig. 14 .3|) we recover the 
unambiguous result of the fermionic SDE - MFT sums only over a subset of the 
diagrams contained in the fermionic SDE. 

The mass-shift diagram in the partially bosonized language has only fermionic 
external legs. This has prompted us to look for similar purely fermionic contributions 
in our partially bosonized RG description. The box diagrams (Fig. 15. 1|) generate 
new four-fermion interactions and contribute to the same order as the mass-shift 
and vertex corrections. Those four-fermion interactions are included in the adapted 

1 Actually, 7 is negative and therefore outside the range of strict validity of MFT. 
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bosonic RG discussed in Sect. 15.21 Here the relation between the bosonic composite 
fields and the fermion bilinears becomes scale dependent. This formulation is well 
adapted to the basic idea of renormalization where only effective degrees of freedom 
at a certain scale k and their effective couplings should matter for physics associated 
to momenta q 2 <k 2 . The system should loose all memory of the detailed microscopic 
physics. In particular, the choice of an optimal bosonic field for the long distance 
physics should not involve the parameters of the microscopic theory, but rather 
the renormalized parameters at the scale k. In this formulation it has also become 
apparent that the distinction between "fundamental degrees of freedom" and "bound 
states" becomes a matter of scale |42T |. The adapted bosonic RG reproduces in our 
crude approximation the results of the fermionic RG. We argue that for precision 
estimates in the partially bosonized approach the "adaption" of the definition of the 
composite field seems mandatory. 

For an improvement of the truncation two possibilities come to mind immedi- 
ately. One is to enlarge the bosonic potential beyond the mass term, the other is to 
add kinetic terms for the bosons. Yet, our discussion of Sect. 15.31 shows that includ- 
ing only a bosonic potential without kinetic terms for the bosons does not help us 
in deciding which type of boson will condense in the SSB phase. The basic reason 
for this is, that the corresponding interactions in the fermionic language are still 
pointlike and can be Fierz transformed in many ways. Kinetic terms correspond to 
a momentum dependence of the interactions in the fermionic language, greatly re- 
ducing the freedom to Fierz transform. For this reason we have turned to the second 
possibility. To cut a long story short, for a consequent inclusion of all terms with 
up to two derivatives on the bosonic fields, we have to take into account not only 
the momentum dependence of the mass-shift diagram (cf. Fig. 13.1)1 . but in addition 
part of the momentum dependence of the vertex correction and the box diagrams 
(cf. Figs. El EH) - 

Comparing the numerical values (cf. Tabs. 15. 1[ 15.2)1 with those of the pointlike 
approximations, we find that the inclusion of the kinetic terms affects the critical 
coupling on a 10% level. Moreover, the contributions from the different diagrams 
are of similar size, confirming once more that an adaption of the flow is necessary to 
achieve a high precision. The 7 dependence is small, as we would expect for a sys- 
tematic enlargement of the truncation. Nevertheless, it does not completely vanish 
since we have been forced to make approximations when including the momentum 
dependence of the vertex corrections and box diagrams. Moreover, if the "right" 
FT is known, inclusion of the pointlike contributions seems already sufficient for 
reducing the 7 dependence (cf. the second to last line in Tab. 15.2)1 . In view of the 
high algebraic complexity involved with including the full momentum dependence 
we would like to suggest this approximation for more moderate demands on accu- 
racy. Finally, let us point out the following important feature (shared by all of the 
three considered approximations which include kinetic terms): with kinetic terms for 
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the bosons we can decide which boson will condense. For values X a (slightly) larger 
than A CT , only the renormalized mass of the scalar boson turns negative, while the 
renormalized masses of the vector and axial vector bosons remain positive. This 
confirms that we have a phase where chiral symmetry (and nothing else) is broken. 

In summary, the FA can be used as one possible test of approximation errors for 
some of the methods proposed to deal with strongly interacting fermionic systems. 
The spread of the results within an acceptable range of the unphysical parameter 7 
should be considered as a lower bound for the systematic uncertainty within a given 
approximation. We find that MFT can have a very substantial ambiguity which 
should then be reduced by systematic improvements. On the other hand, the FA is, 
of course, not the only source of error. Several methods such as SDE or the fermionic 
RG have no such ambiguity by construction. This holds similarly for the adapted 
bosonic RG (with or without kinetic terms) which is constructed to reduce the Fierz 
ambiguity. Hence, in this cases the uncertainty due to the Fierz ambiguity is small 
compared to other uncertainties which can roughly be estimated from the spread of 
the results between the different approximations. However, improving the truncation 
we should keep an eye on the FA as it is likely to increase, when our "improvement" 
neglects something essential. Comparing the RG and SD approaches the adapted 
bosonic RG sums over a larger class of diagrams and therefore seems more reliable. 
Moreover, the RG accounts for the fact that physics at the scale k should involve only 
renormalized parameters at this scale, while the SDE involves "bare" couplings. We 
think that with the adapted bosonic RG we have reached a promising starting point 
for future investigations along the lines of 0; • In particular, a more elaborate 
bosonic potential would allow us to proceed into the SSB phase. 

Removing the FA in a partially bosonized setting has turned out to be quite 
tedious. For a first investigation of the phase diagram, lowest order SDE seem to 
be a viable alternative which allow for a description of the SSB phase without 
auxiliary fields. However, at first sight there is one drawback: SDE determine only 
the derivatives of the action, not the action itself, making it difficult to compare 
different solutions with regard to stability. This shortcoming is cured by the use of 
the 2PI effective action, or its simplified form at vanishing fermionic sources, the 
bosonic effective action (BEA). In this context the SDE is the field equation of the 
BEA. 

Integrating the lowest-order SDE for a multi-fermion interaction we have ob- 
tained the bosonic effective action at "one-vertex" level. Within this approximation 
we computed a simple "one-loop" expression for the BEA and the SDE. In this form 
the BEA appears very similar to MFT but does not suffer from its ambiguities and, 
as we have already seen, sums over a larger class of diagrams. 

We have applied the BEA at one- vertex level to a six-fermion interaction resem- 
bling the instanton vertex for three colors and flavors. For vanishing current quark 
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masses we find a first-order phase transition. It turns out that the value Ccrit for 
the onset of non-trivial solutions for the SDE is not necessarily equal to the value 
Cssb for the onset of SSB. To calculate the latter we need the value of the effective 
action in addition to the solution of the SDE. For nonzero current quark masses and 
positive a coupling constant we find a phase transition only for current quark masses 
below a m° cvit w 0.076. At this point the gap in the effective mass vanishes and we 
expect a second-order phase transition. We have also searched for solutions with a 
non-vanishing color-octet condensate. Although they are not excluded by symme- 
try we have not found a stable solution of the field equations with non-vanishing 
octet condensate \ ^ 0. Nevertheless, we would like to point out that in our simple 
approximation the BEA is unbounded from below in several directions including 
some with x ^ 0. Hence, some sort of cutoff mechanism which renders the potential 
bounded might yield additional solutions. Therefore, a more detailed investigation 
of the instanton interaction in a color-octet background would be of great value. 

Leaving aside the more technical aspects of this work, we have turned to an 
intriguing speculation, namely that the Higgs is not an elementary scalar boson, but 
a bound state of fermions, more precisely a top-antitop bound state BBS 3- To 
motivate a thorough investigation we have briefly reviewed the "Hierarchy Problem" 
and the "Triviality" of the 4 -theory which cause trouble for the Standard Model 
(SM) with its elementary scalar boson. In particular, the latter problem may prevent 
the SM from being renormalizable in a strict sense. In a model with a bound state 
Higgs, renormalizability may be provided by a non-Gaussian fixed point in the RG. A 
toy model with Nf fermions interacting via a four-fermion interaction and a S77(N C ) 
gauge interaction implied that there might be plenty of those. On the other hand 
it became apparent that a solution to the hierarchy problem is beyond a simple 
approximation to a pointlike four-fermion interaction, as in this setting the required 
fine-tuning to achieve a separation of scales is more or less as bad as for an elementary 
scalar. However, already in this simple truncation another interesting effect turned 
up: at a non-Gaussian fixed point, the flow of the gauge coupling might be crucially 
influenced by the four-fermion-interactions. In particular, this aspect seems to be 
interesting for future investigations. 

It may well be that there are fundamental scalar fields, most candidates for 
unifying theories have plenty of them, but so far not even one has been detected 
with certainty. Hence, it might also be that nature has simulated once more an 
elementary scalar with a bound state. Upcoming experiments (e.g. at Tevatron and 
LHC) will help deciding this issue. But, independent of this, bound states of fermions 
are still abundant in nature and we hope that our formalism with scale-dependent 
degrees of freedom may be of help in understanding some of them. 



Appendix A 



Conventions, Abbreviations and 
Symbols 

A.l Conventions 



• We use Euclidean conventions in flat spacetime, i.e. the metric is the d dimen- 
sional unit matrix. 



• Greek indices /x, v . . . — 0, . . . , d denote spacetime indices. 



• Latin indices a,b . . . = 1, . . . , N f are flavor indices, = 1, . . . , N c color 

indices. 



• Our conventions for the Fourier transform are 




(A.l) 



I 



d d q 



■ij}(q) exp(ipx), ip(x) 



I 



d q — 



ip(q) exp(-ipx). 



(2%) d 



(2%) d 
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A. 2 Mathematical Symbols 





Proportional to 




approximately equal to 


® 


Tensor product 


t 


Hermitian conjugation 


C 


Charge Conjugation (operator) 


P 


Parity (operator) 


R 


Reflection (operator) in Euclidean spacetime 


T 


Time Reversal (operator) 



A. 3 Abbreviations 



Ad. 


Adapted 


BBS 


Bosonic Bound State 


BEA 


Bosonic Effective Action 


cf. 


confer 


Chap. 


Chapter 


ERGE 


Exact Renormalization Group Equation(s) 


Eq. 


Equation 


FA 


Fierz Ambiguity 


FP 


Fixed Point 


FT 


Fierz Transformation 


IR 


Infrared 


LHS 


Left Hand Side 


LPA 


Local Potential Approximation 


MF 


Mean Field 


MFT 


Mean Field Theory 


RG 


Renormalization Group 


RHS 


Right Hand Side 


s. 


see 


SD 


Schwinger-Dyson 


SDE 


Schwinger-Dyson Equation(s) 


Sect. 


Section 


SSB 


Spontaneous Symmetry Breaking 


uv 


Ultraviolet 


WFR 


Wave Function Renormalization(s) 



Appendix B 



Fermion Conventions, Fierz 
Identities 



B.l Dirac- Algebra in 4 Dimensions 

Throughout this work we use an Euclidean metric, g^ v = 5^ and 
\x\ 2 = + x\ + ■ • ■ + x 2 d _ x . With the exception of one tiny excursion in Chap. 03 
the number of spacetime dimensions will be d = 4. 

Accordingly the Dirac-algebra is 

= 7Y + 7Y = 2n, (B.i) 

7 5 = tSVt , 

^ = ^f,-f\ = ^(7V-7V). 

We use a chiral basis ^ = ( ^ L J , ^ = (ijj R , tp L ) , with the projection operators 

\ yr J 

Pl,k — on the chiral components. An explicit representation is then given by, 

rA° "TV "A (b.2) 



with cr M = (zl,er l ). a* are the standard Pauli-matrices 

•'-(Si)- '-(J-J- (B - 3) 

and 1 is here the 2x2 unit matrix. 
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Using Eq. (|B.1J) it is quite easy to derive relations to simplify expressions con- 
taining several 7 m atrice s, e.g. ^^ v ^ = —27^. This can be automated and we use 
the Tracer-package 119l | for Mathematica to do this. 



B.2 Dirac- Algebra in 3 Dimensions 

In Chap. El we use the Gross-Neveu model in 3 dimensions to demonstrate a short- 
coming of partial bosonization. The Dirac-algebra in odd dimensions is somewhat 
different from the case of even dimensions. 

Nevertheless, it is quite easy to find an explicit realization as the Pauli matrices 
already fulfill the requirements for the Dirac-algebra, 

{a\ a j } = 25 ij . (B.4) 

Consequently, we can use spinors with only two components. Since we do not need 
any more subtle properties, let us leave it at that and return to the case of four 
dimensions. 



B.3 Fierz Identities 



Defining Os = 1, Op = 7 
following Fierz identities 



r> O v = 7 M , A = 7 M 7 5 and T = we obtain the 



fixity Ma) = J2 C M^ a Oy^ b )^ c O Y ij d ) 



;b.s) 



with 



C 



XY 



1 

4 



/ -1 
-1 

-4 

4 

\-6 



-1 
-1 
i 

-4 
-6 



-1 1 



1 

2 
2 




-1 \ 
-1 



2 
2 







2 



/ 



;b.6) 



From the indices a, b, c, d which combine all but spin-indices we can see that we have 
done nothing but exchanged the two ^'s appearing in the four-fermion term. 

For the special case of only one fermionic species we find that the structure 
(ipOv4>) 2 + (iPOa^) 2 is invariant. Moreover, we can use Eqs. IJB.5J) . ()B.6j) to obtain 
the identity (actually this is exactly Eq. ()1.2j) ). 



(ijOvij) 2 - {^0 A ^f = -2[{^O s ^f - (^OpVOL (B.7) 



B.3. Fierz Identities 



95 



which allows us to transform (ipOyip) 2 — (iPOa4>) 2 completely into scalar and pseu- 
doscalar channels. 

If we have more than one fermion species, e.g. several flavors and/or colors, the 
Fierz transformations turns singlets into non-singlets and vice versa. This can be 
used to reduce the number of possible couplings as we do in Chap. 

Finally, let us mention two possible generalizations of the above. First, the same 
idea of permuting the ?/>'s can, of course, also be applied to higher order interactions 
like a 6- or 8-fermion interaction. Second, the four-fermion operators considered 
above are invariant under the discrete transformations C, P and T (charge conjuga- 
tion, parity and time reversal). However, as we know the weak interactions violate 
parity. Therefore, we might also want to consider interactions which are only invari- 
ant under CP and T. The set of parity violating operators can be obtained by mul- 
tiplying one Ox in the four-fermion operator by 7 5 . This yields (ipOx^iipOxl^ip)- 
Of course there is a set of equations similar to Eqs. (|B.5j) . (JB.6J) . 



Appendix C 



Infrared and Ultraviolet 
Regularization 

One of the central building blocks of the flow equation Eqs. 1)2.32)) . 1)2.34)) is the IR 
cutoff Rk- In this appendix we want to give explicit examples of the cutoff functions 
and define the threshold functions which appear in explicit calculations as a conse- 
quence of the trace over momentum space. Furthermore, we briefly comment on some 
technical points concerning the correspondence between UV and IR regularization. 



C.l Cutoff Functions 

To derive the flow equation ()2.34|) in Chap. EJ we added an additional term to the 
effective action providing for an IR regularization, 

AS k [<p] = 1 -j^ T (-p)R k (p)<f ) (p). (C.I) 

The inverse massless average (i.e. regularized) propagator P B for bosons and the 
corresponding squared quantity Pp for fermions are given by (cf. Sect. 12.1.2) Eqs. 

P B = q 2 + Z^ k R k (q) = q 2 (l+r B (q 2 )), (C.2) 
P F = q 2 (l + r F (q 2 )) 2 , 

where r B and rp reflect the presence of the IR cutoff. The dimensionless functions 
r B and r F depend only on y = q 2 /k 2 . 
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Expressed in terms of the r's the linear cutoff, Eqs. (J2HEJ), (|2.18jl ) is as follows, 

r B (y) = Q-lW-i,), (C.3) 
r F (y) = (-!--l)Q(l-y). 



Other typical examples are the sharp momentum cutoff, 



TF{y) 

and the popular exponential cutoff 



6(1 -y) 

i-e(i- y y 



r B {y) = r^" 1 ' ( as ) 

1 - exp(-y) 
V 1 - exp(-y) 



C.2 Threshold Functions 

The (super-)trace in the flow equations (|2.32jl . (|2.34|l includes a momentum space 
integral, reminiscent of a one-loop expression. Typically, the integral kernels are 
products of the IR regularized propagators and their derivatives. In most parts of 
this work, we use standard threshold functions as defined in |H^ . which we evaluate 
below explicitly for a finite UV cutoff A and for the linear cutoff (|C3|) . In Sect. 15.41 we 
need several additional threshold functions, not defined in the standard literature. 
To facilitate the notation we define a new, enlarged set of threshold functions in 
Sect. IC.2.3[ and label the threshold functions somewhat differently. 



C.2.1 Evaluation with Finite UV Cutoff 

To adapt the IR regulator to a finite UV cutoff one can modify the cutoff functions 
by a term which becomes infinite for all y > p - . For our purposes it is simpler to 
restrict the range of x to [0, A 2 ]. This has the same effect. In absence of mass terms 
the threshold functions can only depend on the ratio s = k 2 /A 2 . With 

~ q 2 d[Z^ k r B ] d 2 P F d[Z^ k r F ] d 
1 Z^ k dt 8P B Z^ k l + r F dt dP F [ ' ' 
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we find for bosons (x = q 2 



l{B),d/ \ 



1 

2 
2 
d 



k~ d 
1 - 



dxxi 1 d t \n(P B (x) + uk 2 



(C.7) 



2 _d 

+ ~S 2 

d 



d + 2 
1 - 



-0(1 -s) 



1+U 

(2 + d(l -s- 1 ))^ 
2(d + 2) 



1 + u 



Q(s-l) 



and for fermions 



{v,Vrl>,S) 



w 


2 


Jo 


2 




d 


1 - 

( 



dxxi 1 d t ln(P F (x) + uk 2 ) 



(C.8) 



2 _ d 

a 



1 - 



-6(1 



1 + u 
(d+1 — ds' 1 )^ 
d+1 



1 + u 



e( s -i; 



Higher threshold functions can be obtained simply by differentiating with respect 
to us: 

i d n+1 (u, v , s) = Vi s) . (c.9) 

7i + o„,o aw 

For a finite value of the UV cutoff A the threshold functions are explicitly s- and 
therefore fc-dependent. Taking A — > oo we have s = for any value of k. This renders 
the threshold functions /c-independent. 

In the pointlike truncations of Chaps. OH we neglect the anomalous dimensions 
r)<f>, r)ip and effectively only consider a fermionic cutoff since = 0. Moreover, the 
fermions are massless and we abbreviate for u = 

£(o,o, *) = #(*). 

This yields explicitly 



(CIO) 



4 n4 M 



1 



[e(i-s) + s - 2 e( s -i)]. 



To obtain the perturbative result from the fermionic RG equation we used 

A 2 



dtkn\ 



2,(F)A, 



(C.11) 



(C.12) 



As long as we keep the sharp momentum cutoff at q = A this integral is universal, 
i.e. it does not depend on the precise choice of the IR cutoff. Indeed the universality 
is necessary to reproduce perturbation theory for every choice of the IR cutoff. 



C.2. Threshold Functions 
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C.2.2 Cutoff Independence for Pointlike Truncations 

We have also used other cutoff functions R k different from the linear cutoff. Within 
the local interaction approximation we have found that the value of the critical 
coupling comes out independent of the choice of R k . The basic reason is that a 
multiplicative change of Z^' 4 due to the use of another threshold function can be 
compensated by a rescaling of k (cf. Eq. (J3.17)) ). The rescaling is simply multiplica- 
tive for s < 1, with a suitable generalization for s > 1. Critical values of the flow 
which are defined for k — > oo are not affected by the rescaling. Let us demonstrate 
this for X V - Writing Eq. ()3.17|) in the scale variable k we have 

d k J v ,k = -Av4 n \s)k(\ k + V, fc ) 2 . (C.13) 

Rescaling to 

k(k) = f dkkl[ F)A {s) (C.14) 
Jo 

we find 

d k K, k = -^M\ k + K,k) 2 - (c.i5) 

Due to the universality of Eq. (jC.12|) the domain for k, [0, oo], is now mapped 
to [0, 4r], giving the domain for k independent of the IR cutoff. Having obtained 
identical differential equations for every choice of IR cutoff without any rescaling of 
A establishes the above claim for the critical couplings. 

Note however, that this would not hold if we would start the integration of the 
flow equation at k — A. In this case the domain [0, A] for k is mapped into an 
interval for k that depends on the threshold function and therefore on R k . Actually, 
the R k dependence in this case is not very surprising because different IR cutoffs 
then correspond to different UV regularizations. Since our model is naively non- 
renormalizable results can depend on the choice of UV regularization (cf. Sect. lC.3"]) . 



C.2.3 Threshold Functions for Sect. ET4l 

Similar to Eq. (|C.7|) . (|C.8|) we define our modified threshold functions by integrals 
over x = q 2 . Threshold functions with derivatives of the fermion propagator are 
denoted by a Greek letter, all others by Latin letters. With 



(C.16) 
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and suppressing the arguments (u F , T) F , uji, u>2, 772) = {F,l,2) of the threshold 
functions and the argument (x) of the inverse propagators we write, 



(FBB),d _ 
n,mi ,1712 


2 


n+rai+m2- 


o{FBB),d _ 
r J n,mi ,m-2 


2 


n+mi+m 2 )- 


iAFBB),d _ 
u n,m\ ,1712 


-!*» 

2 


n+mi+m 2 )- 


(FBB),d _ 
in,mi,m,2 


2 


n+rai+1112)- 


(FBB),d _ 
^"n,m\ ,m2 


2 


n+mi+ms)- 


x(FBB),d _ 
u n,m\ ,m 2 


-±*> 

2 


n+mi+m 2 )- 


AFBB),d _ 
u 'n,mi ,vri2 


-!*> 

2 


n+rai+m2)- 


(FBB),d 
n,mi,iri2 


2 


n+mi+m 2 )- 


p {FBB),d _ 
°n,mi ,i7i2 


2 


n+mi+m 2 )- 



I>— m 2 
2 > 



D— m 2 
2 5 



3,-m 2 
2 5 



(C.17) 
(C.18) 
(C.19) 
(C.20) 
(C.21) 

; <9 t y dxx^FF-^PtP-^^P^ 2 , (C.22) 



3-m 2 
2 ) 



3>— m 2 
2 ) 



-(mi+1) p p -(m 2 +l) 



3>— m 2 
2 ' 



(C.23) 
(C.24) 
(C.25) 



If possible, i.e. if in the threshold function t no x-derivatives act on P 2 , respectively 
Pi and P2, we abbreviate 



(FBB),d _ ( FB ),d j .(FB),d _ .(F), 
L n,m,0 ~ b n,m aLLU - °n,0 ~ b n 



(C.26) 



For the linear cutoff Eq. (jC.3|) the integrals can be done explicitly as in Sect. 
IC.2.11 This is useful for numerical computations but not very enlightening. 

Some of the threshold functions defined above are straightforwardly related to 
the standard threshold functions used in the literature. Particularly noteworthy are 
the following relations, 



(C.27) 



The threshold functions defined above are not mutually independent. By partial 
integration one can obtain relations linking some of the threshold functions above. 



C.3. UV Regularization - ERGE Scheme 
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Let us content ourselves with two examples which are useful for a comparison with 
results in the literature, 

g (FBB),d( F i = — a (FBB),d-2 (F i 2) _ <B B) ,d / p , 2) _ ^2,(F B B) ,d / p 2 -.x 

Zilb lb lb 

(C.28) 

In the special case mi = m2 = 0, d = 2n = 4 this reduces to 



/3f ) ' 4 (F)=af ) ' 2 (F). (C.29) 



Moreover, 



+ 1. 2) + m 2 5™ 2 d (F, 2, 1), 

considering the same case as above and using relation (jC.29|) we find 

ef^(F) + xaf >' 2 (F) + y/^(F) = 3 7 f } ' 4 (F) = Smf^F), for (* + „) = 3. 

(C.31) 

Finally, let us come to a property we need when studying the pointlike limit for 
the bosons. For large boson masses the inverse propagator acts like uo 2 . This gives 
us the asymptotic properties 1 , 

«££2£ = - l -k 2( - n+m ^ m ^- d 8 t j dxxiF~ n {x)ujf m \x)uj 2 2m \x) (C.32) 
' -a^ d 



2mi 2m 2 



C™ = -h 2i - n+m ^ +m ^- d d t j dxxiF- n {x)P l {x)uj- 2(mi+l \x)uj- 2m ^ 



- b^ F) ' d 

2(mi+l) 2m 2 n ' 



^1 ^2 



and similar for all other threshold functions. 



C.3 UV Regularization - ERGE Scheme 
C.3.1 Effect of UV Regularization 

UV regularization can be implemented by a sharp cutoff in all integrals over mo- 
mentum space 2 . Yet, this is not the only possibility. Indeed, it is often not the most 

lr rhis argument and the corresponding properties are only valid at fixed anomalous dimension. 

2 When choosing a UV regularization one has to be careful to take one compatible with the 
symmetries of the theory in question. E.g. for gauge theories the sharp momentum cutoff violates 
gauge symmetry and is therefore not a suitable choice. 
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Figure C.l: In a UV regularized theory not all modes contribute completely. This 
plot schematically depicts "how much" each mode contributes. The thick line is 
for the sharp momentum cutoff. All modes with q 2 < A 2 are included completely. 
Other UV regularizations (dashed, thin dashed and thin solid line) typically not 
only include a small fraction of the high momentum modes, q 2 > A 2 , but in addition 
leave out a small fraction of the low momentum modes. 

practical regularization. E.g. in perturbation theory dimensional regularization is 
often much more convenient. But, even at fixed spacetime dimension there are other 
regularization methods. Prominent examples are the Pauli-Villars and Schwinger 
proper time regularization (for details s. e.g. |l20j ). More or less any modification 
(compatible with the symmetries) of the short distance behavior of the propagators 
which renders all Feynman diagrams finite, can be called a UV regularization. In 
general, the modes with q 2 > A 2 are not completely left out, only suppressed, as 
depicted in Fig. IC.ll 

From another point of view, a modification of the propagator can be implemented 
in the action. This gives then an "UV regularized classical action". The function 
describing the suppression of the UV modes (cf. Fig. IC. 1)1 in some way appears in this 
UV regularized action, e.g. the sharp cutoff limit can be implemented by multiplying 

2 

all terms in the action by appropriate 0(1 — -p)-functions. Hence, different UV 
regularizations usually correspond to different "classical actions". Therefore, it is 
no surprise, that different UV regularizations give different results. In particular, 
this is true for the critical coupling in our NJL model Eq. (jl.l)) . as one can see by 
comparing Tabs. IH.l) IH.2I with Tabs. 15.1) 15.2) calculated using a sharp UV cutoff at 
A and a UV regularization by the ERGE scheme (s. below) with the linear cutoff 
Eq. (jC.3j) . respectively. Only in renormalizable theories 3 we can remove the cutoff 
and obtain results independent of the specifics of the UV regularization. 



3 This includes critical behavior at second-order phase transitions (cf. Chap. 0, in particular 
critical exponents and other universal quantities. 
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C.3.2 ERGE Scheme 



In Sect. IC.2.11 we have evaluated the threshold functions for a theory which is UV 
regularized by a sharp cutoff in momentum space. It turned out that the threshold 
functions depended on the ratio s = j% in a rather complicated way. Constant 
threshold functions would be desirable, among other things to simplify numerical 
calculations. 

We already noted, that for A = oo the cutoff functions are indeed constant for 
all k, because s = for all finite k. Moreover, for an IR cutoff decreasing sufficiently 
fast in the UV all threshold functions are finite, even for A = oo. The flow equation 
(J2.32j) . ()2.34j) is UV finite. Now, let us remember that we have chosen the IR cutoff 
such that only the modes around q 2 = k 2 effectively contribute to the flow (Eq. 
(j2H3)). Starting the flow at k = A and integrating to k = we have included 
only modes with q 2 < A 2 . This is exactly what we expect for an UV regularization. 
More precisely, the UV regularization is now implemented in the finite choice for 
the initial conditions at k = ko = A. This is the so called ERGE scheme for the UV 
regularization. 

As flow equations are much simpler with constant threshold functions, this is 
the method of choice to implement UV regularization in ERGE. Nevertheless, this 
is bought at the prize that we cannot any longer compare non-universal quantities 
between different IR regularizations, as they automatically lead to different UV 
regularizations. 

Although it is usually not the simplest method, we can invoke UV regularization 
by the ERGE scheme also in the context of perturbation theory or SDE. This fol- 
lows along the lines indicated in Sect. 12.1.21 for perturbation theory. Typically any 
expression can be written in terms of inverse propagators P, internal momenta q we 
integrate over, and external momenta p we do not integrate over, 



A specific ERGE scheme is specified by the choice of the IR regulator R^. Replacing 
the inverse propagator P by the IR regularized inverse propagator P + we can 
calculate the contribution from each scale k, k~ 1 d t F{P + Rk, q,p)- Integrating over 
all scales from k = A to k = we obtain the UV regularized expression, 




(C.33) 




(C.34) 



Appendix D 

Flow Equations for Sect. 15.4 



In this appendix we list the flow equations for the effective action (jl.4j) generalized 
to include kinetic terms ()4.23|) . Moreover, we calculate the (momentum-dependent) 
field redefinitions necessary to keep a simple form of the action with Yukawa cou- 
plings constant in momentum space and no four-fermion interactions. 



D.l Flow Equations at Fixed Fields 

The first diagram that we have to evaluate is Fig. 13.11 giving the self-energy fi 2 (p). 
Momentum conservation implies that it depends only on one momentum variable p. 
Expanding for small values of p 2 and abbreviating the arguments of the threshold 
functions (cf. App. IC.2~3|) as F = (uf,t)f) we find, 

dt/4 = 8v A k 2 hla{ F)A (F), (D.l) 
d t /4 = 8v 4 k 2 h 2 v a { 2 F)A (F), dtn\ = 8v 4 k 2 h 2 A a { 2 F )A (F) 



and 



^ = -^h 2 A F)A (F), (D.2) 

9tZ V 16 2 (F),4/ |-a 9 t Z A 16 2 (F),A/ jp\ 



3 v ,z v n Z A 3 



In general, the vertex correction depicted in Fig. 14.21 depends on two momentum 
variables. As discussed in Sect. 15 .4121 we perform the evaluation for the configuration 
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(pi,P2,Ps) = (p, oP> tP)- Witn a = (v<nV<r), v = (^v,Vv), A = (uja,Va) this yields, 



(0) 



ft 



0g 



(0) 



—2v4hv 
2v±h,A 



hla?f>\F,V)-h\a^\F,A) 
hlaif^F, a) + 2/^a™ V, V) + 2^a£ B) ' 4 (F, A) 
fc^'V, -) - 2h 2 v a£ B) '\F, V) - 2h\aZ B) '\F, A) 



(D.3) 



and 



(2) 



01 



B {2) 
Ph A 



Av^hcr 



h 2 v(4 F i B) ' 2 (F, V) + 2P^(F, V) - riT^F, V) + e^{F, V)) (D.4) 
-(V^A) 

hl(3aZ B) > 2 (F, a) + 3/^ 4 (F, a) - 2 7 ™ 4 (F, a) + 2e[f^ a)) 
+2hU3aZ B) ' 2 (F, V) + 3{3Z B)A (F, V) - 2r$*>*{F, V) + 2%*>*(F, V)) 

- hl(3aZ B) ' 2 (F, a) + 3/3™ 4 (F, a) - 2 7 ™ 4 (F, a) + 2e^>\F, a)) 
+2^(3a^' 2 (F, V) + 3(3Z B)A (F, V) - 2 7 ™ 4 (F, V) + 2e^' 4 (F, V))) 
+(V^A) 



? (FB),4 



(FB),4, 



In order to compare the momentum dependence on s and t in the vicinity of 
(s,t) = (0,0) we have evaluated both momentum configurations discussed in Sect. 

For the configuration (pi, P2, P3, Pi) — \{PiPiPtP) with t = and s = p 2 we find 

-Sv^hlaZT^iF.a.V), (D.5) 
24t» 4 ftuft^og B 1 B, ' 4 (F,l/,'l), 
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v 4 h 2 a h 2 A 



P. 



(2) 



3 {2) 



2b^(F, a, A) - 2c™' 4 (F, a, A) + a, a) (D.6) 

+ei F 1 f ) '\F,a,A) + (a^A) 
9b^\F, V, A) - Wc^>\F, V, A) + 5d^(F, V, A) 
+ei F 1 f )A (F,V,A) + (V~A) 
44(96f 2 B) ' 4 (F, V) - 5c£ B) ' 4 (F, V) + 5e£ B) ' 4 (F, V)) + {V - A) 
(-3ig^(F, a) + cg^F, a) - e^(F, a) 



After the appropriate Fierz transformation the configuration {pi,P2,P3,P4) 
\(p-> ~Pi ~P,p) with s = and t = p 2 yields, 



0. 



(0) 



/* 4 a™ 4 (F, a) + 24/^4™>' 4 (F, V, A) - U^ag^F, V) (D.7) 



(0) 
Av 



-U 4 



~/i 4 a^ 4 (F, a) + tt^C^. ^> A ) - Mh 2 v h 2 A aZT l \F, V, A) 

4 (FB),4/ p T ,x „,4 (FS),4, 



^A A 



-t> 4 



,2,2 (FBB),4 



2 ,2 (FBB), 4 I 



-6hU^ BU (F, V) - 6h A a^(F, A) 

K(a^(F, a) - 4/^> 4 (F, a) + ^>\F, a) - e^ )A (F, a)) (D. 
-8h 2 v h 2 A (5aiZ B) ' 2 (F, V, A) + 4/3™' 4 (F, V, A) - V, A) 

+3eiT B)A (F, V, A)) - \{A - V) - \{V - A)] , 

, 4 ^(3ag^ a (F, a) + 3/^ 4 (F, a) - 2 7 ™ 4 (F, a) + 2 e £ B) ' 4 (F, a)) 



4 (FB),4, 



*2? 



r^4 



(2) 



12 



+4h 2 X(3aZT' 2 (F, a, V) + ^^(F, a, V) + <^-(F, a, V)) 
+4^^(3a£ B 1 B )' 2 (F, a, A) + 6/3™' 4 (F, a, A) - 5 72 ™>< 4 (F, A) 



(FBB),4, 



_(FBB),4. 
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+34™ ) ' 4 (F,a,A)) 
+8^^(9a£f } ' 2 (F, V, A) - V, A) + 16 7 (™ } ' 4 (F, V > A ) 

-AeiZ B) '\F, V, A)) + \{A - V) + \{V - A) 

?4(3a™ 2 (F, a) + S^'V, a) - 2 7 ™ 4 (F, a) + 2 e ™ 4 (F, a)) 

-AhlhU3a^(F, a, V) + 7 £^' 4 (F, *0 + e™ 4 (F, a, V)) 
-4/4/4(34™ } ' 2 (F, a, A) + a, A) - S^f^F, A) 

+3 e ™ 4 (F,a,A)) 
+8/£/£(9a<™ ] '\F, V, A) - 21^™ ) ' 4 (F, V, A) + 16 7 ^ B) ' 4 (F, V, A) 

-AeiZ B) '\F, V, A)) + \{A - V) + \{V - A) 



D.2 Field Redefinitions 

Before we start, let us remark that in this section we write down explicit factors of 
the wave function renormalization Z. To obtain the expressions in the renormalized 
couplings we simply have to set Z — 1. To keep the form of the effective action 
simple (more precisely to retain Yukawa couplings constant in momentum space 
and \{p) = 0) we allow for momentum-dependent field redefinitions, 

dtfa) = + 0(<7)<W?), (D.9) 

dt^(q) = (^(^^)^) (^ CT (g)+0*(g)^M<7), 

dtW(q) = { : h^){q)d t u Jv {q) + V^q)d t a v {q), 

d t A»(q) = ^ 1 ^){q)d t u A {q)+A^q)d t a A {q). 

Evaluating 

(D.10) 

we find 

d t K{q) = d t Xa(q) | -h a d t u> a (q), (D.ll) 

d t K(q) = d t K(q) \ + (fi 2 a + Z ' a q 2 )d t uJ < a (q) + h a d t a a (q), 

dtnl = dtfil | +2fild t a a (0), 

d t ZM) = d t Z a (q) | +2^ dta ° iq) ~^ (0) +2Z a d t a a (q), 
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d t \v(q) = d t X v (q) I +2h v d t u v (q), 

d t h v {q) = d t h v (q) | -(//y + Z v q 2 )d t u v (q) + h v d t a v {q), 



(D.12) 



dt/iy = d t /iy | +2jj 1 yd t a v (q), 



d t Z v (q) = d t Z v (q) | +2fi\ 



dtOLvjq) - d t a v (0) 
q 2 



+ 2Z v d t a v (q), 



and an analogous equation with (V — > A). As the expressions for the axial vector 
boson can always be obtained by this replacement, we write in the following only 
the expression for the vector boson. 



Imposing the condition 

determines the functions d t u>, 
d t uj a (q) = 



d t X(q) = 

d t \a(q) 



d t u v (q) 



d t X v (q) 



2h 



v 



The requirement for a constant Yukawa coupling reads, 



K{q) - MO) 



o, a 



h v {q)-h v (Q) 



q j \ q 

This fixes the functions d t a(q) up to a constant in q, 



0. 



(D.13) 



(D.14) 



(D.15) 



d t a a (q) - d t a ff (0) 



1 

K 



K{q) ~ M°) 



q 1 



(D.16) 



+ f4 df f K(q) ~ 



tin 



d t a v (q) - d t a v (0) 



1 

h v 



d, 



q z 



h v (q) ~ h v (Q) 



/4 a fx v (q)-x v (o) 



2h 



v 



q s 



7 

~rd t Xn-(q) 



+ ^ L d t Xv(q) 



2h 



v 



The remaining constant can be fixed by requiring that our fields are always renor- 
malized if they are so at the beginning, i.e. 



d t Z(0) = 0. 

Using our definition of the anomalous dimension, 77 = 2dt<y(°\ we find, 



(D.17) 



r)v = 2d t a { y ) 



d t Z a (0) I + 2£ 



d t Z v (0) I | 2/i 



Zfjhfj 



Zyhy 



MP 1 +^d t xl 2) I +^d t xf 1 
d t h^ I +S L d t x$ I +^d t x 



(D.18) 



(0) 



2h 



2h 



v 
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Inserting this into Eqs. (jD.lljl . (|D.12|) we obtain our final form for the flow equations, 

d t h 2 a = V X + 2h a d t h a \+2£d t \W\, (D.19) 
d t h 2 v = Vv h 2 v + 2h v d t h v \+^ v d t \^\, 



d t i4= VtT i4 + d tf il\, (D.20) 
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